Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 851,027 T→C 5.2% intergenic (‑607/‑198) pcs ← / → hemA phosphatidylcholine synthase/glutamyl‑tRNA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947851,0270TC 5.2% 545.9 / 11.6 171intergenic (‑607/‑198)pcs/hemAphosphatidylcholine synthase/glutamyl‑tRNA reductase
Reads supporting (aligned to +/- strand):  ref base T (93/69);  new base C (2/7);  total (95/76)
Fisher's exact test for biased strand distribution p-value = 7.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 175 total aligned reads displayed.

AAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACGGTGGTTGTCCAAAAACC  >  NC_002947/850878‑851176
                                                                                                                                                     |                                                                                                                                                     
aaGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGt                                                                                                                                                       >  2:3095822/1‑150 (MQ=255)
  gcgcgcGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTc                                                                                                                                                     <  1:1409462/150‑1 (MQ=255)
    gcgcgGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCgg                                                                                                                                                   <  2:4189291/150‑1 (MQ=255)
         gTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGcctcc                                                                                                                                              >  1:3557536/1‑150 (MQ=255)
          ttGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCCTCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGcctcct                                                                                                                                             <  2:5120736/150‑1 (MQ=255)
             aTCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGaa                                                                                                                                          <  2:3571132/150‑1 (MQ=255)
              tCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAg                                                                                                                                         <  1:250295/150‑1 (MQ=255)
                  gcgcTGCAGATTGCGACTTCAGGAGGCCGAGGGTAGGGGTTGCTCCCGGGAGGACCGCGTTGCCGGCCCCATCCTCGGGCACAGCGCTTGCCGGGCATCCATGCCCGTCACCTCCCTGCTCAAGGCCTGCGTTCGGCCTCCTGAAGTcgc                                                                                                                                     <  1:396300/150‑1 (MQ=21)
                  gcgcTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTcgc                                                                                                                                     >  2:109259/1‑150 (MQ=255)
                     cTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaa                                                                                                                                  <  2:2957383/150‑1 (MQ=255)
                            ttGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaag                                                                                                                           <  2:5383561/150‑1 (MQ=255)
                                gACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAagatcaagatca                                                                                                                       >  2:3467050/1‑150 (MQ=255)
                                 aCTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAg                                                                                                                      <  2:5369663/150‑1 (MQ=255)
                                  cTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGGTCAgg                                                                                                                     >  1:3662575/1‑150 (MQ=255)
                                   ttCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGa                                                                                                                    <  1:4000257/150‑1 (MQ=255)
                                         aggCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTGGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCaaaa                                                                                                              <  2:2929175/150‑1 (MQ=255)
                                             cGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCa                                                                                                          <  1:3163431/150‑1 (MQ=255)
                                             cGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCa                                                                                                          <  2:1434265/150‑1 (MQ=255)
                                                  gTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCAAGATCa                                                                                                     >  1:4375036/1‑150 (MQ=255)
                                                   tAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCaa                                                                                                    >  1:1455508/1‑150 (MQ=255)
                                                    aGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCACCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCaac                                                                                                   <  2:3557536/150‑2 (MQ=255)
                                                      ggCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAg                                                                                                 >  1:100799/1‑150 (MQ=255)
                                                       gCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGc                                                                                                >  2:3522157/1‑150 (MQ=255)
                                                       gCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGc                                                                                                >  1:3289522/1‑150 (MQ=255)
                                                       gCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGAGCAAGATCAGGATCAAAAGCCAGATCAAAAGc                                                                                                >  1:2243174/1‑150 (MQ=255)
                                                        cTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGCCGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGcc                                                                                               <  1:2170609/150‑1 (MQ=255)
                                                        cTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGcc                                                                                               <  1:2170671/150‑1 (MQ=255)
                                                        cTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGcc                                                                                               >  2:1450676/1‑150 (MQ=255)
                                                          tGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGAACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCgg                                                                                             <  2:2073713/150‑1 (MQ=255)
                                                             tCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCaa                                                                                          <  2:93256/150‑1 (MQ=255)
                                                                 gggAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCACAAGCCGGCAAGAcc                                                                                      >  2:2867391/1‑150 (MQ=255)
                                                                 gggAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGAcc                                                                                      >  1:2508299/1‑150 (MQ=255)
                                                                    aggaCCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGa                                                                                   <  1:109259/150‑1 (MQ=255)
                                                                      tgacgggcatggatgcccgtcaagcgctgaggccccaggatggggccttcggcgcggtcctgcctccgcaaggCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatc                                                                                                                  >  2:4915015‑M1/74‑117 (MQ=255)
                                                                      tgacgggcatggatgcccgtcaagcgctgaggccccaggatggggccttcggcgcggtcctgcctccgcaaggCCTGCGTTCGGCCTCCTGAAGTCGCGaagatcaagatcaagatc                                                                                                                  <  1:4915015‑M1/44‑1 (MQ=255)
                                                                      gaCCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGaaa                                                                                 >  1:4008814/1‑150 (MQ=255)
                                                                       aCCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGaaaa                                                                                >  2:230587/1‑150 (MQ=255)
                                                                         cgcgCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAaca                                                                              >  1:1476172/1‑150 (MQ=255)
                                                                         cgcgCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAaca                                                                              >  2:4486563/1‑150 (MQ=255)
                                                                         cgcgCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAaca                                                                              >  1:2830174/1‑150 (MQ=255)
                                                                          gcgcTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAtcac                                                                             <  1:3078544/150‑1 (MQ=255)
                                                                          gcgcTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAacac                                                                             <  2:4423666/150‑1 (MQ=255)
                                                                           cgcTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACAcc                                                                            >  2:1963781/1‑150 (MQ=255)
                                                                                ggatgcccgtcaagcgctgtgcccgaggatggggcctgcagcgcggtcctcccgggagcaagcCCTACGCTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAg                                                                       <  1:271317‑M1/87‑1 (MQ=255)
                                                                                ggatgcccgtcaagcgctgtgcccgaggatggggcctgcagcgcggtcctcccgggagcaagcCCTACGCTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAg                                                                       <  1:4124974‑M1/87‑1 (MQ=255)
                                                                                 aGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAgg                                                                      <  2:3432723/150‑1 (MQ=255)
                                                                                  ggCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAggg                                                                     >  1:2725390/1‑150 (MQ=255)
                                                                                     cccATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGcc                                                                                               <  2:1968263/121‑1 (MQ=255)
                                                                                     cccATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGcc                                                                                               >  1:1968263/1‑121 (MQ=255)
                                                                                         tCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCa                                                                                                          <  2:1607931/106‑1 (MQ=255)
                                                                                          ccTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTAc                                                             <  1:5184379/150‑1 (MQ=255)
                                                                                            tCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGa                                                           >  2:2191108/1‑150 (MQ=255)
                                                                                               ggCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGaaaa                                                        >  1:3184577/1‑150 (MQ=255)
                                                                                                gCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAg                                                       >  1:1603437/1‑150 (MQ=255)
                                                                                                gCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAg                                                       <  1:4430799/150‑1 (MQ=255)
                                                                                                  acaGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGc                                                     >  1:5156177/1‑150 (MQ=255)
                                                                                                  acaGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGc                                                     >  2:4153316/1‑150 (MQ=255)
                                                                                                    aGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGa                                                   >  1:4498029/1‑150 (MQ=255)
                                                                                                    aGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGa                                                   >  1:1799884/1‑150 (MQ=255)
                                                                                                    aGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGa                                                   >  1:2349199/1‑150 (MQ=255)
                                                                                                     gcgcTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGAt                                                  >  2:4044774/1‑150 (MQ=255)
                                                                                                        gaggatggggcctgcagcgcggtcctcccgggagcaagcCCTACGCTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGc                                               <  2:3092156‑M1/111‑1 (MQ=255)
                                                                                                        cTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGGCCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGc                                               <  1:986646/150‑1 (MQ=255)
                                                                                                           gACGGGCATGGATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAATGCCAGATCAAAAGCCGGCAAGTCCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGc                                            >  2:2252835/1‑150 (MQ=255)
                                                                                                            aCGGGCATGGATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCt                                           >  2:1367961/1‑150 (MQ=255)
                                                                                                              ggggcctgcagcgcggtcctcccgggagcaagcCCTACGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTgc                                         <  1:4833041‑M1/117‑1 (MQ=255)
                                                                                                              gggCATGGATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTgc                                         >  2:1239976/1‑150 (MQ=255)
                                                                                                              gggCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTgc                                         >  1:388630/1‑150 (MQ=255)
                                                                                                              gggCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTgc                                         >  1:247755/1‑150 (MQ=255)
                                                                                                                gCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGTGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTgcgc                                       >  1:4188739/1‑150 (MQ=255)
                                                                                                                gCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGTGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTgcgc                                       >  1:4188740/1‑150 (MQ=255)
                                                                                                                gCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGTAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTgcgc                                       >  1:818543/1‑150 (MQ=255)
                                                                                                                    ggaTGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATc                                   >  2:4470426/3‑150 (MQ=255)
                                                                                                                    ccATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTAt                                    >  2:1277424/1‑149 (MQ=255)
                                                                                                                    ccATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTAt                                    <  1:1277424/149‑1 (MQ=255)
                                                                                                                     cATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCa                                  >  1:1437599/1‑150 (MQ=255)
                                                                                                                      aTGCCCGTCACCTCCCGGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAg                                 >  2:2984488/1‑150 (MQ=255)
                                                                                                                             tCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATa                          >  2:3655621/1‑150 (MQ=255)
                                                                                                                                ccTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGca                               >  1:4512589/1‑142 (MQ=255)
                                                                                                                                ccTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCg                       >  2:5084159/1‑150 (MQ=255)
                                                                                                                                 cTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGAt                                                                                                       >  2:4392074/1‑69 (MQ=255)
                                                                                                                                  tCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTCGGGCAGGCTACGAAAAGGCGATAGCCGCGGCGCGAGCAGCACCATAGCGag                     >  2:3174801/1‑149 (MQ=38)
                                                                                                                                  tCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGaa                     <  2:2134111/150‑1 (MQ=255)
                                                                                                                                  tCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGaa                     <  2:388630/150‑1 (MQ=255)
                                                                                                                                    cgggagcaagcCCTACGCTCGGCCTCCTGAAGTCGCAAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACg                   <  1:1525929‑M1/139‑1 (MQ=255)
                                                                                                                                    ccTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCg                                                                           >  2:3636339/1‑94 (MQ=255)
                                                                                                                                    ccTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCAACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACg                   >  1:714050/1‑150 (MQ=255)
                                                                                                                                       gcgcATGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACggtg                <  2:4596004/150‑1 (MQ=255)
                                                                                                                                       gcgcAAGGCCTGGGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACggtg                <  1:1144924/150‑1 (MQ=255)
                                                                                                                                       gcgcAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACggtg                <  2:1788047/150‑1 (MQ=255)
                                                                                                                                       gcgcAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACggtg                <  1:3095822/150‑1 (MQ=255)
                                                                                                                                         gcAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACggtggt              <  1:4617462/150‑1 (MQ=255)
                                                                                                                                            aGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACGGTGGTTGt           <  2:1437977/150‑1 (MQ=255)
                                                                                                                                             ggCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACGGTGGTTGTc          <  2:4825131/150‑1 (MQ=255)
                                                                                                                                             ggCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACGGTGGTTGTc          <  1:3596204/150‑1 (MQ=255)
                                                                                                                                              gCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACGGTGGTTGTcc         >  1:3201269/1‑150 (MQ=255)
                                                                                                                                                CTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACGGTGGTTGTCCaa       <  2:4577104‑M1/150‑1 (MQ=255)
                                                                                                                                                 tGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAATCCGGCAATACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGAGATAGCCGCTTCGCTATCAGCACAATAGAGAACGGTGGTTGTCCaaa      >  2:813108/1‑150 (MQ=38)
                                                                                                                                                 tGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACGGTGGTTGTCCaaa      >  2:2523795/1‑150 (MQ=255)
                                                                                                                                                    gTTCGGCCTCCTTATGTCGCGAAGATAAAGATCAGGATCAAAAGCCACATCAAAAGCCGGAAAGACCTTAAAACGCCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACGGTGGTTGTCCAAAAAc   <  1:4729458/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAGCGCGCGGTTGATCAGGCGCTGCAGATTGCGACTTCAGGAGGCCGAGCGTAGGGCTTGCTCCCGGGAGGACCGCGCTGCAGGCCCCATCCTCGGGCACAGCGCTTGACGGGCATCCATGCCCGTCACCTCCCTGCGCAAGGCCTGCGTTCGGCCTCCTGAAGTCGCGAAGATCAAGATCAGGATCAAAAGCCAGATCAAAAGCCGGCAAGACCTGAAAACACCGCTAGGGCAGGCTACGAAAAGGCGATAGCCGCTGCGCTATCAGCACAATAGCGAACGGTGGTTGTCCAAAAACC  >  NC_002947/850878‑851176

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