Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,926,529 C→T 5.1% G84D (GGC→GAC)  peaB ← quinohemoprotein amine dehydrogenase modification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,926,5290CT 5.1% 699.1 / 16.8 214G84D (GGC→GAC) peaBquinohemoprotein amine dehydrogenase modification protein
Reads supporting (aligned to +/- strand):  ref base C (111/92);  new base T (3/8);  total (114/100)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 222 total aligned reads displayed.

TCATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCA  >  NC_002947/3926380‑3926678
                                                                                                                                                     |                                                                                                                                                     
tCATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGc                                                                                                                                                       >  2:2792076/1‑150 (MQ=255)
tCATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGc                                                                                                                                                       >  2:1173688/1‑150 (MQ=255)
  aTCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCg                                                                                                                                                     >  2:418294/1‑150 (MQ=255)
  aTCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCg                                                                                                                                                     >  2:2895667/1‑150 (MQ=255)
  aTCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCg                                                                                                                                                     <  2:3701830/150‑1 (MQ=255)
  aTCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCg                                                                                                                                                     >  1:3701830/1‑150 (MQ=255)
     ttcttGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCg                                                                                                                                                  <  1:2879422/150‑1 (MQ=255)
         tGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGa                                                                                                                                              <  2:2647106/150‑1 (MQ=255)
         tGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGACGGCTCGCCGCTGa                                                                                                                                              <  1:137695/150‑1 (MQ=255)
         tGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGACGGCTCGCCGCTGa                                                                                                                                              <  1:528738/150‑1 (MQ=255)
          gCCGGCTGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGACGGCTCGCCGCTGAc                                                                                                                                             <  1:3261657/150‑1 (MQ=255)
          gCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGAc                                                                                                                                             <  2:3460974/150‑1 (MQ=255)
            cggcggAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGGAACGGGTGTTGACGTTGAGCACCTCGGTGTTGTGCGCGGTGCGCTCGTCTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCGGTCGCTgtgca                                                                                                                                           <  2:1678463/150‑5 (MQ=255)
             ggcggAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGACGGCTCGCCGCTGAcccc                                                                                                                                          <  2:1717130/150‑3 (MQ=255)
              gcggAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGc                                                                                                                                         >  2:1681051/1‑150 (MQ=255)
                ggAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTc                                                                                                                                       >  1:2747864/1‑150 (MQ=255)
                   aGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAg                                                                                                                                    >  2:4056897/1‑150 (MQ=255)
                     gCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAgcg                                                                                                                                  >  2:4143315/1‑150 (MQ=255)
                     gCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAgcg                                                                                                                                  >  1:3952397/1‑150 (MQ=255)
                       ttGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGcgcg                                                                                                                                <  2:1516227/150‑1 (MQ=255)
                       ttGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGcgcg                                                                                                                                >  2:4501548/1‑150 (MQ=255)
                       ttGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGcgcg                                                                                                                                >  2:2040742/1‑150 (MQ=255)
                         gTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGAt                                                                                                                              <  1:2598002/150‑1 (MQ=255)
                          tCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATc                                                                                                                             >  1:2798004/1‑150 (MQ=255)
                           ccAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCa                                                                                                                            >  1:3193781/1‑150 (MQ=255)
                           ccAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCa                                                                                                                            >  1:1221728/1‑150 (MQ=255)
                             aGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGt                                                                                                                          >  1:2236148/1‑150 (MQ=255)
                              ggTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGtt                                                                                                                         >  1:2472427/1‑150 (MQ=255)
                               gTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTc                                                                                                                        <  2:2211412/150‑1 (MQ=255)
                                tCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTcg                                                                                                                       <  2:3072340/150‑1 (MQ=255)
                                    ccTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCTGCGCGAACTGATCgcgct                                                                                                                   <  1:4230725/150‑2 (MQ=255)
                                      ttGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCTAGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAgg                                                                                                                 >  2:1391317/1‑150 (MQ=255)
                                      ttGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGTGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAgg                                                                                                                 >  1:44358/1‑150 (MQ=255)
                                      ttGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAgg                                                                                                                 <  1:3276373/150‑1 (MQ=255)
                                      ttGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAgg                                                                                                                 <  1:807329/150‑1 (MQ=255)
                                      ttGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAgg                                                                                                                 >  2:1671556/1‑150 (MQ=255)
                                      ttGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAgg                                                                                                                 >  2:3801764/1‑150 (MQ=255)
                                      ttGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAgg                                                                                                                 >  1:2131122/1‑150 (MQ=255)
                                       tGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCTGCGCGAACTGATCGCGCTgga                                                                                                                <  1:4356611/150‑2 (MQ=255)
                                       tGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGAACTGATCGCGCTgga                                                                                                                <  1:1252789/150‑2 (MQ=255)
                                           gCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAggagctg                                                                                                            <  2:2845022/150‑6 (MQ=255)
                                            cAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTcc                                                                                                           <  2:1965734/150‑1 (MQ=255)
                                               tAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTgg                                                                                                        >  2:2327007/1‑150 (MQ=255)
                                               tAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTgg                                                                                                        <  1:2511115/150‑1 (MQ=255)
                                                       gCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTc                                                                                                >  2:781596/1‑150 (MQ=255)
                                                        cTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCg                                                                                               >  1:2260103/1‑150 (MQ=255)
                                                          gAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGcc                                                                                             <  1:673425/150‑1 (MQ=255)
                                                             gTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGc                                                                                          >  1:2409791/1‑150 (MQ=255)
                                                             gTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGc                                                                                          >  1:2982542/1‑150 (MQ=255)
                                                                 cAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAAc                                                                                      >  1:2388886/1‑150 (MQ=255)
                                                                 cAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAAc                                                                                      >  1:3869625/1‑150 (MQ=255)
                                                                    ccGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATc                                                                                                                             <  1:1685064/108‑1 (MQ=255)
                                                                    ccGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAActgc                                                                                   >  1:811722/1‑150 (MQ=255)
                                                                     cGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAActgct                                                                                  >  2:1243314/1‑150 (MQ=255)
                                                                      ggTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAActgctg                                                                                 <  2:599791/150‑1 (MQ=255)
                                                                       gtgtTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTgg                                                                                <  1:1014048/150‑1 (MQ=255)
                                                                        tgtTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGc                                                                               <  1:2654587/150‑1 (MQ=255)
                                                                         gtTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGcc                                                                              >  2:260142/1‑150 (MQ=255)
                                                                             aCGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAgg                                                                          >  2:3570722/1‑150 (MQ=255)
                                                                                  gAGCACCACGGTGTTGAGCGCGGTGCTCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCAAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCATCGCCGGCAAACTGCTGGCCAAGGCGCTg                                                                     >  1:376571/1‑150 (MQ=255)
                                                                                        cacGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATctg                                                               >  2:3739553/1‑148 (MQ=255)
                                                                                               ttGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCg                                                        >  2:1706850/1‑150 (MQ=255)
                                                                                                tGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCgg                                                       >  2:210996/1‑150 (MQ=255)
                                                                                                   gcgcGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCg                                                    >  2:1053159/1‑150 (MQ=255)
                                                                                                    cgcgGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGa                                                   <  2:1327140/150‑1 (MQ=255)
                                                                                                     gcgGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGag                                                  <  1:940266/150‑1 (MQ=255)
                                                                                                     gcgGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGag                                                  >  1:2047289/1‑150 (MQ=255)
                                                                                                       ggTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGAcc                                                                                                                                            <  2:441613/58‑1 (MQ=255)
                                                                                                       ggTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGagca                                                <  2:654242/150‑1 (MQ=255)
                                                                                                         tGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGagcagc                                              <  1:117819/150‑1 (MQ=255)
                                                                                                            gcTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCtt                                           >  1:993107/1‑150 (MQ=255)
                                                                                                             cTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTg                                          >  2:1944659/1‑150 (MQ=255)
                                                                                                              tCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGc                                         <  2:2423269/150‑1 (MQ=255)
                                                                                                              tCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGc                                         >  1:4330896/1‑150 (MQ=255)
                                                                                                                gACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTg                                       >  1:630721/1‑150 (MQ=255)
                                                                                                                gACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTg                                       <  2:325072/150‑1 (MQ=255)
                                                                                                                   ttGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGgcg                                    >  1:1946995/1‑150 (MQ=255)
                                                                                                                      ccTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCaa                                 <  1:1609079/150‑1 (MQ=255)
                                                                                                                      ccTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCaa                                 >  2:201546/1‑150 (MQ=255)
                                                                                                                       cTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAg                                <  2:2484542/150‑1 (MQ=255)
                                                                                                                       cTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAg                                <  1:4049889/150‑1 (MQ=255)
                                                                                                                       cTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAg                                <  2:3415885/150‑1 (MQ=255)
                                                                                                                       cTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAg                                <  2:1682910/150‑1 (MQ=255)
                                                                                                                         gAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGgt                              <  1:2899609/150‑1 (MQ=255)
                                                                                                                               gATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCgat                        >  2:2249608/1‑150 (MQ=255)
                                                                                                                               gATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCgat                        <  2:3952397/150‑1 (MQ=255)
                                                                                                                                 tttCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCgatga                      >  2:2519821/1‑150 (MQ=255)
                                                                                                                                 tggagctggTCAGCGACGGCTCGCCGCTGACCCGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCgatga                      >  1:2784113/8‑150 (MQ=255)
                                                                                                                                   tCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGatgatg                    <  1:2535895/150‑1 (MQ=255)
                                                                                                                                       gggTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATgccgc                <  1:1681051/150‑1 (MQ=255)
                                                                                                                                         gTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGccgccg              <  2:4350207/150‑1 (MQ=255)
                                                                                                                                          tCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCgg             <  1:3993473/150‑1 (MQ=255)
                                                                                                                                            aGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTc           >  1:3379656/1‑150 (MQ=255)
                                                                                                                                            aGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTc           >  2:218957/1‑150 (MQ=255)
                                                                                                                                             gcggcgAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCa          >  2:574224/1‑150 (MQ=255)
                                                                                                                                               ggcgAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGTCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCAcc        <  1:418294/150‑1 (MQ=255)
                                                                                                                                                   aGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGt    <  2:814045/150‑1 (MQ=255)
                                                                                                                                                    gCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCACACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTc   >  2:4063548/1‑150 (MQ=255)
                                                                                                                                                     ccGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCa  <  2:496202/150‑1 (MQ=255)
                                                                                                                                                     ccGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCa  <  2:2751321/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCATCTTCTTGCCGGCGGAAGGCTTGTCCAGGTCTTCCTTGTAGCAATAGGTGCAGCTGAGGTTGCAACCGGTGTTGACGTTGAGCACCACGGTGTTGAGCGCGGTGCGCTCGACTTGCCTGAGGGCGATTTCCGGGGTCAGCGGCGAGCCGTCGCTGACCAGCTCCAGCGCGATCAGTTCGCGCAGGGTGTCCTGGATATCCTCGCCGGCAAACTGCTGGCCAAGGCGCTGCATCAGCTCCTCGGCCGAGCAGCCTTGCTGGCGCAAGGTGTCGATGATGCCGCCGGTCACCGCGTCA  >  NC_002947/3926380‑3926678

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