Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,995,958 A→C 6.6% intergenic (‑204/+81) PP_mr26 ← / ← rfbC C4 ncRNA/dTDP‑4‑dehydrorhamnose 3,5‑epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,995,9580AC 6.6% 234.4 / 7.3 76intergenic (‑204/+81)PP_mr26/rfbCC4 ncRNA/dTDP‑4‑dehydrorhamnose 3,5‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (36/35);  new base C (2/3);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

AGCGACGGGTCGAGACTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACG  >  NC_002947/1995811‑1996105
                                                                                                                                                   |                                                                                                                                                   
aGCGACGGGTCGAGACTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGaaaaa                                                                                                                                                    >  4:282022/1‑149 (MQ=255)
 gCGACGGGTCGAGACTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAAc                                                                                                                                                   <  4:269690/149‑1 (MQ=255)
     cGGGTCGAGACTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCa                                                                                                                                               <  5:157839/149‑1 (MQ=255)
     cGGGTCGAGACTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCa                                                                                                                                               <  6:259791/149‑1 (MQ=255)
     cGGGTCGAGACTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCa                                                                                                                                               <  5:62515/149‑1 (MQ=255)
         tCGAGACTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTTTTTCACACGCCACTaa                                                                                                                                           <  3:285411/149‑2 (MQ=38)
         tCGAGACTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTTTTTCACACGCCACTaa                                                                                                                                           <  5:23087/149‑2 (MQ=38)
         tCGAGACTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTAt                                                                                                                                           >  5:200787/1‑149 (MQ=255)
                tAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCGTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGc                                                                                                                                       >  8:226758/1‑146 (MQ=255)
                  aGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGAc                                                                                                                                  >  5:287649/1‑149 (MQ=255)
                   gTCTCACGTTTAACGTGTGCAATGGAAAACAGTCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGAcc                                                                                                                                 <  6:287649/149‑1 (MQ=255)
                        aCGTTTAACGTGTGCAATGTAAAACCGGCGCAGGGGAATCTTTGGTAAATGGCCTACAAGGATGGGGCAAATTCTTATTATTTGCCAGGCCATGCAGCGGTTACTTCTTAGTGGCGTGTGAAAAACGCTACTATCGGCAAGACCCGGt                                                                                                                             <  1:289164/148‑1 (MQ=255)
                        aCGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCg                                                                                                                               <  3:358803/146‑1 (MQ=255)
                        aCGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTa                                                                                                                            >  6:208270/1‑149 (MQ=255)
                        aCGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTa                                                                                                                            >  7:28497/1‑149 (MQ=255)
                            ttAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATgg                                                                                                                         >  5:144805/1‑148 (MQ=255)
                              aaCGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCAt                                                                                                                      >  6:231897/1‑149 (MQ=255)
                                 gtgtGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCAtgtg                                                                                                                   >  7:108927/1‑149 (MQ=255)
                                 gtgtGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCAtgtg                                                                                                                   >  4:330882/1‑149 (MQ=255)
                                 gtgtGCAATGGAAAACAGGCGAAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCAtgtg                                                                                                                   >  7:297204/1‑149 (MQ=255)
                                    tgCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTc                                                                                                                <  5:231897/149‑1 (MQ=255)
                                       aaTGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCaga                                                                                                             <  6:198677/149‑1 (MQ=255)
                                         tGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAgaga                                                                                                           <  2:62421/149‑1 (MQ=255)
                                           gAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTTTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGAc                                                                                                          <  8:28497/148‑1 (MQ=255)
                                              aaCAGGCGCAGGGGTATCTTTGGGAAATGGCATACAATGAAGTGGCAAATTCTTATTATTTGCAAGGCCATGCAGATGTTTCTTCTTAGTGGCGTGTGAAAAAAGCAACTATCGTAAAGAAACGGTATCGCATGTGTTAATAGAAAAg                                                                                                       >  1:152853/1‑148 (MQ=255)
                                                 aGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTg                                                                                                   >  3:160563/1‑149 (MQ=255)
                                                  ggCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTgc                                                                                                  <  6:200787/149‑1 (MQ=255)
                                                   gcgcAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTc                                                                                                                <  8:125352/134‑1 (MQ=255)
                                                   gcgcAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTc                                                                                                                >  7:125352/1‑134 (MQ=255)
                                                     gcAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTgcgc                                                                                                >  3:185619/1‑148 (MQ=255)
                                                       aGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTAc                                                                                             <  3:330882/149‑1 (MQ=255)
                                                           gTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCg                                                                                           >  7:317819/1‑147 (MQ=255)
                                                           gTATATTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACGCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCa                                                                                         >  7:266310/1‑149 (MQ=255)
                                                            tATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAg                                                                                        >  7:11740/1‑149 (MQ=255)
                                                              tCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGcc                                                                                      >  1:302360/1‑149 (MQ=255)
                                                               cTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGccc                                                                                     >  4:370735/1‑149 (MQ=255)
                                                               cTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGccc                                                                                     >  4:157847/1‑149 (MQ=255)
                                                                tttGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCt                                                                                    >  4:113416/1‑149 (MQ=255)
                                                                  tGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCtt                                                                                   <  8:108927/148‑1 (MQ=255)
                                                                    ggAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGGTAGCGCAGccttgtct                                                                                <  6:320619/149‑7 (MQ=255)
                                                                    ggAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTtctc                                                                                >  8:152568/1‑149 (MQ=255)
                                                                    ggAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTtctc                                                                                >  3:250532/1‑149 (MQ=255)
                                                                    ggAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTtct                                                                                 <  6:102629/148‑1 (MQ=255)
                                                                            ccTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTaaaa                                                                        >  5:367617/1‑149 (MQ=255)
                                                                              tACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAAcc                                                                      >  3:287928/1‑149 (MQ=255)
                                                                                cAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGt                                                                    <  3:100334/149‑1 (MQ=255)
                                                                                  aggaaggGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTa                                                                  <  2:1840/149‑1 (MQ=255)
                                                                                    gaaggGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTATCGCAGCCCTTCTCCGCTAAAACCGTTAGt                                                                <  5:208270/149‑1 (MQ=255)
                                                                                      aggGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCg                                                              <  3:370735/149‑1 (MQ=255)
                                                                                        gggCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGaa                                                            <  3:282022/149‑1 (MQ=255)
                                                                                         tgCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGTCAAGACCCTGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGaa                                                            <  5:59865/147‑1 (MQ=255)
                                                                                         ggCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGaa                                                            >  1:295640/1‑148 (MQ=255)
                                                                                          gCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGc                                                                                          <  3:391008/117‑1 (MQ=255)
                                                                                          gCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGc                                                                                          >  4:391008/1‑117 (MQ=255)
                                                                                          gCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGaa                                                            >  8:290526/1‑147 (MQ=255)
                                                                                          gCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACa                                                          <  7:152568/149‑1 (MQ=255)
                                                                                              aTTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTc                                                      >  2:216221/1‑149 (MQ=255)
                                                                                                   gATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGt                                                 >  3:353926/1‑149 (MQ=255)
                                                                                                       attTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGc                                             >  5:382700/1‑149 (MQ=255)
                                                                                                             cAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGc                                       >  8:120015/1‑149 (MQ=255)
                                                                                                              aGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCg                                      >  6:311426/1‑149 (MQ=255)
                                                                                                                gCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGAc                                    <  4:160563/149‑1 (MQ=255)
                                                                                                                  cATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGAccc                                  <  7:120015/149‑1 (MQ=255)
                                                                                                                    tGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTgcgc                                                                                                <  2:205771/85‑1 (MQ=255)
                                                                                                                    tGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTgcgc                                                                                                >  1:205771/1‑85 (MQ=255)
                                                                                                                      cAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAg                              <  8:317819/149‑1 (MQ=255)
                                                                                                                      cAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAg                              <  1:216221/149‑1 (MQ=255)
                                                                                                                      cAGAGGTTGCTTCTTAGTGTCGTTTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCatg                              <  8:266310/149‑3 (MQ=255)
                                                                                                                       aGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGc                             <  3:113416/149‑1 (MQ=255)
                                                                                                                               cttcttAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCCAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCtt                     <  8:297204/149‑1 (MQ=255)
                                                                                                                                 tcttAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGAcccc                                 >  2:128801/1‑135 (MQ=255)
                                                                                                                                 tcttAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCACAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGAcccc                                 <  1:128801/135‑1 (MQ=255)
                                                                                                                                 cgatAGTGGCGTTTTTCACACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGc                   >  1:251455/4‑149 (MQ=255)
                                                                                                                                 cgatAGTGGCGTTTTTCACACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGc                   <  5:279993/146‑1 (MQ=255)
                                                                                                                                 cgatAGTGGCGTTTTTCACACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGc                   >  6:279993/4‑149 (MQ=255)
                                                                                                                                  cttAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCtt                     <  3:346124/146‑1 (MQ=255)
                                                                                                                                              tgAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACACCAAGCTGGTCCTTGCCAGACAAGCTAgg      >  7:334693/1‑149 (MQ=255)
                                                                                                                                                aaaaaCGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTa    >  5:51677/1‑149 (MQ=255)
                                                                                                                                                 aaaaCGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACCAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTAc   >  4:248027/1‑149 (MQ=255)
                                                                                                                                                 aaaaCGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTAc   >  4:253695/1‑149 (MQ=255)
                                                                                                                                                  aaaCGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACg  <  4:250532/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGCGACGGGTCGAGACTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGCCAGACAAGCTTGGTACG  >  NC_002947/1995811‑1996105

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: