Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,370,708 G→C 9.9% P110R (CCG→CGG) ‡ nanE ← putative N‑acetylmannosamine‑6‑phosphate epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,370,7080GC 9.9% 204.2 / 11.4 68P110R (CCG→CGG) ‡nanEputative N‑acetylmannosamine‑6‑phosphate epimerase
Reads supporting (aligned to +/- strand):  ref base G (31/30);  new base C (4/3);  total (35/33)
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.01e-01

GCGTGGTATAGCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTG  >  NC_000913/3370559‑3370857
                                                                                                                                                     |                                                                                                                                                     
gCGTGGTATAGCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGAcgcg                                                                                                                                                       <  1:530807/150‑1 (MQ=255)
gCGTGGTATAGCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGAcgcg                                                                                                                                                       >  1:312413/1‑150 (MQ=255)
gCGTGGTATAGCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGAcgcg                                                                                                                                                       >  2:1892800/1‑150 (MQ=255)
  gTGGTATAGCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACgcggg                                                                                                                                                     <  2:444264/150‑1 (MQ=255)
          gCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGcggtcggt                                                                                                                                             >  2:1034375/1‑150 (MQ=255)
            cAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGc                                                                                                                                           >  1:979915/1‑150 (MQ=255)
             aGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGACCGCCCGCGTCCGGTGcc                                                                                                                                          <  1:1014280/150‑1 (MQ=255)
              gAAAGCGTAGTGCCAATAATTTCGGCTCCCAGATTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCg                                                                                                                                         >  2:452442/1‑150 (MQ=255)
                aaGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTc                                                                                                                                       >  1:111771/1‑150 (MQ=255)
                aaGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGACCGCCCGCGTCCGGTGcctgt                                                                                                                                       <  2:110165/150‑4 (MQ=255)
                 aGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCa                                                                                                                                      <  1:288367/150‑1 (MQ=255)
                      aGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGc                                                                                                                                 >  2:306687/1‑150 (MQ=255)
                        tGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGa                                                                                                                               <  1:905856/150‑1 (MQ=255)
                           cAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATaa                                                                                                                            >  1:188261/1‑150 (MQ=255)
                           cAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATaa                                                                                                                            <  1:316390/150‑1 (MQ=255)
                           cAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATaa                                                                                                                            >  2:314349/1‑150 (MQ=255)
                                atTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGAGATGGTGAATACGTGCCAGCAGCCTTTCAACACGCACCGGACGCCGCCCGTCGGTGCCGTCAATGGCGATAATGTcc                                                                                                                       >  2:1663460/1‑150 (MQ=255)
                                 tttCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCcg                                                                                                                      <  1:1781691/150‑1 (MQ=255)
                                 tttCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCcg                                                                                                                      <  1:933978/150‑1 (MQ=255)
                                     ggCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGcgcc                                                                                                                  <  1:1034375/150‑1 (MQ=255)
                                      gCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGcgccc                                                                                                                 >  2:1962623/1‑150 (MQ=255)
                                         cccAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGcgcccgcc                                                                                                              >  2:1610138/1‑150 (MQ=255)
                                            aGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTgc                                                                                                           >  1:53814/1‑150 (MQ=255)
                                             gCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTgcg                                                                                                          <  1:1899300/150‑1 (MQ=255)
                                                tttGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCa                                                                                                       >  1:1062666/1‑150 (MQ=255)
                                                   ggCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAgcg                                                                                                    <  2:117096/150‑1 (MQ=255)
                                                    gCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATATGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAgcgc                                                                                                   >  2:1962661/1‑150 (MQ=255)
                                                    gCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAgcgc                                                                                                   >  2:780635/1‑150 (MQ=255)
                                                     cATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCa                                                                                                  >  2:1171503/1‑150 (MQ=255)
                                                       tGCCAGGCCGTCTTCCGGAGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGACAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCAta                                                                                                >  2:1230177/1‑149 (MQ=255)
                                                           aGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACa                                                                                            <  2:187537/150‑1 (MQ=255)
                                                            ggCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACAt                                                                                           >  2:1169883/1‑150 (MQ=255)
                                                                 tCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCa                                                                                      <  1:1134475/150‑1 (MQ=255)
                                                                 tCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGACAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCa                                                                                      >  2:1210358/1‑150 (MQ=255)
                                                                  cTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCaa                                                                                     <  1:1892800/150‑1 (MQ=255)
                                                                  cTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCaa                                                                                     >  1:1873062/1‑150 (MQ=255)
                                                                     ccGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAatat                                                                                  >  1:1500492/1‑150 (MQ=255)
                                                                               gCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGAt                                                                        <  1:306687/150‑1 (MQ=255)
                                                                                  gTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCg                                                                     >  2:594604/1‑150 (MQ=255)
                                                                                      gTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTAcc                                                                 >  2:1445445/1‑150 (MQ=255)
                                                                                        cATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCgg                                                               >  2:627185/1‑150 (MQ=255)
                                                                                          tCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGgag                                                             <  1:314349/150‑1 (MQ=255)
                                                                                            gCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGaa                                                           <  1:594604/150‑1 (MQ=255)
                                                                                               aGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAAtcc                                                        <  2:1500492/150‑1 (MQ=255)
                                                                                                 tAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAAtcctc                                                      >  1:1735237/1‑150 (MQ=255)
                                                                                                   aaCCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCa                                                    <  2:53814/150‑1 (MQ=255)
                                                                                                            ggTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGtt                                           >  1:654870/1‑150 (MQ=255)
                                                                                                            ggTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGtt                                           <  1:1610138/150‑1 (MQ=255)
                                                                                                            ggTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGtt                                           >  1:1027138/1‑150 (MQ=255)
                                                                                                              tGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTc                                         <  1:1663460/150‑1 (MQ=255)
                                                                                                                aaTACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTcac                                       <  2:931283/150‑1 (MQ=255)
                                                                                                                     gTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCAAAAtta                                  >  1:984953/1‑149 (MQ=255)
                                                                                                                      tGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTcc                                 <  1:1171503/150‑1 (MQ=255)
                                                                                                                      tGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTcc                                 <  1:1947372/150‑1 (MQ=255)
                                                                                                                       gCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCa                                <  1:1582367/150‑1 (MQ=255)
                                                                                                                        ccAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCaa                               <  2:312413/150‑1 (MQ=255)
                                                                                                                          agcagcGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCTCCAGAGCATCAACATCTACAATATAGACAGTGATGCGTACCGGAGAAACCACCAGATCCCGTTTCACAAATACaata                             >  2:1683314/1‑150 (MQ=255)
                                                                                                                             agcGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATc                          <  1:1664401/150‑1 (MQ=255)
                                                                                                                              gcGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCg                         >  2:1392234/1‑150 (MQ=255)
                                                                                                                                gTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGc                       >  2:1420441/1‑150 (MQ=255)
                                                                                                                                   tCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACg                    <  1:344788/150‑1 (MQ=255)
                                                                                                                                   attgacggCACCGACCGCCCGCGTCCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACg                    >  1:1448050/7‑150 (MQ=255)
                                                                                                                                       cAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTca                <  1:1230177/150‑1 (MQ=255)
                                                                                                                                       acggCACCGACCGCCCGCGTCCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCa                                <  2:1969570/132‑1 (MQ=255)
                                                                                                                                       acggCACCGACCGCCCGCGTCCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCa                                >  1:1969570/3‑134 (MQ=255)
                                                                                                                                       acggCACCGACCGCCCGCGTCCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTca                >  2:850033/3‑150 (MQ=255)
                                                                                                                                         ggCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTcacc              <  2:1027138/150‑1 (MQ=255)
                                                                                                                                         ggCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTcacc              >  2:1519864/1‑150 (MQ=255)
                                                                                                                                           cACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTcaccac            >  2:1749534/1‑150 (MQ=255)
                                                                                                                                            aCCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCaccacc           <  2:825638/150‑1 (MQ=255)
                                                                                                                                               ggACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCa        <  2:111771/150‑1 (MQ=255)
                                                                                                                                               ggACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCa        <  2:1712008/150‑1 (MQ=255)
                                                                                                                                                    cgggcggTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTAACCACCGCACGCGt   >  2:1210844/1‑150 (MQ=255)
                                                                                                                                                     gggcggTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTg  <  1:1445445/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGTGGTATAGCCAGAAAGCGTAGTGCCAATAATTTCGGCTCCCAGCTTTTGGCATGCCAGGCCGTCTTCCGGCGTTGAGCAGTCGGTCATCGCCAGTAAACCGTGATGGTGAATACGTGCCAGCAGCGTTTCAACAGGCACCGGACGCGGGCGGTCGGTGCCGTCAATGGCGATAATGTCCGCGCCCGCCTGCGCCAGCGCATCAACATCTTCAATATAGGCCGTGATGCGTACCGGAGAATCCTCCAGATCGCGTTTCACAATTCCAATAATCGGCACGCTCACCACCGCACGCGTG  >  NC_000913/3370559‑3370857

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: