Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 709,861 C→G 5.9% A201A (GCC→GCG ribF → bifunctional riboflavin kinase/FMN adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947709,8610CG 5.9% 373.0 / 11.6 119A201A (GCC→GCGribFbifunctional 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.90e-01

Only 100 of 133 total aligned reads displayed.

CCAACACGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGTGATGG  >  NC_002947/709716‑710010
                                                                                                                                                 |                                                                                                                                                     
ccAACACGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGcc                                                                                                                                                   >  1:2234829/1‑150 (MQ=255)
   acacGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCGGGCCAGc                                                                                                                                                >  1:333351/1‑150 (MQ=255)
   acacGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCACGTGCTGCACGGCCAGAAGCTGGCCCGCCAGc                                                                                                                                                >  2:2234329/1‑150 (MQ=255)
    cacGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCGGGCCAGCt                                                                                                                                               <  2:333351/150‑1 (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)
                    ggCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGcc                                                                                                                               >  1:2310007/1‑150 (MQ=255)
                      ccTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCaa                                                                                                                             <  1:108944/150‑1 (MQ=255)
                      ccTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGATCGGCACACCAACAGCCaa                                                                                                                             <  2:2233597/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)
                            ggTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGcc                                                                                                                                                       <  2:223630/118‑1 (MQ=255)
                            ggTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGcc                                                                                                                                                       >  1:223630/1‑118 (MQ=255)
                                agcagcACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTc                                                                                                                   <  1:1630705/150‑1 (MQ=255)
                                agcagcACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTc                                                                                                                   <  1:922586/150‑1 (MQ=255)
                                 gcagcaCCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCa                                                                                                                  >  2:540417/1‑150 (MQ=255)
                                    gcaCCGAAGTGCGCCAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAgc                                                                                                               <  1:713642/150‑1 (MQ=255)
                                    gcaCCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGcaca                                                                                                                                        <  2:956617/125‑1 (MQ=255)
                                    gcaCCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAgc                                                                                                               <  2:1757009/150‑1 (MQ=255)
                                     caCCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGACCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAgcg                                                                                                              >  1:391879/1‑150 (MQ=255)
                                         gAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCCGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGcgccgc                                                                                                          <  1:267653/150‑1 (MQ=255)
                                                  aaGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCg                                                                                                 >  1:4968/1‑150 (MQ=255)
                                                      ccTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGt                                                                                             >  2:125001/1‑150 (MQ=255)
                                                       cTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTc                                                                                            >  2:224859/1‑150 (MQ=255)
                                                        ttGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCGGGCCAGCttctgg                                                                                                                                          >  1:396284/1‑98 (MQ=255)
                                                        ttGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCGGGCCAGCttctgg                                                                                                                                          <  2:396284/103‑6 (MQ=255)
                                                        ttGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTcc                                                                                           >  2:1420117/1‑150 (MQ=255)
                                                           tCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCggg                                                                                        <  2:2446985/150‑1 (MQ=255)
                                                           tCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCggg                                                                                        <  2:290789/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                                                                                                                                             <  1:750259/85‑3 (MQ=255)
                                                                       ttCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCGGGCCAGCttc                                                                                                                                             >  2:750259/1‑83 (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)
                                                                                    aaCACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCacac                                                                                                                                       >  2:2030204/1‑78 (MQ=255)
                                                                                    aaCACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCacac                                                                                                                                       <  1:2030204/78‑1 (MQ=255)
                                                                                    aaCACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACTTGGCCAGCATCGTGATCg                                                               >  2:179235/1‑150 (MQ=255)
                                                                                     acacCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTg                                                                               >  1:812944/1‑133 (MQ=255)
                                                                                     acacCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTg                                                                               <  2:812944/133‑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)
                                                                                                 ccgccCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTg                                                  <  1:2443721/150‑1 (MQ=255)
                                                                                                 ccgccCGTACCGCATTAAAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTg                                                  >  2:1408886/1‑150 (MQ=255)
                                                                                                   gccCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCCACCCCCAGCGCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGc                                                >  1:2122638/1‑150 (MQ=255)
                                                                                                    cccGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGcc                                               >  1:2171985/1‑150 (MQ=255)
                                                                                                      cGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCgg                                             <  1:342731/150‑1 (MQ=255)
                                                                                                      cGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCgg                                             <  1:1420117/150‑1 (MQ=255)
                                                                                                       gTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCGTGGCCggg                                            >  1:378525/1‑150 (MQ=255)
                                                                                                        tACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGTTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGgt                                           >  2:1291377/1‑150 (MQ=255)
                                                                                                           cGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACg                                                            >  2:264652/1‑130 (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)
                                                                                                                   aGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAatat                                <  1:286433/150‑1 (MQ=255)
                                                                                                                     gCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGAGCGACGGCAAAGCTTGGCCGGCTGTCGGCAATATCg                              >  2:153740/1‑150 (MQ=255)
                                                                                                                       cgcnTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGc                            <  1:1200672/150‑1 (MQ=255)
                                                                                                                       cgcgTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGc                            <  2:353031/150‑1 (MQ=255)
                                                                                                                       cgcgTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGc                            >  1:454017/1‑150 (MQ=255)
                                                                                                                           tgctgcACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTgc                        >  2:2152270/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)
                                                                                                                                       aGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGTGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTg            >  1:784687/1‑150 (MQ=255)
                                                                                                                                       aGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTg            >  1:1887008/1‑150 (MQ=255)
                                                                                                                                        gAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGc           >  1:1612146/1‑150 (MQ=255)
                                                                                                                                         ttcagcgCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGGTGcc          >  2:2124799/7‑150 (MQ=255)
                                                                                                                                          aGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCg         <  2:2166506/150‑1 (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)
                                                                                                                                                 ccGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTgc                        <  1:235184/128‑1 (MQ=255)
                                                                                                                                                 ccGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTgc                        >  2:232559/1‑128 (MQ=255)
                                                                                                                                                 ccGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTgc                        >  2:235184/1‑128 (MQ=255)
                                                                                                                                                 ccGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGTGATgg  >  2:172970/1‑150 (MQ=255)
                                                                                                                                                 ccGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGTGATgg  >  1:293998/1‑150 (MQ=255)
                                                                                                                                                 ccGCCAGCCCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGTGATgg  <  2:4968/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CCAACACGGTCATCCAGGATGGCCTGCGGGTCAGCAGCACCGAAGTGCGCAAGGCCTTGTCCGAGGGCGACTTCGAACTGGCCGAACACCTGCTGGGCCGCCCGTACCGCATTACAGGCCGCGTGCTGCACGGCCAGAAGCTGGCCCGCCAGCTCGGCACACCAACCGCCAACATCCAGCTCAAGCGCCGCCGCGTGCCGCTGTCCGGGGTATACCTGGCCAGCATCGTGATCGACGGCAAAGCTTGGCCGGGTGTCGGCAATATCGGCGTGCGCCCCACCGTTGCCGGTGATGG  >  NC_002947/709716‑710010

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