Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,316,713 T→G 10.3% intergenic (‑226/+52) PP_1144 ← / ← rapA membrane protein/RNA polymerase‑associated protein RapA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,316,7130TG10.3% 1281.6 / 12.2 460intergenic (‑226/+52)PP_1144/rapAmembrane protein/RNA polymerase‑associated protein RapA
Reads supporting (aligned to +/- strand):  ref base T (151/223);  new base G (49/12);  total (219/241)
Fisher's exact test for biased strand distribution p-value = 4.55e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 473 total aligned reads displayed.

TACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGTCGATTT  >  NC_002947/1316566‑1316862
                                                                                                                                                     |                                                                                                                                                     
tactgctgCCTTCGTGACTACTGGCACAATGTTTTCGCCGTTTTGGTTCGCACTCCAATCGCAGTGAGGAATACATTTGGGGCCCTACTGCCTGCCGGTGGGAGCGGGCGCGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCttt                                                                                                                                                     <  1:4642774/143‑1 (MQ=255)
  cTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGc                                                                                                                                                   <  2:1346099/150‑1 (MQ=255)
   ttAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGGAGCAAACACAGCGGTGGATGCCACCCGCGTTGcc                                                                                                                                                  >  1:3543149/1‑150 (MQ=255)
    tAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCg                                                                                                                                                 <  2:3412233/150‑1 (MQ=255)
    tAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCg                                                                                                                                                 <  2:6411413/150‑1 (MQ=255)
     aaGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGGCCGCGAAGCATGCACTTCGGTGGTTTGCACCGGCTTTGGCgg                                                                                                                                                >  2:6733404/1‑150 (MQ=255)
     aaGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCgg                                                                                                                                                <  2:6929992/150‑1 (MQ=255)
       gCCTTCCTGACTAATGGTACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCAAAGCCGgtt                                                                                                                                              <  2:3688191/150‑2 (MQ=255)
       gCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGTAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGCGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGgtg                                                                                                                                              <  1:4024563/150‑1 (MQ=255)
       gCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTCAGGCCTTCCGGTGGGCGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGgtg                                                                                                                                              <  1:4231099/150‑1 (MQ=255)
        ccttccttACTAATGGCACAATGATTTCTACGATTTGGTTATCACTCCATTCGCTTTGAGTAATAATTTTGGGTCCCTAATGCCTTCCTGTGGTAGCGGGCGAGCCCGCGAAGCATCCACTGCGTTGGATGCCACCGGCTTTGCCGgtgt                                                                                                                                             <  1:1007003/150‑1 (MQ=38)
        ccTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGgtgt                                                                                                                                             <  2:1113794/150‑1 (MQ=255)
        ccTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGGGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCAAAGCCGgtgg                                                                                                                                             <  1:7544077/150‑2 (MQ=255)
         cttcttGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCAAAGCCGgtggc                                                                                                                                            <  2:7516627/150‑3 (MQ=255)
          ttCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATTACTCCCATCGCAGTTAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCTCGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTc                                                                                                                                           <  2:8744893/150‑1 (MQ=38)
          ttCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTc                                                                                                                                           <  2:7230750/150‑1 (MQ=255)
          ttCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACAGCGGTGGATACCACCCGCGTTGCCGGGGTTc                                                                                                                                           >  2:3884346/1‑150 (MQ=255)
          ttCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCAACCACAGCGGTGGGTGCCACCCGCGTTGCCGGTGTTc                                                                                                                                           >  2:3525094/1‑150 (MQ=255)
           tCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGGCCGCGTATCAAACAATGCGGTGGTTTGCACCGGCTTTTGCCGTGttgg                                                                                                                                          >  2:3064529/1‑148 (MQ=255)
              tGACTAATGGCACAATGATTTCTACGTTTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCgg                                                                                                                                       <  1:6682193/150‑1 (MQ=255)
              tGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGGCCGCGGAGCAACCACAGGCGGGGGTGGCACCCGCGTTGGCGGTGTTCGCgg                                                                                                                                       >  1:5703780/1‑150 (MQ=255)
              tGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCgg                                                                                                                                       <  2:5274601/150‑1 (MQ=255)
                  tAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGCAGCAAACACAACGGGGGGTGCCCACCCCGTTTGCCGGGTTCGCgggggg                                                                                                                                   >  1:7794672/1‑147 (MQ=255)
                    aTGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgc                                                                                                                                 >  2:7282154‑M2/1‑148 (MQ=255)
                    aTGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCATTGCCGGTGTTAGCGGGCGTgc                                                                                                                                 <  2:4300127‑M2/150‑3 (MQ=255)
                      ggCACAATGATTTCTACGTTTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgccc                                                                                                                               >  1:1486265‑M2/1‑146 (MQ=255)
                       gCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGGCCGCCCACCAAACACAGCGGGGGGTGGCACCCGCGTTTGCCGGGTTGGGGGGCggggccc                                                                                                                              >  1:1899888/1‑144 (MQ=255)
                         acaATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCGCACCAACCACAGCGGTGGGTTCCCCCCGCGTTTCCGGTGTTAGCGGGCGTgccccca                                                                                                                            <  2:3222912‑M2/150‑8 (MQ=255)
                         acaATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCGAAGCAAACAAAGCGGTGGATTCCACCCGCATTGCCCGTGTTCGCCGGCCGGccccca                                                                                                                            >  1:8815625/1‑147 (MQ=255)
                           aaTGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCCAAACAACCACAGCGGGGGCTGCCACCCGCGTTTCCGGTGTTCGCGGGCGGgcccgcacc                                                                                                                          <  2:2936529‑M2/150‑10 (MQ=255)
                             tGATTTCTACGTTTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctccca                                                                                                                        >  1:1319074‑M2/1‑139 (MQ=255)
                             tGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGGCCCCGCATCAACCACAGCGGTGGTTGGCACCCGCGTTTGCGGGGTTTGCGGGGGGGGCCGCGCCca                                                                                                                        >  2:2264931/1‑150 (MQ=255)
                                    tACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCCTCCACTGCGGGGGGAGCCACCCGCATTTCCGGGGTTGGCGGGGGGGGCCCCACCCACAAgggg                                                                                                                 >  1:2136236/1‑150 (MQ=255)
                                         ttgggTTATAACGCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCgg                                                                                                            <  2:8662755/147‑1 (MQ=255)
                                            ggTTATACCTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGc                                                                                                         <  2:7759930/150‑1 (MQ=255)
                                             gTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATACACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTAGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGcc                                                                                                        >  2:7066622/1‑150 (MQ=255)
                                                  aaCTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCGTTGCCGGTGTTCGCGGGCGTGCCCGCACCCACCAGGGGCGCGGGGCCTgcgc                                                                                                   >  2:5008112/1‑150 (MQ=255)
                                                     tCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGGCCCCCCAACAAACAAAGCGGTGGTTGGCCCCCGCCGTTCCGGGGTTGGGGGGCGGGGCCCCTCCCACCAAGGGGGCGGGGGCGGCGCGag                                                                                                >  2:1001572/1‑150 (MQ=255)
                                                     tCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCCCGAAGCAAACACAGCGGTGGATGCCACCCGCGTTGCCCGTGTTCGCGGGCGTGCCCCCACCCACAAAGGGCCCGGGGCCTGCGCTag                                                                                                >  2:1974495/1‑150 (MQ=255)
                                                     tCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGAGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTag                                                                                                >  2:7846282/1‑150 (MQ=255)
                                                      ccACTCGCAGTGAGTAATAAATTTGGGGCCCTAATGCCTTCCGTTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTagc                                                                                               <  1:2027307/150‑1 (MQ=255)
                                                       cAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCGGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGAGCGCAGGCCCCgcgcccctt                                                                                              <  2:8313335/150‑6 (MQ=37)
                                                        aaTCGCAGTGTGTAATAATTTTGGGGCCCTAATGCCTTCGTGTGGGAGCGGGCGCGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTagcca                                                                                             <  1:6209695/150‑1 (MQ=37)
                                                        aaTCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTagcca                                                                                             <  1:6559621/150‑1 (MQ=255)
                                                             cAGTGAGGAATAAATTTGGGGCCCTACTGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGccac                                                                                        <  1:577769/150‑1 (MQ=255)
                                                                         aaTTGGGGGCCCTAATGCCTTCCTGGGGGCGGGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAAt                                                                            <  1:3751506/150‑1 (MQ=255)
                                                                            ttGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGc                                                                         <  2:6582133/150‑1 (MQ=255)
                                                                              ggggCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGCAGCATCCCCAGCGGTGGCTGCCACCCGCTTTGCCGGGGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCGAGCCAGCCACCAGCCCCCGCAAGGCtt                                                                       >  2:7813858/1‑150 (MQ=37)
                                                                                 gCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCa                                                                    <  2:1569724/150‑1 (MQ=255)
                                                                                    cTACTGCCTTCCTGTGGGCGCGGGCGCGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGgc                                                                 <  1:2518548/150‑1 (MQ=255)
                                                                                       aggCCTTCCGGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCa                                                              <  1:7108809/148‑1 (MQ=255)
                                                                                         gCCTTCCTGTGGGAGCGGGCGAGGCCGGGAAGCATACAAAGCGGTGGATTGCACCCGCGTTTGCGGTGTTGGCGGGCGTGGCCGCGCCCCCCACGGGCGCGGGGGCGGCGCGAGCCCGCCCCCACCACCCCCATGGCGTCCCGGCGCacg                                                            >  2:8207295/1‑148 (MQ=255)
                                                                                           cTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggc                                                          <  2:4416531/150‑1 (MQ=255)
                                                                                            ttCGTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcg                                                         <  2:5422746/150‑1 (MQ=255)
                                                                                              cgtgtgGGCGGGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGGGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcggc                                                       <  1:6889933/148‑1 (MQ=255)
                                                                                               cTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGcc                                                      <  2:2635226/150‑1 (MQ=255)
                                                                                                    ggAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaggaaggcattagggccccaaatttattcctcactgcgattggagttataaccaaatcgtagaaatcatt                                                 <  1:8674817‑M2/150‑83 (MQ=255)
                                                                                                    ggAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTtctc                                                 <  2:3793335/150‑1 (MQ=255)
                                                                                                     gAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTcc                                                <  2:1684257/150‑1 (MQ=255)
                                                                                                       gCGGGCGAGCCCGCGAAGCATACACTGCGGTGGATTCCACCGGCTTTTCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAg                                              >  2:7960936/1‑150 (MQ=255)
                                                                                                          ggCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccccaggaaggccttagggcccccaaggtattcctcacggcgattggagtgataacccaatcggagaaatcattgtgcca                                           <  2:4855137‑M2/150‑89 (MQ=255)
                                                                                                           gCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCTCgcccgcgcccacaggaaggcattagggccccaaaattattactaacagcgattggggttataaccaaaaaggagaaataattgttccat                                          <  1:3219919‑M2/150‑90 (MQ=255)
                                                                                                           gCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaggaaggcattagggccccaaatttattcctcactgcgattggagttataaccaaatcgtagaaatcattgtgccat                                          <  1:3619589‑M2/150‑90 (MQ=255)
                                                                                                           gCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaggaaggcattagggccccaaatttattcctcactgcgattggagttataaccaaatcgtagaaatcattgtgccat                                          <  2:5669513‑M2/150‑90 (MQ=255)
                                                                                                           gCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaggaaggcattagg                                                                                                         <  1:529434‑M2/87‑27 (MQ=255)
                                                                                                           gCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCAtt                                          <  2:6107585/150‑1 (MQ=255)
                                                                                                           gCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCAtt                                          <  1:5917829/150‑1 (MQ=255)
                                                                                                             gAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGc                                        <  2:1817656/150‑1 (MQ=255)
                                                                                                             gAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGc                                        >  1:2599290/1‑150 (MQ=255)
                                                                                                                cccGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAg                                     <  2:2540383/150‑1 (MQ=255)
                                                                                                                cccGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAg                                     <  1:8024208/150‑1 (MQ=255)
                                                                                                                 ccGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTgcccgctcccacaagaaggcattagggccccaaatttattcctcactgcgattggagttataaccaaatcgtagaaatcattgtgccattagtca                                    <  1:2778509‑M2/150‑96 (MQ=255)
                                                                                                                 ccGCGAAGCATCCACTGCGGTGGATGCCACCGGCGTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAgg                                    >  2:4191186/1‑150 (MQ=255)
                                                                                                                 ccGCGAAACATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCTCgcccgctcccacaggaaggcagtagggccccaaatttattcctcacggcgattggagttataaacaaatcgtagaaatcattgtgccattagtca                                    <  1:2325135‑M2/150‑96 (MQ=255)
                                                                                                                   gcgAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGcc                                  >  1:7372909/1‑150 (MQ=255)
                                                                                                                   cggtggCATCCACCGCAGTGGATGCTACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGcc                                  >  1:6592209/6‑150 (MQ=255)
                                                                                                                     gAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCAt                                <  2:3933708/150‑1 (MQ=255)
                                                                                                                      tggCATCCACCGCAGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTagcc                                                                                              <  2:2515474/85‑1 (MQ=255)
                                                                                                                      tggCATCCACCGCAGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATc                               >  2:5044772/3‑150 (MQ=255)
                                                                                                                        gCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCtt                             <  1:5053279/150‑1 (MQ=255)
                                                                                                                         cATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTc                            >  2:2726254/1‑150 (MQ=255)
                                                                                                                          aTCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAgg                                                                  <  2:2803377/111‑1 (MQ=255)
                                                                                                                          aTCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTcg                           >  2:2075685/1‑150 (MQ=255)
                                                                                                                           tCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTcgc                          >  1:1088880/1‑150 (MQ=255)
                                                                                                                           tCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTcgc                          <  2:5844313/150‑1 (MQ=255)
                                                                                                                           tCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTcgc                          >  2:5539488/1‑150 (MQ=255)
                                                                                                                              aCTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGcgcc                       >  1:48673/1‑150 (MQ=255)
                                                                                                                                tGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGcgccgc                     >  2:8080538/1‑150 (MQ=255)
                                                                                                                                   ggTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTg                  <  1:8719520/150‑1 (MQ=255)
                                                                                                                                       gATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCa              <  1:3024830/150‑1 (MQ=255)
                                                                                                                                         tGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAgg            <  1:8150766/150‑1 (MQ=255)
                                                                                                                                          gCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAggg           >  1:4988491/1‑150 (MQ=255)
                                                                                                                                          gCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAggg           >  2:1504674/1‑150 (MQ=255)
                                                                                                                                           aacACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCAGTGCCAGGCCATCTACGCGCCGCTTGCGCAggg           >  2:4630438/3‑150 (MQ=255)
                                                                                                                                              aCCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGt        >  2:1181739/1‑150 (MQ=255)
                                                                                                                                               ccGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGTc       >  2:6250388/1‑150 (MQ=255)
                                                                                                                                               ccGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGTc       <  1:1666080/150‑1 (MQ=255)
                                                                                                                                               ccGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGTc       <  2:3611415/150‑1 (MQ=255)
                                                                                                                                                 ggCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGTCGa     <  1:4300127/150‑1 (MQ=255)
                                                                                                                                                     tttGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGTCGAttt  >  1:7241233/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGCTTGCGCAGGGCGTCGATTT  >  NC_002947/1316566‑1316862

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