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)
             gTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGtata                                                                                 <  1:507027/100‑1 (MQ=255)
             gTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACCGCAGACTGGCAGGCCGCTTACGGTACGCTGGG‑‑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)
                 tCACGTAAGCCAGGGTGATACCTCGAATAATTGGTACACGGGCAGCTCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTa                                                                             <  1:871501/100‑1 (MQ=255)
                   aCGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACg                                                                           >  1:1941544/1‑100 (MQ=255)
                    cGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGAAGGCCGCTTACGGTACGAT‑‑GG‑‑GGGGCGGGTATAACTACGAc                                                                         >  1:319361/1‑100 (MQ=255)
                      tAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGAcc                                                                        >  2:1152134/1‑100 (MQ=255)
                         tcaattgTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGa                                                                    <  2:195713/94‑1 (MQ=255)
                         gCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTCCGACcgcg                                                                     >  2:2499698/1‑100 (MQ=255)
                          ccAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACcgcg                                                                     >  2:2146929/1‑99 (MQ=255)
                           cAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGAt                                                                   <  2:661837/100‑1 (MQ=255)
                           cAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGAt                                                                   >  2:1291977/1‑100 (MQ=255)
                            aGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATc                                                                  <  1:1201039/100‑1 (MQ=255)
                            aGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGAGATc                                                                  <  2:524031/100‑1 (MQ=255)
                             gggTGATACCTCGAATGATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCa                                                                 <  2:1972319/100‑1 (MQ=255)
                              ggTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCaa                                                                >  1:1825781/1‑100 (MQ=255)
                              ggTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCaa                                                                <  2:2480861/100‑1 (MQ=255)
                               gttaTACCTTGAATATTGGGTACTCGGGCAGCGCCACTGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAAc                                                               <  1:652990/97‑1 (MQ=255)
                               gTGATATCTCGAATAATGGGTACACGGGGAGCGCCACGGCAAACTGGCAGGCCGATTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAAc                                                               <  1:1280086/99‑1 (MQ=255)
                               gTGATACCTCGAATAATGGGTACACGGGCTGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAAc                                                               <  2:1592035/100‑1 (MQ=255)
                                 gATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACAt                                                             >  2:94328/1‑100 (MQ=255)
                                  tctccTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATg                                                            <  1:2532133/97‑1 (MQ=255)
                                  aTACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATg                                                            <  2:1634950/100‑1 (MQ=255)
                                  aTACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATg                                                            <  1:1294015/100‑1 (MQ=255)
                                  aTACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑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)
                                     ccTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGaga                                                         >  2:2532103/1‑98 (MQ=255)
                                     ccTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACg                                                         >  2:1555325/1‑100 (MQ=255)
                                     ccTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACg                                                         >  1:888002/1‑100 (MQ=255)
                                     ccTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGt                                                        <  2:2596661/100‑1 (MQ=255)
                                      cTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGATGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACg                                                         >  1:2393390/1‑99 (MQ=255)
                                       tCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGtt                                                       >  1:817152/1‑100 (MQ=255)
                                       tCGAATAATGAGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑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)
                                            taatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTg                                                  <  2:1740036/100‑1 (MQ=255)
                                             aatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTg                                                  >  1:1243652/1‑99 (MQ=255)
                                              atGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTgg                                                 >  1:2377362/1‑100 (MQ=255)
                                               tGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACTGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAg                                              <  2:1689093/100‑1 (MQ=255)
                                               tGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCa                                               >  1:522644/1‑100 (MQ=255)
                                               tGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCa                                               >  1:1999544/1‑100 (MQ=255)
                                               tGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCa                                               <  1:2240493/100‑1 (MQ=255)
                                                gggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAg                                              >  2:2127394/1‑100 (MQ=255)
                                                gggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGc                                             <  2:2502230/100‑1 (MQ=255)
                                                 ggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGc                                             >  1:1023279/1‑100 (MQ=255)
                                                   tACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTg                                           <  1:657423/100‑1 (MQ=255)
                                                    acacGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGt                                          >  1:382046/1‑100 (MQ=255)
                                                     cacGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTc                                         >  2:621389/1‑100 (MQ=255)
                                                     cacGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTc                                         <  1:679746/100‑1 (MQ=255)
                                                     cacGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GA‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCt                                        >  2:1471677/1‑100 (MQ=255)
                                                       cGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTg                                       >  1:1501305/1‑100 (MQ=255)
                                                       cGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTg                                       >  2:1823232/1‑100 (MQ=255)
                                                       cGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTAAGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTg                                       >  2:498345/1‑100 (MQ=255)
                                                        gggCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTgg                                      <  1:1212448/100‑1 (MQ=255)
                                                           cAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCgg                                   <  1:1056797/100‑1 (MQ=255)
                                                            aGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGgtg                                 <  1:288403/100‑1 (MQ=255)
                                                              cgcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGgtgt                                >  2:1803511/1‑100 (MQ=255)
                                                               gcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGtgtg                               >  2:1021301/1‑100 (MQ=255)
                                                                ccACGGCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGt                             >  2:1276419/1‑100 (MQ=255)
                                                                   cGGCAAACTGTCTGGCCGCTTACGGTACGCT‑‑GGGTGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTc                            <  1:505942/100‑1 (MQ=255)
                                                                     gCAAACTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCggg                         <  2:2399803/100‑1 (MQ=255)
                                                                     gCAAACTGGCAGGCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGc                        <  2:557868/100‑1 (MQ=255)
                                                                          cTGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGAGTGGTCGGGCATGa                    >  2:2441520/1‑100 (MQ=255)
                                                                           tGGCAGGCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGAATGaa                   <  2:754337/100‑1 (MQ=255)
                                                                               cttccgatctcGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATg               <  1:1209521/90‑1 (MQ=255)
                                                                                ggCCGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATgg              >  1:1239318/1‑100 (MQ=255)
                                                                                 gCCGCTTACGGTACGCTGGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATgg              <  1:1842771/100‑1 (MQ=255)
                                                                                 gCCGCTTACGGTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCa            <  1:2165622/100‑1 (MQ=255)
                                                                                  ccGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCa            <  2:1006928/100‑1 (MQ=255)
                                                                                  ccGCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCa            <  2:1409248/100‑1 (MQ=255)
                                                                                    gCTTACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATa          >  1:1577181/1‑100 (MQ=255)
                                                                                       tACGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCACGATCAACATGACGGTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACg       >  2:1184273/1‑100 (MQ=255)
                                                                                         cGGTACGCTGGG‑‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATgagatcggaag     >  2:1209509/1‑90 (MQ=255)
                                                                                           gTACGCT‑‑GG‑‑GGGGCGGGTATAACTACGACCGCGATCAAAATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGag  <  2:884950/100‑1 (MQ=255)
                                                                                            tACGCTGGG‑‑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: