Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,504,067 T→A 7.5% intergenic (+34/‑229) rpsI → / → petA 30S ribosomal protein S9/ubiquinol‑cytochrome c reductase iron‑sulfur subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,504,0670TA 7.5% 219.9 / 7.7 80intergenic (+34/‑229)rpsI/petA30S ribosomal protein S9/ubiquinol‑cytochrome c reductase iron‑sulfur subunit
Reads supporting (aligned to +/- strand):  ref base T (34/39);  new base A (4/2);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.78e-01

CTGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCATGA  >  NC_002947/1503920‑1504214
                                                                                                                                                   |                                                                                                                                                    
cTGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCtt                                                                                                                                                     >  4:252816/1‑149 (MQ=255)
        cgcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTTAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATTCCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCgg                                                                                                                                             <  2:85747/149‑1 (MQ=255)
        cgcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCgg                                                                                                                                             <  7:314719/149‑1 (MQ=255)
         gcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCggg                                                                                                                                            <  5:378892/149‑1 (MQ=255)
         gcTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCCTGGAACCggg                                                                                                                                            <  8:337007/149‑1 (MQ=255)
          ctctGCGTCGTGCTGGCTACGTCACCCGCGACGATCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTTAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGc                                                                                                                                           <  7:189622/149‑1 (MQ=255)
               cgtcgtGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGTAGTCGAAAAAAGCCCGGTTCCAAGGAACCGGGCtt                                                                                                                                        <  3:388748/146‑1 (MQ=255)
               cgtcgtGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCttttt                                                                                                                                     >  3:59079/1‑149 (MQ=255)
                gtcgtGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCtttttt                                                                                                                                    >  5:102260/1‑149 (MQ=255)
                    tGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCtttttttttt                                                                                                                                <  2:83427/149‑2 (MQ=255)
                      cgggCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGGCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGc                                                                                                                              <  2:362569/147‑1 (MQ=255)
                      cTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCTTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAg                                                                                                                               <  2:16679/148‑1 (MQ=255)
                       tGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGcc                                                                                                                             >  6:211011/1‑149 (MQ=255)
                       tGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACAGGGCTTTTTTTTTAGcc                                                                                                                             >  6:26921/1‑149 (MQ=255)
                         gCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGcc                                                                                                                             <  2:24796/147‑1 (MQ=255)
                         gCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGcc                                                                                                                             >  1:24796/1‑147 (MQ=255)
                             cGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAAc                                                                                                                       <  2:78994/149‑1 (MQ=255)
                             cGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAAc                                                                                                                       <  4:374653/149‑1 (MQ=255)
                             cGTCACCCGCGACGCTCGTGAAGTTGAGCCTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTTAGCCTTg                                                                                                                          <  3:79994/147‑1 (MQ=255)
                                cACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTaa                                                                                                                    >  6:328694/1‑149 (MQ=255)
                                 aCCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCtt                                                                                                                           >  2:392608/1‑141 (MQ=255)
                                 aCCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCtt                                                                                                                           <  1:392608/141‑1 (MQ=255)
                                 aCCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAAc                                                                                                                   <  3:252816/149‑1 (MQ=255)
                                  cccGCGACGCTCGTGAAGTTGAGCGTAATAACGTGTGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAAcc                                                                                                                  <  7:233624/149‑1 (MQ=255)
                                   ccGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTg                                                                                                                          <  4:365977/140‑1 (MQ=255)
                                   ccGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTg                                                                                                                          <  4:371765/140‑1 (MQ=255)
                                   ccGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTg                                                                                                                          >  3:371765/1‑140 (MQ=255)
                                   ccGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTg                                                                                                                          >  3:365977/1‑140 (MQ=255)
                                   ccGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCt                                                                                                                 <  8:180491/149‑1 (MQ=255)
                                   ccGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCt                                                                                                                 >  7:170425/1‑149 (MQ=255)
                                    cgcgACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTg                                                                                                                >  2:121659/1‑149 (MQ=255)
                                    cgcgACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTg                                                                                                                >  2:109708/1‑149 (MQ=255)
                                       gACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCTTAAGCGTCCTCAGTACTCCAAGCGTTAAAACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTca                                                                                                             <  7:349465/149‑1 (MQ=255)
                                                 aaGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTg                                                                                                   >  5:121965/1‑149 (MQ=255)
                                                      gAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGc                                                                                              >  7:121493/1‑149 (MQ=255)
                                                       aGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGc                                                                                              <  4:59079/148‑1 (MQ=255)
                                                           tAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCTAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTc                                                                                         >  6:267161/1‑149 (MQ=255)
                                                           tAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTc                                                                                         >  5:257350/1‑149 (MQ=255)
                                                              gAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCa                                                                                      <  2:376805/149‑1 (MQ=255)
                                                                aaGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCa                                                                                      >  1:44239/1‑147 (MQ=255)
                                                                aaGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATa                                                                                    <  2:241928/149‑1 (MQ=255)
                                                                     ggTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGt                                                                               <  8:27020/149‑1 (MQ=255)
                                                                      gTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACg                                                                                >  8:193010/1‑147 (MQ=255)
                                                                      gTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTa                                                                              <  6:257350/149‑1 (MQ=255)
                                                                          gCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGAt                                                                           >  1:364397/1‑148 (MQ=255)
                                                                                     gTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGt                                                               <  5:328694/149‑1 (MQ=255)
                                                                                       aaGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTgg                                                             >  5:204544/1‑149 (MQ=255)
                                                                                              cTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCg                                                      >  4:114590/1‑149 (MQ=255)
                                                                                                 aGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGccc                                                   <  2:44239/149‑1 (MQ=255)
                                                                                                 aGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGcc                                                    <  2:364397/148‑1 (MQ=255)
                                                                                                   tACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCtt                                                 >  2:197691/1‑149 (MQ=255)
                                                                                                      tCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGcc                                              >  3:379574/1‑149 (MQ=255)
                                                                                                        cAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCgg                                            >  8:199920/1‑149 (MQ=255)
                                                                                                         aaGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCggg                                           <  8:170425/149‑1 (MQ=255)
                                                                                                           gCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCgggg                                          >  4:189038/1‑148 (MQ=255)
                                                                                                              ttAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGAAAAGCGTGGGTTGTCGCCCTTGCCGGGGTaa                                       >  5:87133/1‑148 (MQ=255)
                                                                                                                aaTACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATc                                     <  6:102260/148‑1 (MQ=255)
                                                                                                                 aTACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGAGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCAc                                   >  2:94582/1‑149 (MQ=255)
                                                                                                                    ccGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACAGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCt                                 <  5:26921/148‑1 (MQ=255)
                                                                                                                          cGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTgg                                                                                                  >  2:123985/1‑78 (MQ=255)
                                                                                                                           ggCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAgt                          <  6:121965/148‑1 (MQ=255)
                                                                                                                           ggCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAgt                          <  1:197691/148‑1 (MQ=255)
                                                                                                                           ggCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATAAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGGAATCACCGTGTCAgtg                         >  4:284577/1‑149 (MQ=255)
                                                                                                                            gCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAgtgt                        <  8:121493/149‑1 (MQ=255)
                                                                                                                            gCAGTCGAAAAAAGCCCGGGTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAgtgt                        <  5:211011/149‑1 (MQ=255)
                                                                                                                                tCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCgggg                                          >  8:289928/1‑127 (MQ=255)
                                                                                                                                tCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCgggg                                          <  7:289928/127‑1 (MQ=255)
                                                                                                                                  aaaaaaaGCCCGGTTCCAAGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGGAGATCAAGCAAAGCGTGGGTTGTCGCGCTTGCCGGGGTAATCAGCTTGTCTGTGTGTGGaa                  >  2:205130/2‑148 (MQ=39)
                                                                                                                                      aaaGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGt              >  2:203770/1‑149 (MQ=255)
                                                                                                                                      aaaGCCCGGTTCCAAGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGt              >  3:9020/1‑149 (MQ=255)
                                                                                                                                        aGCCCGGTTCCAAGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAg                                                                                   <  8:33533/78‑1 (MQ=255)
                                                                                                                                        aGCCCGGTTCCAAGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAg                                                                                   >  7:33533/1‑78 (MQ=255)
                                                                                                                                         gCCCGGTTCCAAGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTc           >  2:355698/1‑149 (MQ=255)
                                                                                                                                           ccGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTTGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCtt         <  1:121659/149‑1 (MQ=255)
                                                                                                                                           ccGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCtt         <  1:109708/149‑1 (MQ=255)
                                                                                                                                           ccGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCt          >  5:224734/1‑148 (MQ=255)
                                                                                                                                            cGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTa        >  6:120744/1‑149 (MQ=255)
                                                                                                                                              gTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTAcc      <  7:199920/149‑1 (MQ=255)
                                                                                                                                               ttCCTTGGCACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCa     <  4:37650/149‑1 (MQ=255)
                                                                                                                                                tCCTTGGAACCGGGCTTTTTTTTTAGCATGGAACTAACCTGACGGAAGTCAGTGGCGGATTGACGCATAGAGGTAGAGCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACGTTTACAGTGGGTGGATTTGTTTCttccaat    >  4:306046/1‑144 (MQ=255)
                                                                                                                                                  cTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCATGa  <  5:267161/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTGCGTGGCGCTCTGCGTCGTGCTGGCTACGTCACCCGCGACGCTCGTGAAGTTGAGCGTAAGAAAGTGGGTCTGCGTAAAGCGCGTAAGCGTCCTCAGTACTCCAAGCGTTAATACCGCTTCGGCAGTCGAAAAAAGCCCGGTTCCTTGGAACCGGGCTTTTTTTTTAGCCTTGAACTAACCTGTCAGCAGTCAGTGGCCGCTTGTCGCATAGACGTAGATCAAGCAAAGCGTGGGTTGTCGCCCTTGCCGGGGTAATCACCTTGTCAGTGTGTGGATTTGTTTCTTACCATGA  >  NC_002947/1503920‑1504214

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