Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,154,434 T→G 5.1% intergenic (‑423/+25) PP_mr51 ← / ← PP_4535 C4 ncRNA/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,154,4340TG 5.1% 1243.2 / 36.5 361intergenic (‑423/+25)PP_mr51/PP_4535C4 ncRNA/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (201/141);  new base G (16/3);  total (217/144)
Fisher's exact test for biased strand distribution p-value = 3.03e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 370 total aligned reads displayed.

TAACTATGCGCGCGTCTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCTTGACGCCGGTTT  >  NC_002947/5154285‑5154583
                                                                                                                                                     |                                                                                                                                                     
tAACTATGCGCGCGTCTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCt                                                                                                                                                       >  1:2583430/1‑150 (MQ=255)
  aCTATGCGCGCGTCTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTAt                                                                                                                                                     >  1:5935927/1‑150 (MQ=255)
  aCTATGCGCGCGTCTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTAt                                                                                                                                                     <  1:4589237/150‑1 (MQ=255)
   cTATGCGCGCGTCTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATc                                                                                                                                                    <  2:7285254/150‑1 (MQ=255)
    tATGCGCGCGTCTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATcc                                                                                                                                                   >  1:2251455/1‑150 (MQ=255)
        cgcgcgTCTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATcccccc                                                                                                                                               <  1:3958762/150‑1 (MQ=255)
         gcgcgTCTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATccccccc                                                                                                                                              >  2:8637917/1‑150 (MQ=255)
           gcgTCTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTg                                                                                                                                            >  2:7076261/1‑150 (MQ=255)
             gTCTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGt                                                                                                                                          >  2:5714388/1‑150 (MQ=255)
              tCTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGtt                                                                                                                                         <  2:8784131/150‑1 (MQ=255)
              tCTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGtt                                                                                                                                         <  2:4513608/150‑1 (MQ=255)
              tCTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGtt                                                                                                                                         <  1:762301/150‑1 (MQ=255)
               cTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTg                                                                                                                                        <  2:149575/150‑1 (MQ=255)
                ttGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTgg                                                                                                                                       <  1:7127772/150‑1 (MQ=255)
                  gAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGaa                                                                                                                                     <  2:3985871/150‑1 (MQ=255)
                     tCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACtt                                                                                                                                  >  2:5368655/1‑150 (MQ=255)
                      cGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACttg                                                                                                                                 <  2:2570660/150‑2 (MQ=255)
                        cTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAg                                                                                                                               <  2:2029930/150‑1 (MQ=255)
                               ttCTGAATGACCACGAGCACTCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCAGATCTTGGCTTCAACTTCAGGGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCg                                                                                                                        <  2:7719611/150‑1 (MQ=255)
                               ttCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCg                                                                                                                        <  1:7310018/150‑1 (MQ=255)
                                tCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCgg                                                                                                                       <  1:2807167/150‑1 (MQ=255)
                                 cTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGc                                                                                                                      <  1:1024861/150‑1 (MQ=255)
                                  tGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGcc                                                                                                                     >  1:7256114/1‑149 (MQ=255)
                                  tGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCa                                                                                                                     >  2:8211144/1‑150 (MQ=255)
                                   gAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCaa                                                                                                                    >  2:3627549/1‑150 (MQ=255)
                                   gAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCaa                                                                                                                    <  2:6041802/150‑1 (MQ=255)
                                        aCCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGa                                                                                                               <  1:2095359/150‑1 (MQ=255)
                                         ccACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGAc                                                                                                              >  1:8034321/1‑150 (MQ=255)
                                         ccACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGAc                                                                                                              >  2:2235298/1‑150 (MQ=255)
                                          cACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGAcc                                                                                                             >  2:2729179/1‑150 (MQ=255)
                                          cACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGAcc                                                                                                             >  1:2641868/1‑150 (MQ=255)
                                                 aaTCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGtt                                                                                                      >  1:1397211/1‑150 (MQ=255)
                                                     aGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTcc                                                                                                  >  1:3654925/1‑150 (MQ=255)
                                                          cATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGa                                                                                             <  2:1075212/150‑1 (MQ=255)
                                                           ataataCATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAg                                                                                            >  1:4330144/1‑150 (MQ=255)
                                                               taCATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGcgcg                                                                                        >  2:4683499/1‑150 (MQ=255)
                                                                aCATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCgg                                                                                       <  1:513395/150‑1 (MQ=255)
                                                                 cATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCACGGTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGa                                                                                      <  2:5220030/150‑1 (MQ=255)
                                                                  aTAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGAc                                                                                     <  1:4861631/150‑1 (MQ=255)
                                                                  aTAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGAc                                                                                     <  1:792861/150‑1 (MQ=255)
                                                                   tAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACg                                                                                    >  2:7804860/1‑150 (MQ=255)
                                                                     aTGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGaa                                                                                  <  2:7631659/150‑1 (MQ=255)
                                                                     aTGGGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCACGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACgca                                                                                  <  1:4730169/150‑3 (MQ=255)
                                                                       gtgtCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGt                                                                                >  1:3361747/1‑150 (MQ=255)
                                                                          tCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCt                                                                             >  1:5719072/1‑150 (MQ=255)
                                                                           ccACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTc                                                                            <  2:807795/150‑1 (MQ=255)
                                                                              cgccgcCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGaa                                                                         <  2:6265248/150‑1 (MQ=255)
                                                                              cgccgcCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGaa                                                                         <  1:739505/150‑1 (MQ=255)
                                                                              cgccgcCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGaa                                                                         <  2:4925677/150‑1 (MQ=255)
                                                                                ccgcccgcCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGa                                                                       >  1:6741452/1‑150 (MQ=255)
                                                                                       ccATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGtt                                                                >  1:3710119/1‑150 (MQ=255)
                                                                                       ccATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGtt                                                                <  2:1387830/150‑1 (MQ=255)
                                                                                       ccATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGtt                                                                <  2:7480753/150‑1 (MQ=255)
                                                                                          tatCTTGTCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGcc                                                             <  1:3163/150‑1 (MQ=255)
                                                                                          tatCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGGTGcc                                                             <  1:8427235/150‑1 (MQ=255)
                                                                                            tCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCACGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAAATTCAGTTACTCAGCAATCAGACCTCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATTGTTGCCCt                                                           <  2:6077069/150‑1 (MQ=255)
                                                                                             cTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTg                                                          <  1:4264198/150‑1 (MQ=255)
                                                                                                  cTTCAACTTCAGTGTCCAAGCCAATCATCATCAAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAAATTCAGTTACTCGGCAATAAGACCGCGCCTTCTCAACCGAGCGCGGACGAAGTCCTCGAAGACATAGTAGCCCTGCAACg                                                     >  1:8357293/1‑150 (MQ=255)
                                                                                                   ttCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACgg                                                    >  1:8561476/1‑150 (MQ=255)
                                                                                                     cAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCa                                                  >  1:5257475/1‑150 (MQ=255)
                                                                                                      aaCTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCACCGGCAg                                                 >  1:4152279/1‑150 (MQ=255)
                                                                                                      aaCTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAg                                                 >  1:587577/1‑150 (MQ=255)
                                                                                                       ccTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAgg                                                <  1:7544471/149‑1 (MQ=255)
                                                                                                              tgtCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCcg                                         >  1:6784493/1‑150 (MQ=255)
                                                                                                              tgtCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCcg                                         >  1:1595664/1‑150 (MQ=255)
                                                                                                               gtCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCcgc                                        <  2:1193033/150‑1 (MQ=255)
                                                                                                               gtCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCcgc                                        <  1:951033/150‑1 (MQ=255)
                                                                                                                   aaGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGa                                    >  1:3517800/1‑150 (MQ=255)
                                                                                                                   aaGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGa                                    >  2:5544345/1‑150 (MQ=255)
                                                                                                                   aaGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGa                                    >  1:18079/1‑150 (MQ=255)
                                                                                                                   aaGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGACAGGCCATCCGCGTGa                                    >  2:844147/1‑150 (MQ=255)
                                                                                                                      ccAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCg                                 <  1:6679771/150‑1 (MQ=255)
                                                                                                                        aaTCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCa                               >  1:4991632/1‑150 (MQ=255)
                                                                                                                          tCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGACCTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGa                             >  2:7417402/1‑150 (MQ=255)
                                                                                                                            cTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCACCGGCAGGCCATCCGCGTGATCGCAGACa                           >  1:1918775/1‑150 (MQ=255)
                                                                                                                             tCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACaa                          >  2:3416205/1‑150 (MQ=255)
                                                                                                                               aTCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATg                        >  2:3231675/1‑150 (MQ=255)
                                                                                                                               aTCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATg                        >  2:545809/1‑150 (MQ=255)
                                                                                                                                  cAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGAtt                     >  2:5276860/1‑150 (MQ=255)
                                                                                                                                  cAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGAtt                     >  2:1689468/1‑150 (MQ=255)
                                                                                                                                  cAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGAtt                     >  2:3686693/1‑150 (MQ=255)
                                                                                                                                      ttCCAACCAGGGGGGGATAGCCTCCTGGTTGTAAATTCAGTTACTCGGCAATCAGACCGCGCGTTATCCACCGAGCGAGGACGAAGTCCTCGAAGACATAGTAGCAATGCAACGGCAGGCCAACCGCGTGATCGCAGACAATGATTTGCg                 >  1:2477363/1‑150 (MQ=255)
                                                                                                                                      ttCCAACCAGGGGGGGATAGCCTCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCg                 >  1:5514220/1‑150 (MQ=39)
                                                                                                                                      ttCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCg                 <  1:2964955/150‑1 (MQ=255)
                                                                                                                                      ttCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCg                 >  1:2154697/1‑150 (MQ=255)
                                                                                                                                       tCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCgg                <  1:6447134/150‑1 (MQ=255)
                                                                                                                                       tCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCgg                >  2:6885667/1‑150 (MQ=255)
                                                                                                                                         cAACCAGGGGGGGATAGCCTCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCt              >  2:6802331/1‑150 (MQ=255)
                                                                                                                                         cAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCt              >  2:5860894/1‑150 (MQ=255)
                                                                                                                                         cAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCt              >  1:8759046/1‑150 (MQ=255)
                                                                                                                                         cAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCt              >  1:7007967/1‑150 (MQ=255)
                                                                                                                                          aaCCAGGGGGGGATAGCCTCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACaa                          >  1:5363045/1‑137 (MQ=39)
                                                                                                                                          aaCCAGGGGGGGATAGCCTCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCtt             >  2:1766616/1‑150 (MQ=255)
                                                                                                                                          aaCCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCtt             >  1:122374/1‑150 (MQ=255)
                                                                                                                                          aaCCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCtt             >  2:2170944/1‑150 (MQ=255)
                                                                                                                                            ccAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCTTGa           >  2:513332/1‑150 (MQ=255)
                                                                                                                                             cAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATGGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCTTGAc          >  1:60672/1‑150 (MQ=255)
                                                                                                                                               ggaggCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCTTGACGc        >  1:8299706/1‑150 (MQ=255)
                                                                                                                                                     tATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCTTGACGCCGGttt  >  2:4411863/1‑150 (MQ=255)
                                                                                                                                                     tATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCTTGACGCCGGttt  >  2:4837623/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TAACTATGCGCGCGTCTTGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCTTGACGCCGGTTT  >  NC_002947/5154285‑5154583

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: