Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 561,696:1 (G)6→7 100% coding (2000/2604 nt) sfmD → putative outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913561,6901.G69.0% 191.5 / 110.3 87coding (1994/2604 nt)sfmDputative outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base . (12/15);  new base G (32/27);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.02e-01

Only 100 of 107 total aligned reads displayed.

GGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAG  >  NC_000913/561593‑561783
                                                                                                  |                                                                                                
ggCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑gg                                                                                               <  2:18795/100‑1 (MQ=255)
ggCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑gg                                                                                               >  1:964993/1‑100 (MQ=255)
 gCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑gg‑‑g                                                                                            >  2:1754198/1‑100 (MQ=255)
 gCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑gg‑‑g                                                                                            >  1:557872/1‑100 (MQ=255)
 gCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑gg‑‑g                                                                                            >  2:735251/1‑100 (MQ=255)
 gCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑gg                                                                                               >  2:2099576/1‑99 (MQ=255)
  ccACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑gg‑‑gg                                                                                           <  1:2049302/100‑1 (MQ=255)
  ccACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑gg‑‑gg                                                                                           <  1:436310/100‑1 (MQ=255)
  ccACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑gg‑‑gg                                                                                           <  1:739924/100‑1 (MQ=255)
  ccACAACCTGAGTTATCACGTAAGCCAGGGAGATAACTCGAAAAATGGGTACACGGGCAGCGCCACGGCAACCTGGCAGGCCGCTTACGGACCGCT‑‑gg‑‑g                                                                                            >  2:811347/1‑99 (MQ=255)
   cacaACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑gg‑‑ggg                                                                                          >  1:1752359/1‑100 (MQ=255)
       aCCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCAGCTTACGGTACGCT‑‑gg‑‑gggggcg                                                                                      >  1:1564387/1‑97 (MQ=255)
        ccTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑gg‑‑gggggcgg                                                                                     >  1:1006938/1‑96 (MQ=255)
        ccTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑gg‑‑gggggcgg                                                                                     >  1:1985635/1‑96 (MQ=255)
        ccTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCggg                                                                                     <  1:572081/100‑1 (MQ=255)
         cTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑gg‑‑gggggcggg                                                                                    <  2:246781/100‑6 (MQ=255)
         cTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑gg‑‑gggggcggg                                                                                    <  2:1530081/100‑6 (MQ=255)
         cTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑gg‑‑gggggcggg                                                                                    <  2:900/100‑6 (MQ=255)
         cTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGt                                                                                    <  2:1320041/100‑1 (MQ=255)
         cTGAGTTATCACGTAAGCCAGAGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGTAGGCGGCTTACGGTACACT‑‑gg‑‑tggggc                                                                                       >  2:913278/1‑95 (MQ=255)
           gAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGtat                                                                                  >  2:2507258/1‑100 (MQ=255)
            aGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGtata                                                                                 >  2:2622955/1‑100 (MQ=255)
             gTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGtata                                                                                 <  1:507027/100‑1 (MQ=255)
             gTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACCGCAGACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGtata                                                                                 <  2:1252017/100‑1 (MQ=255)
               tATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACt                                                                              >  1:2348980/1‑100 (MQ=255)
               tATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACt                                                                              >  1:2502260/1‑100 (MQ=255)
                 tCACGTAAGCCAGGGTGATACCTCGAATAATTGGTACACGGGCAGCTCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTa                                                                             <  1:871501/100‑1 (MQ=255)
                   aCGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACg                                                                           >  1:1941544/1‑100 (MQ=255)
                    cGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGAAGGCCGCTTACGGTACGAT‑‑GG‑‑GGGGCGGGTATAACTACGAc                                                                         >  1:319361/1‑100 (MQ=255)
                      tAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGAcc                                                                        >  2:1152134/1‑100 (MQ=255)
                         tcaattgTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGa                                                                    <  2:195713/94‑1 (MQ=255)
                         gCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTCCGACcgcg                                                                     >  2:2499698/1‑100 (MQ=255)
                          ccAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACcgcg                                                                     >  2:2146929/1‑99 (MQ=255)
                           cAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGAt                                                                   <  2:661837/100‑1 (MQ=255)
                           cAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGAt                                                                   >  2:1291977/1‑100 (MQ=255)
                            aGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATc                                                                  <  1:1201039/100‑1 (MQ=255)
                            aGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGAGATc                                                                  <  2:524031/100‑1 (MQ=255)
                             gggTGATACCTCGAATGATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCa                                                                 <  2:1972319/100‑1 (MQ=255)
                              ggTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCaa                                                                >  1:1825781/1‑100 (MQ=255)
                              ggTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCaa                                                                <  2:2480861/100‑1 (MQ=255)
                               gttaTACCTTGAATATTGGGTACTCGGGCAGCGCCACTGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAAc                                                               <  1:652990/97‑1 (MQ=255)
                               gTGATATCTCGAATAATGGGTACACGGGGAGCGCCACGGCAAACTGGCAGGCCGATTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAAc                                                               <  1:1280086/99‑1 (MQ=255)
                               gTGATACCTCGAATAATGGGTACACGGGCTGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAAc                                                               <  2:1592035/100‑1 (MQ=255)
                                 gATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACAt                                                             >  2:94328/1‑100 (MQ=255)
                                  tctccTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATg                                                            <  1:2532133/97‑1 (MQ=255)
                                  aTACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATg                                                            <  2:1634950/100‑1 (MQ=255)
                                  aTACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATg                                                            <  1:1294015/100‑1 (MQ=255)
                                  aTACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATg                                                            <  2:2256336/100‑1 (MQ=255)
                                  aTACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGa                                                           <  1:205857/100‑1 (MQ=255)
                                    aCCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACg                                                         <  1:1628008/100‑1 (MQ=255)
                                    aCCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACg                                                         >  1:1545381/1‑100 (MQ=255)
                                     ccTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGaga                                                         >  2:2532103/1‑98 (MQ=255)
                                     ccTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACg                                                         >  2:1555325/1‑100 (MQ=255)
                                     ccTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACg                                                         >  1:888002/1‑100 (MQ=255)
                                     ccTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGt                                                        <  2:2596661/100‑1 (MQ=255)
                                      cTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGATG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACg                                                         >  1:2393390/1‑99 (MQ=255)
                                       tCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGtt                                                       >  1:817152/1‑100 (MQ=255)
                                       tCGAATAATGAGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGtt                                                       >  2:731183/1‑100 (MQ=255)
                                         gAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAAc                                                    <  2:1745329/100‑1 (MQ=255)
                                         gAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAAc                                                    >  2:1890839/1‑100 (MQ=255)
                                            taatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTg                                                  <  2:1740036/100‑1 (MQ=255)
                                             aatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTg                                                  >  1:1243652/1‑99 (MQ=255)
                                              atGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTgg                                                 >  1:2377362/1‑100 (MQ=255)
                                               tGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACTGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAg                                              <  2:1689093/100‑1 (MQ=255)
                                               tGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCa                                               >  1:522644/1‑100 (MQ=255)
                                               tGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCa                                               >  1:1999544/1‑100 (MQ=255)
                                               tGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCa                                               <  1:2240493/100‑1 (MQ=255)
                                                gggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAg                                              >  2:2127394/1‑100 (MQ=255)
                                                gggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGc                                             <  2:2502230/100‑1 (MQ=255)
                                                 ggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGc                                             >  1:1023279/1‑100 (MQ=255)
                                                   tACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTg                                           <  1:657423/100‑1 (MQ=255)
                                                    acacGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGt                                          >  1:382046/1‑100 (MQ=255)
                                                     cacGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTc                                         >  2:621389/1‑100 (MQ=255)
                                                     cacGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTc                                         <  1:679746/100‑1 (MQ=255)
                                                     cacGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GA‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCt                                        >  2:1471677/1‑100 (MQ=255)
                                                       cGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTg                                       >  1:1501305/1‑100 (MQ=255)
                                                       cGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTg                                       >  2:1823232/1‑100 (MQ=255)
                                                       cGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTAAGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTg                                       >  2:498345/1‑100 (MQ=255)
                                                        gggCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTgg                                      <  1:1212448/100‑1 (MQ=255)
                                                           cAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCgg                                   <  1:1056797/100‑1 (MQ=255)
                                                            aGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGgtg                                 <  1:288403/100‑1 (MQ=255)
                                                              cgcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGgtgt                                >  2:1803511/1‑100 (MQ=255)
                                                               gcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGtgtg                               >  2:1021301/1‑100 (MQ=255)
                                                                ccACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGt                             >  2:1276419/1‑100 (MQ=255)
                                                                   cGGCAAACTGTCTGGCCGCTTACGGTACGCT‑‑GGGTGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTc                            <  1:505942/100‑1 (MQ=255)
                                                                     gCAAACTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCggg                         <  2:2399803/100‑1 (MQ=255)
                                                                     gCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGc                        <  2:557868/100‑1 (MQ=255)
                                                                          cTGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGAGTGGTCGGGCATGa                    >  2:2441520/1‑100 (MQ=255)
                                                                           tGGCAGGCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGAATGaa                   <  2:754337/100‑1 (MQ=255)
                                                                               cttccgatctcGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATg               <  1:1209521/90‑1 (MQ=255)
                                                                                ggCCGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATgg              >  1:1239318/1‑100 (MQ=255)
                                                                                 gCCGCTTACGGTACGCTGGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATgg              <  1:1842771/100‑1 (MQ=255)
                                                                                 gCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCa            <  1:2165622/100‑1 (MQ=255)
                                                                                  ccGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCa            <  2:1006928/100‑1 (MQ=255)
                                                                                  ccGCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCa            <  2:1409248/100‑1 (MQ=255)
                                                                                    gCTTACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATa          >  1:1577181/1‑100 (MQ=255)
                                                                                       tACGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCACGATCAACATGACGGTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACg       >  2:1184273/1‑100 (MQ=255)
                                                                                         cGGTACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATgagatcggaag     >  2:1209509/1‑90 (MQ=255)
                                                                                           gTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAAAATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGag  <  2:884950/100‑1 (MQ=255)
                                                                                            tACGCTG‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGag  <  2:61887/100‑1 (MQ=255)
                                                                                                  |                                                                                                
GGCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAG  >  NC_000913/561593‑561783

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑