Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,083,094 T→C 5.9% T148A (ACG→GCG)  purE ← 5‑(carboxyamino)imidazole ribonucleotide mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,083,0940TC 5.9% 242.3 / 3.0 83T148A (ACG→GCG) purE5‑(carboxyamino)imidazole ribonucleotide mutase
Reads supporting (aligned to +/- strand):  ref base T (37/41);  new base C (3/2);  total (40/43)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

AGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACGGCCGATGGCCAGGGTGGCAACCGGCACACCGGCTGGCATCTGCACGATCGA  >  NC_002947/6082948‑6083235
                                                                                                                                                   |                                                                                                                                             
aGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGc                                                                                                                                             <  6:314353/149‑1 (MQ=255)
aGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGc                                                                                                                                             <  8:328426/149‑1 (MQ=255)
  gAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCgg                                                                                                                                           <  3:22475/149‑1 (MQ=255)
     ggCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAAc                                                                                                                                        >  1:275500/1‑149 (MQ=255)
     ggCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAAc                                                                                                                                        >  4:7953/1‑149 (MQ=255)
     ggCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAAc                                                                                                                                        <  3:372571/149‑1 (MQ=255)
       cGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGCACGGAAcag                                                                                                                                      >  8:125762/1‑147 (MQ=255)
        gAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACTCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCAAGCACGGTGTCGGTCTGTTCCGTGCGGAACTg                                                                                                                                      >  1:309119/1‑149 (MQ=255)
        gAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGACCTGc                                                                                                                                     >  5:164925/1‑149 (MQ=255)
          aGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGCACGGAAcagacc                                                                                                                                   <  7:125762/149‑6 (MQ=255)
           gTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTc                                                                                                                                  <  4:356643/149‑1 (MQ=255)
               aTGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAAc                                                                                                                               <  1:11373/148‑1 (MQ=255)
                 gCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGAACAGCTGGACGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGACGGTATGTTACGTGCGGAACTGCTTCAAcgcc                                                                                                                            >  8:179203/1‑149 (MQ=255)
                  cccAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAAcgcc                                                                                                                            <  7:10617/148‑1 (MQ=255)
                   ccAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAAcgccgc                                                                                                                          >  7:251247/1‑149 (MQ=255)
                   ccAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAAcgccgc                                                                                                                          >  2:30367/1‑149 (MQ=255)
                      gCGGGGTTCCCGCCAAGGCCAGGATGCGGCGCAGCTGGCCGCGTCGGATTACACCGATCTTCATTGGTTGAGTCTCATGCCTGGGGCGGGTCTGGGTTGTGCAGCACGGTGTCGGTCTGTTCCGTGGGGAACTGGTTCAACGCCGCGTg                                                                                                                       <  6:179051/149‑1 (MQ=255)
                           gTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGTTTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAAcgccg                                                                                                                           >  4:374513/1‑140 (MQ=255)
                           gTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGTTTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAAcgccg                                                                                                                           <  3:374513/140‑1 (MQ=255)
                              cccGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGc                                                                                                                <  5:97349/148‑1 (MQ=255)
                              cccGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCg                                                                                                               >  7:170215/1‑149 (MQ=255)
                               ccGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCgg                                                                                                              >  6:280485/1‑149 (MQ=255)
                               ccGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCgg                                                                                                              >  4:190597/1‑149 (MQ=255)
                                cGCCAAGCCCAGCATGCGGCCCAGCGGGCCGCCACCGATTACCCCGACCATCCTGCGCTGACCCTCAAGCCGGGCGCGGGCCTGGGTTGTCCAGCACGGTCGCGCTCTGCTCCGTGCGGAACTGCTTCAACGCCGCGTGGCACTGCggg                                                                                                             <  8:170215/149‑1 (MQ=255)
                                   cAAGGCCAGCATGCGGCCCAGCTGGCCGCCAGCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGTGCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTAc                                                                                                          <  7:87263/149‑1 (MQ=255)
                                   cAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAAcgcc                                                                                                                            >  2:367728/1‑131 (MQ=255)
                                   cAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAAcgcc                                                                                                                            <  1:367728/131‑1 (MQ=255)
                                      ggCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTg                                                                                                       >  2:370836/1‑149 (MQ=255)
                                         cAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGgcg                                                                                                    <  5:321387/149‑1 (MQ=255)
                                                cGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGaa                                                                                             <  2:2760/149‑1 (MQ=255)
                                                 ggCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGCGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAAt                                                                                            >  7:244519/1‑149 (MQ=255)
                                                 ggCCCAGCTGGCCGCCACCGATTACACCGAACTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAAt                                                                                            >  1:35080/1‑149 (MQ=255)
                                                 ggCCCAGCTGGCAGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCAGCGTGGTACTGAGGGTACTTGGCGCCCAGAAt                                                                                            >  8:140187/1‑149 (MQ=255)
                                                  gCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATg                                                                                           >  4:113415/1‑149 (MQ=255)
                                                    ccAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCt                                                                                         >  3:203627/1‑149 (MQ=255)
                                                     cAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTg                                                                                        >  8:250138/1‑149 (MQ=255)
                                                      aGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTgg                                                                                       <  1:30367/149‑1 (MQ=255)
                                                       gCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGc                                                                                      >  3:58353/1‑149 (MQ=255)
                                                        cTGGCCGCCACCGCTTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGAACTTGGCGCCCAGAATGCTGGcc                                                                                     <  8:296575/149‑1 (MQ=255)
                                                           gccgccACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGa                                                                                   <  5:111823/148‑1 (MQ=255)
                                                             cgccACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGAcag                                                                                <  3:190597/149‑1 (MQ=255)
                                                                cACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGAcagcag                                                                             <  2:275500/149‑1 (MQ=255)
                                                                    cGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAggg                                                                           >  4:113223/1‑148 (MQ=255)
                                                                    cGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGc                                                                          <  4:176969/149‑1 (MQ=255)
                                                                      aTTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCGGGGCAg                                                                        >  8:24005/1‑149 (MQ=255)
                                                                      aTTACACCGATCTGCATGGGTTGAGCCTCAAGCCTGGCGCGGGAATGGGATGTCCAGCAAGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTAGTGCTGGTACTTGGCGGCCAGAATGCTGGCCGACAGCAGCGCAg                                                                        >  6:57961/1‑149 (MQ=255)
                                                                       ttACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGc                                                                       <  7:201602/149‑1 (MQ=255)
                                                                          cacCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGtt                                                                    <  1:370836/149‑1 (MQ=255)
                                                                          cacCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGt                                                                     >  1:25644/1‑148 (MQ=255)
                                                                          cacCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGt                                                                     <  2:25644/148‑1 (MQ=255)
                                                                           acCGATCTTCATGGGTTGAGCCTCAAGCGTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTg                                                                   >  1:220315/1‑149 (MQ=255)
                                                                           acCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTg                                                                   >  1:306358/1‑149 (MQ=255)
                                                                           acCGATCTTCATGGGTTGAGCCTCAAGACTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTg                                                                   >  6:80434/1‑149 (MQ=255)
                                                                                  ttCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCAc                                                            <  2:309119/149‑1 (MQ=255)
                                                                                    cATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGAc                                                                 <  3:113223/142‑1 (MQ=255)
                                                                                    cATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCa                                                             <  3:253078/146‑1 (MQ=255)
                                                                                    cATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCg                                                          >  8:40838/1‑149 (MQ=255)
                                                                                     aTGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTAAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCgg                                                         >  5:8344/1‑149 (MQ=255)
                                                                                      tGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCgg                                                         >  1:31142/1‑148 (MQ=255)
                                                                                         gTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACg                                                     <  2:31142/149‑1 (MQ=255)
                                                                                               ccTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACGGCCGAt                                               <  3:311375/149‑1 (MQ=255)
                                                                                                cTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACGGCCGATg                                              >  3:38950/1‑149 (MQ=255)
                                                                                                   aaGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACGGCCGATGGcc                                           <  7:140187/149‑1 (MQ=255)
                                                                                                     gCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACGGCCGATGGCCAg                                         >  3:207952/1‑149 (MQ=255)
                                                                                                       cTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCg                                                          >  1:152220/1‑130 (MQ=255)
                                                                                                       cTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCg                                                          <  2:152220/130‑1 (MQ=255)
                                                                                                         ggCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACGGCCGATGGCCAGGGTg                                     >  3:10070/1‑149 (MQ=255)
                                                                                                          gcgcGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACGGCCGATGGCCAGGGTgg                                    <  7:40838/149‑1 (MQ=255)
                                                                                                          gcgcGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACGGCCGATGGCCAGGGTg                                     <  4:10070/148‑1 (MQ=255)
                                                                                                          gcgcGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACGGCCGATGGCCAGGGTg                                     <  3:113415/148‑1 (MQ=255)
                                                                                                                    gggTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACGGCCGATGGCCAGGGTGGCAACCGGca                           >  5:350504/1‑148 (MQ=255)
                                                                                                                          tCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACGGCCGATGGCCAGGGTGGCAACCGGCACACCgg                     <  5:28020/148‑1 (MQ=255)
                                                                                                                          tCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACGGCCGATGGCCAGGGTGGCAACCGGCACACCGGc                    >  4:106611/1‑149 (MQ=255)
                                                                                                                           ccAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACGGCCGATGGCCAGGGTGGCAACCGGCACACCGGCt                   <  3:106611/149‑1 (MQ=255)
                                                                                                                            cAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACGGCCGATGGCCAGGGTGGCAACCGGCACACCGGCTg                  <  6:85314/149‑1 (MQ=255)
                                                                                                                            cAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACGGCCGATGGCCAGGGTGGCAACCGGCACACCGGCTg                  >  7:288554/1‑149 (MQ=255)
                                                                                                                               cACGGTGTCGGTCTGTTTCGCGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACGGCCGATGGCCAGGGTGGCAACCGGCACACCGGCTGGc                <  8:244519/148‑1 (MQ=255)
                                                                                                                                  ggTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACGGCCGATGGCCAGGGTGGCAACCGGCACACCGGCTGGCATCt            <  2:35080/149‑1 (MQ=255)
                                                                                                                                    tgtCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACGGCCGATGGCCAGGGTGGCAACCGGCACACCGGCTGGCATCTGc          <  5:57961/149‑1 (MQ=255)
                                                                                                                                      tCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACGGCCGATGGCCAGGGTGGCAACCGGCACACCGGCTGGCATCTGCa         <  6:164925/148‑1 (MQ=255)
                                                                                                                                           cTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGCCCGCACCGGCACGGCCGATGGCCAGGGTGGCAACCGGCACACCGGCTGGCAt              >  6:4769/1‑138 (MQ=255)
                                                                                                                                           cTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACGGCCGATGGCCAGGGTGGCAACCGGCACACCGGCTGGCAt              <  5:4769/138‑1 (MQ=255)
                                                                                                                                            agTTCCGCACGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACGGCCGATGGCCAGGGTGGCAACCGGCACACCGGCTGGCATCTGCACGATCGa  >  2:349206/2‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
AGGAAGGCGAAGTTCATGCCCAGCGGGGTACCCGCCAAGGCCAGCATGCGGCCCAGCTGGCCGCCACCGATTACACCGATCTTCATGGGTTGAGCCTCAAGCCTGGCGCGGGTCTGGGTTGTCCAGCACGGTGTCGGTCTGTTCCGTGCGGAACTGCTTCAACGCCGCGTGGTACTGCGGGTACTTGGCGCCCAGAATGCTGGCCGACAGCAGGGCAGCGTTGACCGCACCGGCACGGCCGATGGCCAGGGTGGCAACCGGCACACCGGCTGGCATCTGCACGATCGA  >  NC_002947/6082948‑6083235

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: