Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,995,952 G→T 8.2% intergenic (‑198/+87) 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,9520GT 8.2% 237.3 / 5.8 76intergenic (‑198/+87)PP_mr26/rfbCC4 ncRNA/dTDP‑4‑dehydrorhamnose 3,5‑epimerase
Reads supporting (aligned to +/- strand):  ref base G (36/33);  new base T (2/5);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGACGGGTCGAGACTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGC  >  NC_002947/1995811‑1996088
                                                                                                                                             |                                                                                                                                        
aGCGACGGGTCGAGACTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGaaaaa                                                                                                                                   >  4:282022/1‑149 (MQ=255)
 gCGACGGGTCGAGACTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAAc                                                                                                                                  <  4:269690/149‑1 (MQ=255)
     cGGGTCGAGACTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCa                                                                                                                              <  5:62515/149‑1 (MQ=255)
     cGGGTCGAGACTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCa                                                                                                                              <  6:259791/149‑1 (MQ=255)
     cGGGTCGAGACTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCa                                                                                                                              <  5:157839/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                                                                                                           >  7:28497/1‑149 (MQ=255)
                        aCGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTa                                                                                                           >  6:208270/1‑149 (MQ=255)
                            ttAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATgg                                                                                                        >  5:144805/1‑148 (MQ=255)
                              aaCGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCAt                                                                                                     >  6:231897/1‑149 (MQ=255)
                                 gtgtGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCAtgtg                                                                                                  >  4:330882/1‑149 (MQ=255)
                                 gtgtGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCAtgtg                                                                                                  >  7:108927/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                                                                                               >  7:125352/1‑134 (MQ=255)
                                                   gcgcAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTc                                                                                               <  8:125352/134‑1 (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:157847/1‑149 (MQ=255)
                                                               cTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGccc                                                                    >  4:370735/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                                                               >  3:250532/1‑149 (MQ=255)
                                                                    ggAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTtctc                                                               >  8:152568/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                                                                         >  4:391008/1‑117 (MQ=255)
                                                                                          gCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGc                                                                         <  3:391008/117‑1 (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                                                                               >  1:205771/1‑85 (MQ=255)
                                                                                                                    tGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTgcgc                                                                               <  2:205771/85‑1 (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  >  6:279993/4‑149 (MQ=255)
                                                                                                                                 cgatAGTGGCGTTTTTCACACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGc  <  5:279993/146‑1 (MQ=255)
                                                                                                                                  cttAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCtt    <  3:346124/146‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                        
AGCGACGGGTCGAGACTAAGTCTCACGTTTAACGTGTGCAATGGAAAACAGGCGCAGGGGTATCTTTGGGAAATGGCCTACAAGGAAGGGGCAAATTCTGATTATTTGCCAGGCCATGCAGCGGTTGCTTCTTAGTGGCGTGTGAAAAACGCCACTATCGGCAAGACCCGGTATGGCATGTGTTCAGAGACAAGGCTGCGCTACCGCAGCCCTTCTCCGCTAAAACCGTTAGTCGAACAGTTCGGCGTCAGCCAATGCGACCCCAAGCTGGTCCTTGC  >  NC_002947/1995811‑1996088

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: