Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,879,933 C→A 5.6% intergenic (+75/+61) PP_1686 → / ← PP_1687 glutathione peroxidase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,879,9330CA 5.6% 689.1 / 19.2 206intergenic (+75/+61)PP_1686/PP_1687glutathione peroxidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (91/103);  new base A (7/5);  total (98/108)
Fisher's exact test for biased strand distribution p-value = 5.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 213 total aligned reads displayed.

TTCTCCAGCCTGACCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCACCGGCATCGCCAAC  >  NC_002947/1879784‑1880081
                                                                                                                                                     |                                                                                                                                                    
ttCTCCAGCCTGACCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGc                                                                                                                                                      >  1:546451/1‑150 (MQ=255)
ttCTCCAGCCTGACCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGc                                                                                                                                                      >  2:4714915/1‑150 (MQ=255)
ttCTCCAGCCTGACCAAACCGGATGATCCCGAGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGAGACTTCAGGAGGc                                                                                                                                                      >  1:3214868/1‑150 (MQ=255)
      aGCCTGACCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTtct                                                                                                                                                >  2:3443218/1‑150 (MQ=255)
       gCCTGACCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTtctc                                                                                                                                               <  2:1887882/150‑1 (MQ=255)
       gCCTGACCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTtctc                                                                                                                                               <  1:4145796/150‑1 (MQ=255)
           gACCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGa                                                                                                                                           <  2:4194345/150‑1 (MQ=255)
            aCCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGAACGGCCTCCTGaa                                                                                                                                          <  2:1605861/150‑1 (MQ=255)
             ccAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAg                                                                                                                                         <  1:3488839/150‑1 (MQ=255)
              cAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGt                                                                                                                                        >  2:750940/1‑150 (MQ=255)
              cAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGAACGGCCTCCTGAAGt                                                                                                                                        <  1:2814280/150‑1 (MQ=255)
               aaaCCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGtt                                                                                                                                       >  1:440251/1‑150 (MQ=255)
               aaaCCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCAACTTCAGGAGAACGGCCTCCTGAAGtc                                                                                                                                       <  1:4037988/150‑2 (MQ=39)
                aaCCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCCGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGCTTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTg                                                                                                                                      >  1:4829133/1‑150 (MQ=255)
                  ccGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCg                                                                                                                                    <  2:546451/150‑1 (MQ=255)
                   cGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCgg                                                                                                                                   <  1:4542180/150‑1 (MQ=255)
                     gatgatCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGga                                                                                                                                 >  2:3255224/1‑150 (MQ=255)
                       tgatCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatc                                                                                                                               >  2:4716945/1‑150 (MQ=255)
                       tgatCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatc                                                                                                                               >  1:241332/1‑150 (MQ=255)
                       tgatCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatc                                                                                                                               <  1:1173516/150‑1 (MQ=255)
                          tCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatcaag                                                                                                                            >  1:3442077/1‑150 (MQ=255)
                          tCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatcaag                                                                                                                            >  1:3441972/1‑150 (MQ=255)
                           cccGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatcaaga                                                                                                                           <  2:1088359/150‑1 (MQ=255)
                           cccGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatcaaga                                                                                                                           >  1:1088359/1‑150 (MQ=255)
                             cgcgTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCgg                                                                                                                                   <  1:3321602/140‑1 (MQ=255)
                             cgcgTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatcaagatc                                                                                                                         >  2:1862112/1‑150 (MQ=255)
                              gcgTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatcaagatca                                                                                                                        <  2:804925/150‑1 (MQ=255)
                              gcgTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatcaagatca                                                                                                                        <  2:3689328/150‑1 (MQ=255)
                              gcgTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGgatcaagatca                                                                                                                        <  2:1242942/150‑1 (MQ=255)
                                    cgcgCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGcaca                                                                                                                  <  2:4253032/150‑1 (MQ=255)
                                        ctgctATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGaa                                                                                                              <  2:820084/150‑1 (MQ=255)
                                           ctATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGtg                                                                                                           <  1:1946273/150‑1 (MQ=255)
                                            tATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCATGCACAGGAAGtgg                                                                                                          <  2:2121709/150‑1 (MQ=255)
                                            tATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGtgg                                                                                                          >  2:4400985/1‑150 (MQ=255)
                                            tATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGtgg                                                                                                          >  2:1406850/1‑150 (MQ=255)
                                            tATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGtgg                                                                                                          <  1:1862112/150‑1 (MQ=255)
                                            tATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGAGCAAGATCACGCACAGGACGtgg                                                                                                          >  2:4214916/1‑150 (MQ=255)
                                             aTCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGtggc                                                                                                         >  1:692737/1‑150 (MQ=255)
                                                gAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGt                                                                                                                                                   >  1:1911276/1‑105 (MQ=255)
                                                 aaaaGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGtggctggc                                                                                                     >  1:3009539/1‑150 (MQ=255)
                                                     ggCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCg                                                                                                 <  2:2736967/150‑1 (MQ=255)
                                                      gCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCgg                                                                                                <  2:4283241/150‑1 (MQ=255)
                                                               tCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCa                                                                                       >  1:1411578/1‑150 (MQ=255)
                                                               tCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGAACGGCCTCCTGAAg                                                                                                                                         <  1:1515507/100‑1 (MQ=37)
                                                                   aGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCggg                                                                                                                                  <  1:1018010/103‑1 (MQ=255)
                                                                   aGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCggg                                                                                                                                  >  2:1018010/1‑103 (MQ=255)
                                                                   aGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGt                                                                                   >  2:1276139/1‑150 (MQ=255)
                                                                   aGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGt                                                                                   <  1:4883780/150‑1 (MQ=255)
                                                                    gCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGtg                                                                                  <  1:89915/150‑1 (MQ=255)
                                                                     ccGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGtgc                                                                                 <  1:3721656/150‑1 (MQ=255)
                                                                       gTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGtgcat                                                                               <  2:3759220/150‑1 (MQ=255)
                                                                           ccAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCAc                                                                           <  2:2174158/150‑1 (MQ=255)
                                                                           ccAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCAc                                                                           <  2:4235064/150‑1 (MQ=255)
                                                                           ccAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCAc                                                                           <  1:2654683/150‑1 (MQ=255)
                                                                            cAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACg                                                                          <  2:4672412/150‑1 (MQ=255)
                                                                              ggCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGt                                                                        >  2:5204098/1‑150 (MQ=255)
                                                                               gCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTa                                                                       >  2:1579528/1‑150 (MQ=255)
                                                                                 tGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTTACCGTGCATGCACGGtgaa                                                                     <  2:1650152/150‑4 (MQ=255)
                                                                                  gcgcTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATa                                                                    <  1:1120132/150‑1 (MQ=255)
                                                                                     cTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGcc                                                                 >  2:2383520/1‑150 (MQ=255)
                                                                                          tgcGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGtgc                                                                                 <  2:4932559/129‑1 (MQ=255)
                                                                                            cggcgTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTaaagt                                                                   <  2:176183/141‑4 (MQ=255)
                                                                                            cggcgTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTaaagt                                                                   >  1:176183/1‑138 (MQ=255)
                                                                                                gTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGcac                                                      <  1:2740688/150‑1 (MQ=255)
                                                                                                     tttAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGtggctgg                                                                                                      >  2:666441/1‑97 (MQ=255)
                                                                                                     tttAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGcaccacca                                                 >  2:3274277/1‑150 (MQ=255)
                                                                                                      ttAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGcaccaccac                                                >  1:3188440/1‑150 (MQ=255)
                                                                                                          aGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACtt                                                                                         >  2:1112758/1‑105 (MQ=255)
                                                                                                          aGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCAccggc                                            >  2:4061078/1‑150 (MQ=255)
                                                                                                              cgAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCAcc                                                                                     >  2:5269536/1‑105 (MQ=255)
                                                                                                              cgAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGa                                        <  2:4002780/150‑1 (MQ=255)
                                                                                                               gAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCAc                                                                           <  1:1847593/114‑1 (MQ=255)
                                                                                                               gAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCAc                                                                           >  2:1847593/1‑114 (MQ=255)
                                                                                                                 ggCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGacca                                     <  2:2429770/150‑1 (MQ=255)
                                                                                                                    aGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGAccacca                                  <  1:4714915/150‑1 (MQ=255)
                                                                                                                    aGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGAccacca                                  >  1:1260741/1‑150 (MQ=255)
                                                                                                                    aGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGAccacca                                  <  2:241332/150‑1 (MQ=255)
                                                                                                                     gTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAg                                 >  1:1313632/1‑150 (MQ=255)
                                                                                                                           gCGACGCAGATTGCGACTTCAGGATGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGcacaca                           <  2:3600/150‑1 (MQ=255)
                                                                                                                           gCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGcacaca                           >  2:3181694/1‑150 (MQ=255)
                                                                                                                           gCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGcacaca                           <  1:293564/150‑1 (MQ=255)
                                                                                                                            cGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGcacacac                          <  2:1994206/150‑1 (MQ=255)
                                                                                                                             gACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGcacacaca                         >  1:5245290/1‑150 (MQ=255)
                                                                                                                              aCGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGcacacacac                        <  2:4505656/150‑1 (MQ=255)
                                                                                                                               cGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGcacacacaca                       >  1:4319764/1‑150 (MQ=255)
                                                                                                                                gCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACaa                      >  1:2020628/1‑150 (MQ=255)
                                                                                                                                gCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACaa                      <  1:3274277/150‑1 (MQ=255)
                                                                                                                                  aGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAAt                    <  1:4950370/150‑1 (MQ=255)
                                                                                                                                     ttGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCGCCGGCCGGACCACCAGCACACACACAAATGcc                 >  1:624691/1‑150 (MQ=255)
                                                                                                                                      tGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCa                <  1:1680752/150‑1 (MQ=255)
                                                                                                                                      cgcaaCTTCAGGAGAACGGCCTCCTGAAGTCGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCa                >  1:1785438/5‑150 (MQ=255)
                                                                                                                                      cgcaaCTTCAGGAGAACGGCCTCCTGAAGTCGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCa                >  1:254317/5‑150 (MQ=255)
                                                                                                                                       gCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCAc               >  1:2365208/1‑150 (MQ=255)
                                                                                                                                        cGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATCCCAcc              <  2:4326827/150‑1 (MQ=255)
                                                                                                                                         gACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAACTGCCACCg             <  2:4829133/150‑1 (MQ=255)
                                                                                                                                         gACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCACCg             >  1:3628128/1‑150 (MQ=255)
                                                                                                                                         gACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCACCg             <  1:4068256/150‑1 (MQ=255)
                                                                                                                                              cAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCACCGGCATc        <  2:2497998/150‑1 (MQ=255)
                                                                                                                                                 gaggCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCACCGGCATCGcc     <  1:5297680/150‑1 (MQ=255)
                                                                                                                                                    gCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCACCGGCATCGCCAAc  >  2:4404244/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTCTCCAGCCTGACCAAACCGGATGATCCCGCGTTCCGCGCTGCTATCGAAAAGGCCATTGCCTCACAGCCGTGACCAGGCTGCGCTTGCTGCGGCGTTGCTTTAAAGCGCGAGGCAGTTCAGGCGACGCAGATTGCGACTTCAGGAGGCCGTTCTCCTGAAGTTGCGGGATCAAGATCAAGCACAGGAAGTGGCTGGCTGCGGGGCACTTCACCGTGCATGCACGGTAATAGCCACCAACCGCACCACCACCGGCCGGACCACCAGCACACACACAAATGCCACCGGCATCGCCAAC  >  NC_002947/1879784‑1880081

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: