Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 797,982 T→G 11.3% intergenic (+40/‑46) PP_0683 → / → fklB‑I hypothetical protein/FKBP‑type peptidyl‑prolyl cis‑trans isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947797,9820TG11.3% 174.8 / 13.9 62intergenic (+40/‑46)PP_0683/fklB‑Ihypothetical protein/FKBP‑type peptidyl‑prolyl cis‑trans isomerase
Reads supporting (aligned to +/- strand):  ref base T (28/27);  new base G (5/2);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

GTCTCGATGAACCTGGTGCACTGGTCGGTGCAGGACGAAATGTGCGGCATGGGTACGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTGGGTGGCCAGCTGCGCGACAACCCGCCACCAGGCGTTA  >  NC_002947/797835‑798127
                                                                                                                                                   |                                                                                                                                                 
gTCTCGATGAACCTGGTGCACTGGTCGGTGCAGGACGAAATGTGCGGCATGGGTACGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCtttt                                                                                                                                                  >  5:106333/1‑149 (MQ=255)
gTCTCGATGAACCTGGTGCACTGGTCGGTGCAGGACGAAATGTGCGGCATGGGTACGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCttt                                                                                                                                                   <  3:101047/148‑1 (MQ=255)
 tctcGATGAACCTGGTGCACTGGTCGGTGCAGGACGAAATGTGCGGCATGGGTACGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCtttt                                                                                                                                                  >  8:187647/1‑148 (MQ=255)
     gATGAACCTGGTGCACTGGTCGGTGCAGGACGAAATGTGCGGCATGGGTACGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATc                                                                                                                                              >  7:65491/1‑148 (MQ=255)
     gATGAACCTGGTGCACTGGTCGGTGCAGGACGAAATGTGCGGCATGGGTACGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTGTTTATc                                                                                                                                              >  8:154026/1‑148 (MQ=255)
      aTGAACCTGGTGCACTGTTCGGTTCAGGACGAAATGTGCGGCATGGGTACGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCaa                                                                                                                                            <  8:188323/149‑1 (MQ=255)
      aTGAACCTGGTGCACTGGTCGGTGCAGGACGAAATGTGCGGCATGGGTACGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATc                                                                                                                                              >  7:16945/1‑147 (MQ=255)
      aTGAACCTGGTGCACTGGTCGGTGCAGGACGAAATGTGCGGCATGGGTACGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATc                                                                                                                                              >  6:178685/1‑147 (MQ=255)
           ccTGGTGCACTGGTCGGTGCAGGACGAAATGTGCGGCATGGGTACGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGgcg                                                                                                                                       <  8:246520/149‑1 (MQ=255)
           ccTGGTGCACTGGTCGGTGCAGGACGAAATGTGCGGCATGGGTACGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGgcg                                                                                                                                       >  5:238895/1‑149 (MQ=255)
           ccTGGTGCACTGGTCAGTGCAGGACGAAATGTGCGGCATGGGTACGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGgcg                                                                                                                                       >  6:46559/1‑149 (MQ=255)
               gTGCACTGGTCGGTGCAGGACGAAATGTGCGGCATGGGTACGCACCCTGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTGATAAAAAGGgcgagca                                                                                                                                   <  1:159190/149‑5 (MQ=255)
               gTGCACTGGTCGGTGCAGGACGAAATGTGCGGCATGGGTACGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCgg                                                                                                                                    <  8:16945/148‑1 (MQ=255)
               gTGCACTGGTCGGTGCAGGACGAAATGTGCGGCATGGGTACGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGc                                                                                                                                   >  4:254847/1‑149 (MQ=255)
                         cGGTGCAGGACGAAATGTGCGGCATGGGTCCGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCTGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTgcgc                                                                                                                         <  4:245516/149‑1 (MQ=255)
                          ggTGCAGGACGAAATGTGCGGCATGGGTACGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGcc                                                                                                                        <  4:46380/149‑1 (MQ=255)
                          ggTGCAGGACGAAATGTGCGGCATGGGTACGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGcc                                                                                                                        <  4:101248/149‑1 (MQ=255)
                             gCAGGACGAAATGTGCGGCATGGGTACGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTg                                                                                                                     <  6:106333/149‑1 (MQ=255)
                                 gACGAAATGTGCGGCATGGGTACGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGt                                                                                                                 >  6:185350/1‑149 (MQ=255)
                                     aaaTGTGCGGCATGGGTACGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAg                                                                                                             <  8:65491/149‑1 (MQ=255)
                                          tgCGGCATGGGTACGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGa                                                                                                        >  3:13234/1‑149 (MQ=255)
                                          tgCGGCATGGGTACGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGa                                                                                                        <  5:178685/149‑1 (MQ=255)
                                                     tACGCACCCGGAAGCAGACCTGAAGGGGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGaa                                                                                             <  8:237426/149‑1 (MQ=255)
                                                          aCCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCaa                                                                                        >  4:288112/1‑149 (MQ=255)
                                                            ccGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAAcc                                                                                      >  5:148566/1‑149 (MQ=255)
                                                                   aagaCCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTAc                                                                               <  3:11270/149‑1 (MQ=255)
                                                                    agaCCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTAcc                                                                              <  5:185350/149‑1 (MQ=255)
                                                                        cTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACcgacg                                                                          <  7:187647/149‑1 (MQ=255)
                                                                              gTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAAcgc                                                                    >  7:269808/1‑149 (MQ=255)
                                                                                 ccGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTc                                                                                          >  6:45184/1‑124 (MQ=255)
                                                                                 ccGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTc                                                                                          <  5:45184/124‑1 (MQ=255)
                                                                                  cGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAAcgcgcg                                                                 <  3:91839/148‑1 (MQ=255)
                                                                                  cGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAAcgcgcg                                                                 >  4:91839/1‑148 (MQ=255)
                                                                                  cGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGt                                                                >  8:5262/1‑149 (MQ=255)
                                                                                  cGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGt                                                                <  7:5262/149‑1 (MQ=255)
                                                                                          cAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTCCCGACGACACGCGCGTCAGCTACg                                                        >  8:106854/1‑149 (MQ=255)
                                                                                                 cGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCgg                                                 >  2:174432/1‑149 (MQ=255)
                                                                                                     cTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAAcgcgc                                                                  >  5:305608/1‑128 (MQ=255)
                                                                                                     cTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAAcgcgc                                                                  <  6:305608/128‑1 (MQ=255)
                                                                                                      tGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGt                                             >  4:112956/1‑148 (MQ=255)
                                                                                                             tgctCGCCGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGTCCGCAACCTGTCTACCGCCGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTggg                                     >  3:48132/1‑149 (MQ=255)
                                                                                                             tgctCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTgg                                      >  8:173192/1‑148 (MQ=255)
                                                                                                               ctCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTgg                                      <  7:173192/146‑1 (MQ=255)
                                                                                                               ctCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTGGGt                                    <  6:206220/148‑1 (MQ=255)
                                                                                                                  gCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTc                                                               <  4:267492/118‑1 (MQ=255)
                                                                                                                  gCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTc                                                               >  3:267492/1‑118 (MQ=255)
                                                                                                                  gCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTGGGTGGc                                 >  2:74002/1‑148 (MQ=255)
                                                                                                                         cccAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCAGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTGGGTGGCCAGCTgcg                         >  5:129319/1‑149 (MQ=255)
                                                                                                                           cAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTGGGTGGCCAGCTgcgc                        >  1:119402/1‑148 (MQ=255)
                                                                                                                           cAGTGTGCTAGGCTGCTCGCCCTTTTTAACTAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTGGGTGGCCAGCTgcgc                        >  1:183811/1‑148 (MQ=255)
                                                                                                                            aGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTGGGTGGCCAGCTGCGCGa                      <  8:301016/149‑1 (MQ=255)
                                                                                                                                tgCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCATCTACGGTATCGGCCGTCAGCTGGGTGGCCAGCTGCGCGACAAc                  <  5:46559/149‑1 (MQ=255)
                                                                                                                                  cTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTGGGTGGCCAGCTGCGCGACAAcc                 <  1:74002/148‑1 (MQ=255)
                                                                                                                                     ggCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTGGGTGGCCAGCTGCGCGACAACCCGc              >  5:122132/1‑148 (MQ=255)
                                                                                                                                        gccgcGCCCTTGATAAAAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTGGGTGGCCAGCTGCGCGACAACCCGccac           >  6:141040/5‑148 (MQ=255)
                                                                                                                                         gctCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTGGGTGGCCAGCTGCGCGACAACCCGccac           >  5:200937/1‑147 (MQ=255)
                                                                                                                                         ccgcGCCCTTGATAAAAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTGGGTGGCCAGCTGCGCGACAACCCGccac           >  8:218477/4‑147 (MQ=255)
                                                                                                                                         ccgcGCCCTTGATAAAAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTGGGTGGCCAGCTGCGCGACAACCCGccac           >  5:131209/4‑147 (MQ=255)
                                                                                                                                          ctCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTa                                                                                <  4:172933/77‑1 (MQ=255)
                                                                                                                                           gcGCCCTTGATAAAAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTGGGTGGCCAGCTGCGCGACAACCCGCCACCAgg       >  2:37787/2‑149 (MQ=255)
                                                                                                                                           gcGCCCTTGATAAAAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTGGGTGGCCAGCTGCGCGACAACCCGCCACCAgg       <  1:37787/148‑1 (MQ=255)
                                                                                                                                            cGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTGGGTGGCCAGCTGCGCGACAACCCGCCACCAGGc      <  6:200937/149‑1 (MQ=255)
                                                                                                                                            cGCCCTTGATAAAAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTGGGTGGCCAGCTGCGCGACAACCCGCCACCAGGc      >  3:55993/1‑149 (MQ=255)
                                                                                                                                             gCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTGGGTGGCCAGCTGCGCGACAACCCGCCACCAGGc      <  3:39087/148‑1 (MQ=255)
                                                                                                                                                cTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTGGGTGGCCAGCTGCGCGACAACCCGCCACCAGGCGTTa  <  3:112956/149‑1 (MQ=255)
                                                                                                                                                 tttttATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTggg                                     >  5:314296/1‑113 (MQ=255)
                                                                                                                                                 tttttATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTggg                                     <  6:314296/113‑1 (MQ=255)
                                                                                                                                                 tttttATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTGGGTGGCCAGCTgcg                         >  3:216342/1‑125 (MQ=255)
                                                                                                                                                 tttttATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTGGGTGGCCAGCTgcg                         <  4:216342/125‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GTCTCGATGAACCTGGTGCACTGGTCGGTGCAGGACGAAATGTGCGGCATGGGTACGCACCCGGAAGAAGACCTGAAGGTGCCGATTTCCCAGCGTTCGGCCTGGTAGCTGCTCGCAGCGCCCCAGTGTGCTAGGCTGCTCGCCCTTTTTATCAAGGGCGCGGCTGCATTGCGCCCTGCAGTGGAGCACGACAATGTCCGAACTCAACCTGTCTACCGACGAAACGCGCGTCAGCTACGGTATCGGCCGTCAGCTGGGTGGCCAGCTGCGCGACAACCCGCCACCAGGCGTTA  >  NC_002947/797835‑798127

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