Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 709,862 C→G 5.9% R202G (CGC→GGC)  ribF → bifunctional riboflavin kinase/FMN adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947709,8620CG 5.9% 380.6 / 12.3 119R202G (CGC→GGC) ribFbifunctional riboflavin kinase/FMN adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (63/49);  new base G (3/4);  total (66/53)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

Only 100 of 132 total aligned reads displayed.

CCAACACGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGTGATGGG  >  NC_002947/709716‑710011
                                                                                                                                                  |                                                                                                                                                     
ccAACACGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGcc                                                                                                                                                    >  1:2234829/1‑150 (MQ=255)
   acacGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCACGTGCTGCACGGCCAGAAGCTGGCCCGCCAGc                                                                                                                                                 >  2:2234329/1‑150 (MQ=255)
    cacGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCGGGCCAGCt                                                                                                                                                <  2:333351/150‑1 (MQ=255)
     acGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTc                                                                                                                                               >  2:1358581/1‑150 (MQ=255)
      cGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCAACCACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCGGGCCAGCttc                                                                                                                                              <  1:1992140/150‑3 (MQ=255)
       ggTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGATGGCCCGCCAGCTCgg                                                                                                                                             >  1:2294830/1‑150 (MQ=255)
           aTCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGcaca                                                                                                                                         >  1:164317/1‑150 (MQ=255)
            tCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCacac                                                                                                                                        <  2:200514/150‑1 (MQ=255)
              cAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCa                                                                                                                                      <  2:1739470/150‑1 (MQ=255)
                ggATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCATGTCGGCACACCAAc                                                                                                                                    <  2:1594803/150‑1 (MQ=255)
                      ccTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCaa                                                                                                                              <  1:108944/150‑1 (MQ=255)
                      ccTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGGACCGCATTACAGGGCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCGCACCAACCGCCaa                                                                                                                              >  2:1301295/1‑150 (MQ=255)
                          cGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCGCTTCATCCTCCAACATc                                                                                                                          <  2:2434554/150‑1 (MQ=255)
                          cGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATc                                                                                                                          <  2:883010/150‑1 (MQ=255)
                                agcagcACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTc                                                                                                                    <  1:922586/150‑1 (MQ=255)
                                agcagcACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTc                                                                                                                    <  1:1630705/150‑1 (MQ=255)
                                 gcagcaCCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCa                                                                                                                   >  2:540417/1‑150 (MQ=255)
                                    gcaCCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGcaca                                                                                                                                         >  1:956617/1‑125 (MQ=255)
                                    gcaCCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGcaca                                                                                                                                         <  2:956617/125‑1 (MQ=255)
                                    gcaCCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAgc                                                                                                                <  2:1757009/150‑1 (MQ=255)
                                     caCCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGACCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAgcg                                                                                                               >  1:391879/1‑150 (MQ=255)
                                                  aaGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCg                                                                                                  >  1:4968/1‑150 (MQ=255)
                                                      ccTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGt                                                                                              >  2:125001/1‑150 (MQ=255)
                                                       cTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTc                                                                                             >  2:224859/1‑150 (MQ=255)
                                                        ttGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCGGGCCAGCttctgg                                                                                                                                           <  2:396284/103‑6 (MQ=255)
                                                        ttGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCGGGCCAGCttctgg                                                                                                                                           >  1:396284/1‑98 (MQ=255)
                                                        ttGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTcc                                                                                            >  2:1420117/1‑150 (MQ=255)
                                                           tCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCggg                                                                                         <  2:290789/150‑1 (MQ=255)
                                                           tCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCggg                                                                                         <  2:2446985/150‑1 (MQ=255)
                                                             cGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCgggg                                                                                        <  2:1563909/149‑1 (MQ=255)
                                                             cGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCgggg                                                                                        >  1:1563909/1‑149 (MQ=255)
                                                             cGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGt                                                                                       >  1:131277/1‑150 (MQ=255)
                                                             cGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGt                                                                                       >  1:353031/1‑150 (MQ=255)
                                                                 ggCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATAc                                                                                   >  2:659298/1‑150 (MQ=255)
                                                                 ggCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATAc                                                                                   <  1:219628/150‑1 (MQ=255)
                                                                  gCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCAAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATAcc                                                                                  >  2:1650930/1‑150 (MQ=255)
                                                                    gACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTg                                                                                >  2:643547/1‑150 (MQ=255)
                                                                       ttCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCGGGCCAGCttc                                                                                                                                              >  2:750259/1‑83 (MQ=255)
                                                                       ttCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCGGGCCAGCttc                                                                                                                                              <  1:750259/85‑3 (MQ=255)
                                                                       ttCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGcc                                                                             >  1:1507163/1‑150 (MQ=255)
                                                                             cTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATc                                                                       <  1:540417/150‑1 (MQ=255)
                                                                             cTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGATGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATc                                                                       <  2:2294830/150‑1 (MQ=255)
                                                                              tGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGcgccgc                                                                                                           <  2:1889160/113‑1 (MQ=255)
                                                                                  cGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCCCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCGGGCCAGCATCGTGAt                                                                  >  2:319816/1‑150 (MQ=255)
                                                                                   gAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTgg                                                                               >  1:432555/1‑136 (MQ=255)
                                                                                    aaCACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCacac                                                                                                                                        <  1:2030204/78‑1 (MQ=255)
                                                                                     acacCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTg                                                                                <  2:812944/133‑1 (MQ=255)
                                                                                      cacCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGAc                                                              <  2:2234829/150‑1 (MQ=255)
                                                                                      cacCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGAc                                                              <  1:224859/150‑1 (MQ=255)
                                                                                      cacCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGAc                                                              <  1:2313861/150‑1 (MQ=255)
                                                                                        ccTGCTGGGCCGGCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGGCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACgg                                                            <  1:1566189/150‑1 (MQ=255)
                                                                                          tgctgGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCa                                                          >  1:186227/1‑150 (MQ=255)
                                                                                            ctgGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCaaa                                                        <  2:2220732/150‑1 (MQ=255)
                                                                                            ctgGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCCGGGTATACCTGGCCAGCATCGTGATCGACGGCaaa                                                        <  1:788143/150‑1 (MQ=255)
                                                                                                 ccgccCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTg                                                   <  1:2443721/150‑1 (MQ=255)
                                                                                                   gccCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCCACCCCCAGCGCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGc                                                 >  1:2122638/1‑150 (MQ=255)
                                                                                                    cccGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGcc                                                >  1:2171985/1‑150 (MQ=255)
                                                                                                      cGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCACCCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCgg                                              <  2:1463363/150‑1 (MQ=255)
                                                                                                      cGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCgg                                              <  1:342731/150‑1 (MQ=255)
                                                                                                      cGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCgg                                              <  1:1420117/150‑1 (MQ=255)
                                                                                                           cGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACg                                                             >  2:264652/1‑130 (MQ=255)
                                                                                                           cGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACg                                                             <  1:264652/130‑1 (MQ=255)
                                                                                                           cGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTc                                         >  1:1365800/1‑150 (MQ=255)
                                                                                                           cGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTc                                         >  1:1196777/1‑150 (MQ=255)
                                                                                                             cATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATGCCTTGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCgg                                       >  2:2098786/1‑150 (MQ=255)
                                                                                                             cATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCgg                                       <  1:2419402/150‑1 (MQ=255)
                                                                                                               ttANAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCa                                     >  2:1502325/1‑150 (MQ=255)
                                                                                                                tACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTAGACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCaa                                    >  1:877424/1‑150 (MQ=255)
                                                                                                                   aGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAATATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAatat                                 >  1:2282096/1‑150 (MQ=255)
                                                                                                                   aGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAatat                                 <  1:286433/150‑1 (MQ=255)
                                                                                                                     gCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCg                               >  1:554298/1‑150 (MQ=255)
                                                                                                                       cgcnTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGc                             <  1:1200672/150‑1 (MQ=255)
                                                                                                                       cgcgTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGc                             <  2:353031/150‑1 (MQ=255)
                                                                                                                          gTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTg                          >  2:1747995/1‑150 (MQ=255)
                                                                                                                             ctgcACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTgcgc                       <  1:1094237/150‑1 (MQ=255)
                                                                                                                             ctgcACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTgcgc                       >  2:517805/1‑150 (MQ=255)
                                                                                                                                 aCGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCAAGCCCAAGCGCCGCCGCGTGCCGCTGTCCGGGGCATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCCCGCGCCCCa                   >  2:414788/1‑150 (MQ=255)
                                                                                                                                   ggCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTg                                                   >  1:1366734/1‑116 (MQ=255)
                                                                                                                                   ggCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTg                                                   <  2:1366734/116‑1 (MQ=255)
                                                                                                                                     ccAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCCAAGCTACCGCGCATTTCGGCAATATCGGCGTGCGCCCCACCGt               >  2:190198/1‑150 (MQ=255)
                                                                                                                                      cAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTg                                                   >  1:311640/1‑113 (MQ=255)
                                                                                                                                       aGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGTGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTg             >  1:784687/1‑150 (MQ=255)
                                                                                                                                       aGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTg             >  1:1887008/1‑150 (MQ=255)
                                                                                                                                        gAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGc            <  2:2122638/150‑1 (MQ=255)
                                                                                                                                        gAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGc            >  1:1612146/1‑150 (MQ=255)
                                                                                                                                         ttcagcgCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGGTGcc           >  2:2124799/7‑150 (MQ=255)
                                                                                                                                          aGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCg          <  2:2166506/150‑1 (MQ=255)
                                                                                                                                           gCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCgg         >  1:354210/1‑150 (MQ=255)
                                                                                                                                            cTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGt        >  1:1884379/1‑150 (MQ=255)
                                                                                                                                              ttcccTCCATCTCCTCACACCAACAGCCACCGTCAAGCTCAGGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGTGa      <  2:391879/148‑1 (MQ=255)
                                                                                                                                              ggCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGTTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGTGa      <  1:1291377/150‑1 (MQ=255)
                                                                                                                                                 ccGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAATCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGTGATgg   <  1:125001/150‑1 (MQ=255)
                                                                                                                                                 ccGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGc                             >  1:2175475/1‑124 (MQ=255)
                                                                                                                                                 ccGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGc                             <  2:2175475/124‑1 (MQ=255)
                                                                                                                                                 ccGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTgc                         <  1:232559/128‑1 (MQ=255)
                                                                                                                                                 ccGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTgc                         <  1:235184/128‑1 (MQ=255)
                                                                                                                                                 ccGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTgc                         >  2:235184/1‑128 (MQ=255)
                                                                                                                                                 ccGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTgc                         >  2:232559/1‑128 (MQ=255)
                                                                                                                                                 ccGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGTGATgg   >  1:293998/1‑150 (MQ=255)
                                                                                                                                                 ccGCCAGCCCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGTGATgg   <  2:4968/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CCAACACGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGTGATGGG  >  NC_002947/709716‑710011

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