Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_007793 828,711 C→T T925I (ACA→ATA)  SAUSA300_RS04005 → ABC‑ATPase UvrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_007793828,7110CT98.8% 224.3 / ‑5.9 85T925I (ACA→ATA) SAUSA300_RS04005ABC‑ATPase UvrA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (78/4);  total (81/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

TCAAAACAGCAGACTATATTATAGACTTAGGTCCTGAAGGTGGTAGTGGCGGTGGTACTATTGTTGCGACT‑GGCACACCCGAAGATATTGCTCAG‑ACAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGATAAATAAGATTAAAAGAAGTGAAGGATG  >  NC_007793/828615‑828805
                                                                                                  |                                                                                               
acAAAACAGCAGACTATATTATAGACTTAGGTCCTGAAGGTGGTAGTGGCGGTGGTGCCTTTTTTTCGGCT‑GGGCAACCCCAAGATATTTCTCCG‑ATaaa                                                                                              >  1:1470614/2‑100 (MQ=255)
  aaaaCAGCAGACTATATTATAGACTTAGGTCCTGAAGGTGGTAGTGGCGGTGGGTGTAATTTTTCGGCC‑GGCACCCCCCCAGATATTTCTCCG‑GCaaaat                                                                                            >  1:6322521/1‑98 (MQ=255)
   caaCAGCAGACTATATTATAGACTTAGGTCCTGAAGGTGGTAGTGGCGGTGGGACTATTTTTTTGACT‑GGGACACCCCCAGAAATTGTTCCG‑ATAAAgac                                                                                           >  1:182438/2‑98 (MQ=255)
    aaCAGCAGACTATATTATAGACTTAGGTCCTGAAGGTGGTAGTGGCGGTGGGTCTATTTTTTGGACC‑GGCCCACCCCCAAATATTTCTCAG‑ATAAAAtca                                                                                          >  1:3012124/1‑100 (MQ=255)
    aaCAGCAGACTATATTATAGACTTAGGTCCTGAAGGTGGTAGTGGCGGTGGGGCTCTTTCTTCGTCT‑TGGACACCCCAAGATATTTGTCAG‑AGAAAgaaa                                                                                          >  1:5403297/1‑97 (MQ=255)
    aaCAGCAGACTATATTATAGACTTAGGTCCTGAAGGTGGTAGTGGCGGTGGGACTTTTTTTTGGACT‑GGGTCACCCCCCGACATTGCTCCG‑AAAAAGtca                                                                                          >  1:4373159/1‑100 (MQ=255)
    aaCAGCAGACTATATTATAGACTTAGGTCCTGAAGGTGGTAGTGGCGGAGGTACTTTTTTTTTGACT‑TGCCCCCCCCAAGATATTGCTTAG‑ATaaaataa                                                                                          >  1:6104018/1‑96 (MQ=255)
     aCAGCAGACTATATTATAGACTTAGGTCCTGAAGGTGGTAGTGGCGGTGGTAGCTTTTTTTTGACC‑TGCCCACCCCGAGAAATTTTTCAG‑ACAAAAtcat                                                                                         >  1:4710438/1‑100 (MQ=255)
      cagcagACTATATTATAGACTTAGGTCCTGAAGGTGGTAGTGGCGGTGGGAGTTTTGTTTTGGCT‑GGGGCACCCCCAGAAATTTCCCAG‑ATAAAgatatt                                                                                        >  1:3204677/1‑95 (MQ=255)
      cagcagACTATATTATAGACTTAGGTCCTGAAGGTGGTAGTGGCGGCGGGACTATTGTTTCGGCT‑GGCCCACCCGCAGATATTGCTCCG‑ATAAAGtcatc                                                                                        >  1:2566013/1‑100 (MQ=255)
       tgcagACTATATTATAGACTTAGGTCCTGAAGGTGGTAGTGGCGGTGGGACTATTTTTTGGGCT‑GGCACACCCCAAGAAATTGTTCCG‑ATAAAAacatca                                                                                       >  1:5600872/2‑100 (MQ=255)
       cgcagACTATATTATAGACTTAGGTCCTGAAGGTGGTAGTGGCGGTGGGACTATTTTTTTGACT‑GGGAAACCCGCAGAAATTTTTTAG‑ATAAAAacatca                                                                                       >  1:701857/2‑100 (MQ=255)
       agcagACTATATTATAGACTTAGGTCCTGAAGGTGGTAGTGGCGGTGGGACTATTTTTGCGGCT‑GGGCCACCCCCAGAGAATGCTCCG‑ATaaaaacaaca                                                                                       >  1:4307364/1‑93 (MQ=255)
        tcagACTATATTATAGACTTAGGTCCTGAAGGTGGTAGTGGCGGTGGGGCTATTTTTGGGACT‑GGGACACCCGAAGAAATTGCTCAG‑ATAAAAAcatcat                                                                                      >  1:3629206/2‑100 (MQ=255)
         cagACTATATTATAGACTTAGGTCCTGAAGGTGGTAGTGGCGGTGGGACTATTTTTTCGGCT‑GGGACACCCCAAGAAAATGTTCAG‑ATaaaaaaataata                                                                                     >  1:1325406/1‑91 (MQ=255)
          agACTATATTATAGACTTAGGTCCTGAAGGTGGTAGTGGCGGTGGGGCTATTTTTTTGGCT‑GGCACACCCCAAGATATTTCTCTG‑ATTAAAACATCatat                                                                                    >  1:400475/1‑100 (MQ=255)
            gcTATATTATAGACTTAGGTCCTGAAGGTGGTAGTGGCGGAGGTGCTATTTTTTTTGCT‑GGGTCACCCCAAGATGATGCTCTG‑ATAAAAtccatatatta                                                                                  >  1:1742173/2‑91 (MQ=255)
               atatTATAGACTTAGGTCCTGAAGGTGGTAGTGGCGGGGGGGATATTTTTTTGGCC‑GGCCCACCCCAAGATATTGTTCCG‑ATAAAAGAATCATAGTC‑Agg                                                                              >  1:5015347/1‑100 (MQ=255)
                   tataGACTTAGGTCCTGAAGGTGGTAGTGGCGGCGGGACTATTTTTGGGACT‑GGCCCACCCCCAGATATTGCTCAGAAAAAAGTC‑TCATATTC‑AAGAAAg                                                                          >  1:3499297/1‑100 (MQ=255)
                   aataGACTTAGGTCCTGAAGGTGGTAGTGGCGGCGGGACTAATGTTGGGCGC‑GGGGCAACCCAAGAAAATGTTCTG‑AAAAAAGCCTCATATAC‑AGGGAAg                                                                          >  1:462524/2‑100 (MQ=255)
                    ctaGACTTAGGTCCTGAAGGTGGTAGTGGCGGAGGTCCCTTTGTTGGGCCC‑CGCCCACCCCCAGAAAATGCCCAG‑ATAAAATCAATATATAA‑AGGGAAga                                                                         >  1:1795571/2‑99 (MQ=255)
                      aGACTTAGGTCCTGAAGGTGGTAGTGGCGGTGGGGCCATTTTTTGGACT‑CGCACCCCCCCAGATATTTGTCAG‑ATAAAGTCAACATATAA‑AAGaaaagat                                                                       >  1:1004466/1‑96 (MQ=255)
                         cTTAGGTCCTGAAGGTGGTAGTGGCGGTTGTGATATTTTTTTGACTAGGCCCACCC‑CAGATATTGTTCCG‑ATAAAAGTATCATATAC‑AAGGACATATTTa                                                                    >  1:403639/1‑100 (MQ=255)
                            aGGTCCTGAAGGTGGTAGTGGCGGTGGGTCTATTTTTTGGGCC‑GGGGCACCCCCAGAGATTTCTCCG‑ATAAAAAACTCACATAC‑AAGAAAGAAATTaaaa                                                                 >  1:1913756/1‑100 (MQ=255)
                            aGGTCCTGAAGGTGGTAGTGGCGGAGGGACTATTTTTGCGGCC‑GGCACCCCCGCAGAGAATGGTTCG‑GAAAAGACATCATTAAC‑AAGAAAGAAtttttaa                                                                 >  1:292930/1‑96 (MQ=255)
                             ggTCCTGAAGGTGGTAGTGGCGGGGGGGCTATTTTTTTGGCT‑TGCACAACCCGAGAAATTGCTCAG‑ATAAAAGCATTATATAC‑AGGAAAGAATTTaaaaa                                                                >  1:859992/1‑99 (MQ=255)
                                ccTGAAGGTGGTAGTGGCGGGGGGGATAATTTTTTGACT‑GGGACACCCCGAGAAAATTCTCAG‑ATAAAATTATCATATAC‑AAGAAAATAATTAaaaaaaa                                                             >  1:2401911/1‑99 (MQ=255)
                                  tGAAGGTGGTAGTGGCGGTGGGGCTATTTTTTTGACT‑GGGACCCCCCAAGAAAATTCTCAG‑GTAAAAACATCATATACAAGAAAAGAATTTaaaaaagga                                                            >  1:5089229/1‑95 (MQ=255)
                                    gaGGTGGTAGTGGCGGTGGGGCTATTGTTTGGGCT‑GGCGCACCCCCAGAAATTTCTCCG‑ATAAAAGCATCATCTAC‑AAGAAAAAATTTTAaaaaagaaat                                                         >  1:330361/2‑96 (MQ=255)
                                            cttGGCGGTGGTACTATTGTTGCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACgag                                                 >  1:122692/3‑100 (MQ=255)
                                             ttGGCGGTGGTACTATTGTTGCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACgaga                                                >  1:386925/2‑100 (MQ=255)
                                             ctGGCGGTGGTACTATTGTTGCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACgaga                                                >  1:1029464/2‑100 (MQ=255)
                                              tGGCGGTGGTACTATTGTTGCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGAt                                               >  1:1744992/1‑100 (MQ=255)
                                               ggCGGGGGGACTAATGGTGCGACT‑GGCACACCCGAAGATATTGCACAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATa                                              >  1:5629874/1‑100 (MQ=255)
                                                 cGGTGGTACTATTGTTGCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGGAAGTATTTAAAAGAAGTACTTGAACGAGATaaa                                            >  1:3212467/1‑100 (MQ=255)
                                                   gtggtACTATTGTTGCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACa                                          >  1:6204549/1‑100 (MQ=255)
                                                   atggtACTATTGTTGCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACa                                          >  1:3582041/2‑100 (MQ=255)
                                                    tggtACTATTGTTGCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAGGAAGTACTTGAACGAGATAAACaa                                         >  1:5622545/1‑100 (MQ=255)
                                                    tggtACTATTGTTGCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACaa                                         >  1:723552/1‑100 (MQ=255)
                                                    tggtACTATTGTTGCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACaa                                         >  1:1350729/1‑100 (MQ=255)
                                                    tggtACTATTGTTGCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACaa                                         >  1:4121309/1‑100 (MQ=255)
                                                    tggtACTATTGTTGCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTAAACGAGAtcggaag                                         >  1:4894822/1‑94 (MQ=255)
                                                    tggtACTATTGTTGCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTAAACGAGAtcggaag                                         >  1:328310/1‑94 (MQ=255)
                                                    tggtACTATTGTTGCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTCGAACGAGAtcggaag                                         >  1:315836/1‑94 (MQ=255)
                                                     ggtACTATTGTTGCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGGAAGTATTTAAAAGAAGTACTTGAACGAGATAAACaaa                                        >  1:6584644/1‑100 (MQ=255)
                                                     ggtACTATTGTTGCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACaaa                                        >  1:3412049/1‑100 (MQ=255)
                                                     ggtACTATTGTTGCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACaaa                                        >  1:1098278/1‑100 (MQ=255)
                                                         cccttgttgCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAAtgc                                    >  1:5025560/4‑98 (MQ=255)
                                                         cTATTGTTGCGACT‑GGCACACTCGAAGATATTGCCCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATAc                                    >  1:3373196/1‑100 (MQ=255)
                                                           ctcttccgatcT‑GGCATACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTg                                  <  1:269535/90‑1 (MQ=255)
                                                           aTTGTTGCGACT‑GGCACACCCGAGGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTg                                  >  1:3483295/1‑100 (MQ=255)
                                                           aTTGTTGCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTg                                  >  1:99674/1‑100 (MQ=255)
                                                           aTTGTTGCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTg                                  >  1:3049745/1‑100 (MQ=255)
                                                            ttgttgTGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGa                                 >  1:2687570/1‑100 (MQ=255)
                                                            tcttccgatcT‑GGTACGCCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAATATTTAAAAGAAGTACTTGAACGAGATAAACAAAATGCTGa                                 <  1:2548716/91‑1 (MQ=255)
                                                            tcttccgatcT‑GGCATACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACCGa                                 <  1:1363613/91‑1 (MQ=255)
                                                             tgttgCGACT‑GGCACACCCGAAGATATTGCTTAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGaa                                >  1:4228933/1‑100 (MQ=255)
                                                              gttgCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAACACTGAAg                               >  1:5566104/1‑100 (MQ=255)
                                                               ttgCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGa                              >  1:510258/1‑100 (MQ=255)
                                                                 gCGACT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACGAAATACAGAAGatc                            >  1:4849524/1‑99 (MQ=255)
                                                                     cT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGatcggaa                        >  1:3852777/1‑95 (MQ=255)
                                                                     cT‑GGCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAAGACTGAAGatcggaa                        >  1:2850224/1‑95 (MQ=255)
                                                                        ggCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAACGAAGTACTTGAACGAGATAAACAAAATACTGAAGACAAATAAg                      >  1:2263163/1‑100 (MQ=255)
                                                                        ggCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATGCTGAAGatcggaaga                      >  1:1388527/1‑93 (MQ=255)
                                                                        ggCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGaaagatcggaag                      >  1:5217650/1‑90 (MQ=255)
                                                                        ggCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGGAGatcggaaga                      >  1:1172885/1‑93 (MQ=255)
                                                                        ggCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGatcggaaga                      >  1:4126807/1‑93 (MQ=255)
                                                                        ggCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGatcggaaga                      >  1:1852196/1‑93 (MQ=255)
                                                                        ggCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGatcggaaga                      >  1:717823/1‑93 (MQ=255)
                                                                        ggCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGatcggaaga                      >  1:5037326/1‑93 (MQ=255)
                                                                        ggCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGaagatcgga                      >  1:5069440/1‑92 (MQ=255)
                                                                        ggCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGATAAATAAg                      >  1:5259815/1‑100 (MQ=255)
                                                                        ggCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGATAAATAAg                      >  1:2059519/1‑100 (MQ=255)
                                                                        ggCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACAGAAGatcggaaga                      >  1:1254229/1‑93 (MQ=255)
                                                                        ggCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAACACTAAAGatcggaaga                      >  1:2399132/1‑93 (MQ=255)
                                                                        ggCACACCCGAAGATATTGCTCAG‑ATAAAGCCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACAGAAGatcggaaga                      >  1:3188587/1‑93 (MQ=255)
                                                                        ggCACACCCGAAGATATTGCTCAG‑ATAAAGCCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAAGACTGAAGatcggaaga                      >  1:3073371/1‑93 (MQ=255)
                                                                        cacacaCCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGCACTTGAACGAGATAAACAAAACACTGAAGatcggaaga                      >  1:3093732/3‑93 (MQ=255)
                                                                         gCACACCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAACACTGAAGatcggaagag                     >  1:3762043/1‑92 (MQ=255)
                                                                         ccacaCCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAACACTGAAGatcggaagag                     >  1:5776677/2‑92 (MQ=255)
                                                                         atgcacCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGATAAATAAGa                     >  1:3184052/4‑100 (MQ=255)
                                                                          tacacCCGAAGATATTGCTCAG‑TTAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGCAGatcggaagagc                    >  1:2981251/2‑91 (MQ=255)
                                                                          cacaCCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATGCTGAAGatcggaagagc                    >  1:2997973/1‑91 (MQ=255)
                                                                          cacaCCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGatcggaagagc                    >  1:3196716/1‑91 (MQ=255)
                                                                          cacaCCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGatcggaagagc                    >  1:1191496/1‑91 (MQ=255)
                                                                          cacaCCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAAAACTGAAGatcggaagagc                    >  1:3274924/1‑91 (MQ=255)
                                                                          cacaCCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAAAACTGAAGatcggaagagc                    >  1:4658181/1‑91 (MQ=255)
                                                                           ccacCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGATAAATAAGAtt                   >  1:5288664/2‑100 (MQ=255)
                                                                           acacCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGatcggaagagca                   >  1:4637816/1‑90 (MQ=255)
                                                                           acacCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGATAAATAAGAtt                   >  1:4018146/1‑100 (MQ=255)
                                                                           acacCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAAGACTGAAGatcggaagagca                   >  1:341297/1‑90 (MQ=255)
                                                                           acacCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAACACTGAAGATAAATAAGAtt                   >  1:1822999/1‑100 (MQ=255)
                                                                             acCCGAAGATATTGCTCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGATAAATAAGAttga                 >  1:3873534/1‑98 (MQ=255)
                                                                                            tCAG‑ATAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGATAAATAAGATTAAAAGAAGCGagtgcat  <  1:1146787/100‑7 (MQ=255)
                                                                                                  |                                                                                               
TCAAAACAGCAGACTATATTATAGACTTAGGTCCTGAAGGTGGTAGTGGCGGTGGTACTATTGTTGCGACT‑GGCACACCCGAAGATATTGCTCAG‑ACAAAGTCATCATATAC‑AGGAAAGTATTTAAAAGAAGTACTTGAACGAGATAAACAAAATACTGAAGATAAATAAGATTAAAAGAAGTGAAGGATG  >  NC_007793/828615‑828805

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑