Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 537,579 C→G 8.4% intergenic (+109/‑44) PP_mr07 → / → rpoB P26 ncRNA/DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947537,5790CG 8.4% 219.4 / 9.5 71intergenic (+109/‑44)PP_mr07/rpoBP26 ncRNA/DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (31/34);  new base G (2/4);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.24e-01

AGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCTGGATTCG  >  NC_002947/537434‑537724
                                                                                                                                                 |                                                                                                                                                 
aGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCgg                                                                                                                                                <  3:351907/149‑1 (MQ=255)
aGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCgg                                                                                                                                                <  5:292598/149‑1 (MQ=255)
 gCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCggg                                                                                                                                               <  8:229524/149‑1 (MQ=255)
 gCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCgg                                                                                                                                                <  3:53465/148‑1 (MQ=255)
  cGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCgg                                                                                                                                                >  3:110756/1‑147 (MQ=255)
     aaGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAg                                                                                                                                          >  3:328560/1‑149 (MQ=255)
       ggCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTTCTCCCGGCGGGagcagc                                                                                                                                         <  5:232196/149‑1 (MQ=255)
        gCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTTCTCCCGGCGGGagcagc                                                                                                                                         >  2:355661/1‑148 (MQ=255)
        gCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTTCTCCCGGCGGGAGCAGcc                                                                                                                                        <  1:355661/149‑1 (MQ=255)
        gCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTTCTCCCGGCGGGAGCAGcc                                                                                                                                        <  1:125244/149‑1 (MQ=255)
         cTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTTCTCCCGGCGGGAGCAGCCa                                                                                                                                       <  5:203625/149‑1 (MQ=255)
               gCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACaa                                                                                                                                  <  5:7880/148‑1 (MQ=255)
               gCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAg                                                                                                                                 <  5:71222/149‑1 (MQ=255)
               gCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAg                                                                                                                                 <  5:76633/149‑1 (MQ=255)
                cTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGCGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGgaga                                                                                                                                            <  2:115958/137‑1 (MQ=255)
                cTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGCGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGgaga                                                                                                                                            >  1:115958/1‑137 (MQ=255)
                 tggtggCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAg                                                                                                                                 >  2:27277/1‑147 (MQ=255)
                   gtggCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAg                                                                                                                              <  1:332746/148‑1 (MQ=255)
                    tggCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGt                                                                                                                            <  2:229690/149‑1 (MQ=255)
                             ccACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGAc                                                                                                                         >  1:319711/1‑143 (MQ=255)
                             ccACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGAc                                                                                                                         <  2:319711/143‑1 (MQ=255)
                              cACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGaa                                                                                                                                           <  7:286111/124‑1 (MQ=255)
                              cACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGaa                                                                                                                                           >  8:286111/1‑124 (MQ=255)
                              cACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGaa                                                                                                                                           <  2:273984/124‑1 (MQ=255)
                              cACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGaa                                                                                                                                           >  1:273984/1‑124 (MQ=255)
                              cACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTg                                                                                                                  >  1:340649/1‑149 (MQ=255)
                               aCCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTgg                                                                                                                 <  6:121531/149‑1 (MQ=255)
                               aCCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTgg                                                                                                                 <  2:340649/149‑1 (MQ=255)
                               aCCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTgg                                                                                                                 >  3:71186/1‑149 (MQ=255)
                                   gCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGaa                                                                                                             <  6:377846/149‑1 (MQ=255)
                                           ccGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAgg                                                                                                                             >  6:281947/1‑125 (MQ=255)
                                           ccGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAgg                                                                                                                             <  5:281947/125‑1 (MQ=255)
                                                   ggtggCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACt                                                                                              >  5:203607/1‑148 (MQ=255)
                                                        ccGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATaca                                                                                        >  3:119246/1‑149 (MQ=255)
                                                               ggCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTgaga                                                                                  >  2:340629/1‑148 (MQ=255)
                                                               ggCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTgaga                                                                                  >  7:215910/1‑148 (MQ=255)
                                                                gCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCa                                                                                                                                       <  5:20958/94‑1 (MQ=255)
                                                                gCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCa                                                                                                                                       >  6:20958/1‑94 (MQ=255)
                                                                  cGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGaaaa                                                                               >  7:215198/1‑148 (MQ=255)
                                                                   ggTCCCGTTAACGAGCTGCAAGTCTGCGATGTGGGTGGCGCATACCTAGGGGTTTGGACGTTTTTCTGGCTGTTCCCGCCGGGAGATTCCAAACATGCAGGTGACCAAGCTGGGGAACGGTGATGTCTTACTCATACTCTGAGtataaa                                                                             <  4:119246/149‑1 (MQ=255)
                                                                   ggTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGaaaaaa                                                                             <  1:27277/149‑1 (MQ=255)
                                                                       ccGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTa                                                                         >  1:61280/1‑149 (MQ=255)
                                                                         gATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATc                                                                       <  4:110756/149‑1 (MQ=255)
                                                                          aTAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATcc                                                                      <  8:215198/149‑1 (MQ=255)
                                                                                    gCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCaa                                                                                                                      <  6:11215/91‑1 (MQ=255)
                                                                                    gCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCaa                                                                                                                      >  5:11215/1‑91 (MQ=255)
                                                                                     cAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACtt                                                            >  2:74688/1‑148 (MQ=255)
                                                                                     cAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTAGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACtt                                                            <  7:206695/148‑1 (MQ=255)
                                                                                       aGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAg                                                         >  3:32972/1‑149 (MQ=255)
                                                                                        gACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAGGGACTTTAg                                                         >  3:116620/1‑148 (MQ=255)
                                                                                        gACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGc                                                        <  6:203607/149‑1 (MQ=255)
                                                                                          cacCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCaa                                                      <  1:162369/149‑1 (MQ=255)
                                                                                           acCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAg                                                     >  3:69807/1‑149 (MQ=255)
                                                                                              cATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTg                                                  >  7:21871/1‑149 (MQ=255)
                                                                                               aTGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGc                                                 >  3:480/1‑149 (MQ=255)
                                                                                                  tgggtggCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAGGGACTTTAGCAAGTTGCCgg                                              <  4:116620/149‑1 (MQ=255)
                                                                                                          gcAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATg                                      <  4:32972/149‑1 (MQ=255)
                                                                                                                cAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGAtgtg                                >  1:285062/1‑149 (MQ=255)
                                                                                                                 aaGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGc                               <  2:285062/149‑1 (MQ=255)
                                                                                                                   ggggTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCt                             <  8:215910/149‑1 (MQ=255)
                                                                                                                    gggTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGAGGTCATGGATGTGCCt                             >  4:361583/1‑148 (MQ=255)
                                                                                                                     ggTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCAt                                                                                           <  5:38020/85‑1 (MQ=255)
                                                                                                                     ggTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCAt                                                                                           >  6:38020/1‑85 (MQ=255)
                                                                                                                       tttGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTAcc                         <  1:74688/149‑1 (MQ=255)
                                                                                                                           cACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTAcctcct                     >  4:44704/1‑149 (MQ=255)
                                                                                                                           cACGATATTCTCGCTGCTCCCGCCGGGAGACGCCAAACAACCACGTGACCAAGCTGCGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGGCATGGATGTGCCTTAcctcct                     >  6:74622/1‑149 (MQ=255)
                                                                                                                            aCGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTg                    <  1:340629/149‑1 (MQ=255)
                                                                                                                                      ggCTTCTCCCGGCGGGAGCAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGc          >  4:83315/1‑149 (MQ=255)
                                                                                                                                         tgctCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCg                                               >  1:245409/1‑109 (MQ=255)
                                                                                                                                         tgctCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCg                                               <  2:245409/109‑1 (MQ=255)
                                                                                                                                         tgctCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGTCTTACCTCCTGGCCATCCAGCTg        >  1:169543/1‑148 (MQ=255)
                                                                                                                                              ccgccgGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCTGGATTCg  >  3:42975/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
AGCGAAAGGCTGATGGCTGGTGGCTTATGCCACCGGCCTTTTTCCGTTATTGGTGGCCGGTCAGGCCGGTCCCGATAACGAGCTGCAAGACACCCATGTGGGTGGCGCAAACCAAGGGGTTTGCACGATTTTCTGGCTGCTCCCGCCGGGAGAAGCCAAACAAGCAGGTGACCAAGCTGGGGAACGCTGATGGCTTACTCATACACTGAGAAAAAACGTATCCGCAAGGACTTTAGCAAGTTGCCGGACGTCATGGATGTGCCTTACCTCCTGGCCATCCAGCTGGATTCG  >  NC_002947/537434‑537724

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