Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,103,107 T→G 11.3% pseudogene (283/446 nt) wbbL ← pseudogene, lipopolysaccharide biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,103,1070TG11.3% 302.1 / 32.4 97pseudogene (283/446 nt)wbbLpseudogene, lipopolysaccharide biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base T (41/45);  new base G (5/6);  total (46/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.19e-01

Only 100 of 102 total aligned reads displayed.

AATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCTGGCA  >  NC_000913/2102967‑2103254
                                                                                                                                            |                                                                                                                                                   
aaTAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatgat                                                                                                                                               >  4:1537727/1‑147 (MQ=255)
aaTAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatgat                                                                                                                                               <  3:1537727/147‑1 (MQ=255)
aaTAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCAtgatga                                                                                                                                             >  3:1706008/1‑149 (MQ=255)
 aTAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCATgatgat                                                                                                                                            <  2:906121/149‑1 (MQ=255)
 aTAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCATgatgat                                                                                                                                            <  4:1706008/149‑1 (MQ=255)
    aaaTCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATc                                                                                                                                         <  1:289501/149‑1 (MQ=255)
     aaTCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCg                                                                                                                                        >  1:1521433/1‑149 (MQ=255)
        cAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGa                                                                                                                                     <  1:1204335/149‑1 (MQ=255)
        cAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGa                                                                                                                                     <  2:1772962/149‑1 (MQ=255)
          gAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCATgatgat                                                                                                                                            <  2:310780/140‑1 (MQ=255)
          gAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCATgatgat                                                                                                                                            >  1:310780/1‑140 (MQ=255)
               cacaGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCaaaa                                                                                                                              >  1:780578/1‑149 (MQ=255)
                acaGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCtt                                                                                                                                                     <  1:1156808/125‑1 (MQ=255)
                acaGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCtt                                                                                                                                                     >  2:1156808/1‑125 (MQ=255)
                acaGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCaaaaa                                                                                                                             <  2:114253/149‑1 (MQ=255)
                acaGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCAtgatga                                                                                                                                             >  2:676260/1‑133 (MQ=255)
                acaGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCAtgatga                                                                                                                                             <  1:676260/133‑1 (MQ=255)
                   ggAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACaa                                                                                                                          >  3:517621/1‑149 (MQ=255)
                    gAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGa                                                                                                                                     <  3:1471263/137‑1 (MQ=255)
                    gAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGa                                                                                                                                     >  4:1471263/1‑137 (MQ=255)
                    gAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACaaa                                                                                                                         <  2:780578/149‑1 (MQ=255)
                     aaaTTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACaaaa                                                                                                                        <  4:517621/149‑1 (MQ=255)
                          ttCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTaa                                                                                                                   <  2:1418319/149‑1 (MQ=255)
                             ttACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCtt                                                                                                                                                     <  1:1656397/112‑1 (MQ=255)
                             ttACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCtt                                                                                                                                                     >  2:1656397/1‑112 (MQ=255)
                             ttACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGa                                                                                                                                     <  1:1656689/128‑1 (MQ=255)
                             ttACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGAAAATCATCATGCTTCATGATGATATCGGGa                                                                                                                                     >  2:1656689/1‑128 (MQ=255)
                               aCGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAatcatc                                                                                                              >  1:184599/1‑149 (MQ=255)
                                       atCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCgg                                                                                                      <  2:184599/149‑1 (MQ=255)
                                        tCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCggg                                                                                                                                      >  3:1873240/1‑116 (MQ=255)
                                        tCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCggg                                                                                                                                      <  4:1873240/116‑1 (MQ=255)
                                          aTGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTc                                                                                                   >  2:1290021/1‑149 (MQ=255)
                                           tGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCt                                                                                                  <  2:1521433/149‑1 (MQ=255)
                                             tAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatgat                                                                                                                                               <  3:1592770/102‑1 (MQ=255)
                                             tAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatgat                                                                                                                                               >  4:1592770/1‑102 (MQ=255)
                                             tAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCaaa                                                                                                                               <  4:187304/118‑1 (MQ=255)
                                             tAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCaaa                                                                                                                               >  3:187304/1‑118 (MQ=255)
                                             tAAAGATTCCGCTTCATCTCGGAACAGGCATGATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCAtgatga                                                                                                                                             <  4:1552294/104‑1 (MQ=255)
                                             tAAAGATTCCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCAtgatga                                                                                                                                             >  3:1552294/1‑104 (MQ=255)
                                                 gATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCaaaaa                                                                                                                             <  4:1072112/116‑1 (MQ=255)
                                                 gATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCaaaaa                                                                                                                             >  3:1072112/1‑116 (MQ=255)
                                                   tttCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACaaaa                                                                                                                        >  2:1530585/1‑119 (MQ=255)
                                                   tttCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACaaaa                                                                                                                        <  1:1530585/119‑1 (MQ=255)
                                                       gCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGa                                                                                                                                     <  2:133619/102‑1 (MQ=255)
                                                       gCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGa                                                                                                                                     >  1:133619/1‑102 (MQ=255)
                                                         ttCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatga                                                                                                                                                >  3:1832478/1‑89 (MQ=255)
                                                         ttCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatga                                                                                                                                                <  4:1832478/89‑1 (MQ=255)
                                                           cATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGt                                                                                                    <  4:1339893/131‑1 (MQ=255)
                                                           cATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGt                                                                                                    >  3:1339893/1‑131 (MQ=255)
                                                            aTCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAatcatca                                                                                                             >  3:2396/1‑121 (MQ=255)
                                                            aTCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAatcatca                                                                                                             <  4:2396/121‑1 (MQ=255)
                                                             tctcGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGCTATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATa                                                                                >  3:289356/1‑149 (MQ=255)
                                                             tctcGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATc                                                                                                                                         >  2:208328/1‑92 (MQ=255)
                                                             tctcGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATc                                                                                                                                         <  1:208328/92‑1 (MQ=255)
                                                             tctcGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATa                                                                                <  1:1290021/149‑1 (MQ=255)
                                                               tcGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCaaaaa                                                                                                                             >  2:776478/1‑102 (MQ=255)
                                                               tcGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCaaaaa                                                                                                                             <  1:776478/102‑1 (MQ=255)
                                                                          aTAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCATgatgat                                                                                                                                            <  2:396895/76‑1 (MQ=39)
                                                                          aTAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCATgatgat                                                                                                                                            >  1:396895/1‑76 (MQ=38)
                                                                                   tAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATTCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCtt                                                                                      <  2:1345505/121‑1 (MQ=255)
                                                                                   tAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAATAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCtt                                                                                      >  1:1345505/1‑121 (MQ=255)
                                                                                   tAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAatcatc                                                                                                              >  3:566007/1‑97 (MQ=255)
                                                                                   tAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAatcatc                                                                                                              <  4:566007/97‑1 (MQ=255)
                                                                                     aaaGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTaa                                                                                                                   <  1:1160112/90‑1 (MQ=255)
                                                                                     aaaGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTaa                                                                                                                   >  2:1160112/1‑90 (MQ=255)
                                                                                     aaaGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAatcatcatc                                                                                                           <  2:1170947/98‑1 (MQ=255)
                                                                                     aaaGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAatcatcatc                                                                                                           >  1:1170947/1‑98 (MQ=255)
                                                                                     aaaGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATg                                                        >  4:69452/1‑149 (MQ=255)
                                                                                      aaGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGa                                                       <  3:69452/149‑1 (MQ=255)
                                                                                          aTAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAat                                                                                                                  >  4:746726/1‑86 (MQ=255)
                                                                                          aTAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAat                                                                                                                  <  3:746726/86‑1 (MQ=255)
                                                                                           tAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCg                                                                                                                                        >  2:657058/1‑63 (MQ=255)
                                                                                           tAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCg                                                                                                                                        <  1:657058/63‑1 (MQ=255)
                                                                                               gCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAatcatcatc                                                                                                           <  3:294928/88‑1 (MQ=255)
                                                                                               gCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAatcatcatc                                                                                                           >  4:294928/1‑88 (MQ=255)
                                                                                               gCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCg                                              >  1:1566634/1‑149 (MQ=255)
                                                                                                  tACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGtat                                           >  3:390740/1‑149 (MQ=255)
                                                                                                  tACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGtat                                           >  1:1393536/1‑149 (MQ=255)
                                                                                                     tttCGACATATTTAATATCTGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATTTGCGGGTCTatat                                                                                              >  1:1829183/1‑95 (MQ=255)
                                                                                                     tttCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTatat                                                                                              <  2:1829183/95‑1 (MQ=255)
                                                                                                     tttCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATaca                                        >  3:408756/1‑149 (MQ=255)
                                                                                                      ttCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATacac                                       <  4:390740/149‑1 (MQ=255)
                                                                                                      gctgACATATATTAAATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACaaaa                                                                                                                        <  3:946941/65‑1 (MQ=37)
                                                                                                          aCATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACa                                                                                  <  2:786404/102‑1 (MQ=255)
                                                                                                          aCATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACa                                                                                  >  1:786404/1‑102 (MQ=255)
                                                                                                                 tAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAata                            <  1:673636/149‑1 (MQ=255)
                                                                                                                  aaTATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTc                                                                                         >  1:589603/1‑87 (MQ=255)
                                                                                                                  aaTATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTc                                                                                         <  2:589603/87‑1 (MQ=255)
                                                                                                                  aaTATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACa                                                                                  <  4:148508/94‑1 (MQ=255)
                                                                                                                  aaTATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACa                                                                                  >  3:148508/1‑94 (MQ=255)
                                                                                                                   atatatGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAAtata                          <  4:408756/149‑1 (MQ=255)
                                                                                                                   atatatGGCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAAtata                          <  4:289356/149‑1 (MQ=255)
                                                                                                                    tatatGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAg                         <  2:1393536/149‑1 (MQ=255)
                                                                                                                      tatGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTc                       <  2:1566634/149‑1 (MQ=255)
                                                                                                                              cAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCgg                                                                                                      <  3:797506/62‑1 (MQ=255)
                                                                                                                              cAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCgg                                                                                                      >  4:797506/1‑62 (MQ=255)
                                                                                                                                         gCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCTgg    >  4:642738/1‑149 (MQ=255)
                                                                                                                                          cTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCTg     >  2:1330808/1‑147 (MQ=255)
                                                                                                                                          cTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCTg     <  1:1330808/147‑1 (MQ=255)
                                                                                                                                           ttCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCTGGCa  >  4:1050236/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
AATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCTGGCA  >  NC_000913/2102967‑2103254

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: