Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,804,786 G→A 6.1% R200H (CGT→CAT)  rbsR → DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,804,7860GA 6.1% 197.8 / 2.1 65R200H (CGT→CAT) rbsRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (37/24);  new base A (2/2);  total (39/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

TCGATGCCGACCTGGTGCGCATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGATTTTTGCCGGCAACGACATGATCGGCTTT  >  NC_002947/2804639‑2804934
                                                                                                                                                   |                                                                                                                                                    
tCGATGCCGACCTGGTGCGCATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGt                                                                                                                                                     <  5:355616/149‑1 (MQ=255)
tCGATGCCGACCTGGTGCGCATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGt                                                                                                                                                     <  2:370912/149‑1 (MQ=255)
      ccGACCTGGTGCGCATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCTCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGt                                                                                                                                                     <  6:270819/143‑1 (MQ=255)
      ccGACCTGGTGCGCATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCACCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGt                                                                                                                                                     >  5:270819/1‑143 (MQ=255)
       cGACCTGGTGCGCATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATgg                                                                                                                                              <  8:247593/149‑1 (MQ=255)
        gACCTGGTGCGCATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGc                                                                                                                                             >  5:271133/1‑149 (MQ=255)
         aCCTGGTGCGCATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGcc                                                                                                                                            <  5:111412/149‑1 (MQ=255)
           cTGGTGCGCATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGa                                                                                                                                          <  5:314081/149‑1 (MQ=255)
              gTGCGCATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGaa                                                                                                                                         <  8:306005/147‑1 (MQ=255)
                   cATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTg                                                                                                                                  <  3:282796/149‑1 (MQ=255)
                   cATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTg                                                                                                                                  >  3:315582/1‑149 (MQ=255)
                     ttGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCGGCCATGGCACGGCGGAAGCCgctcag                                                                                                                                <  6:151301/149‑6 (MQ=255)
                       gACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGTTgcgc                                                                                                                               <  4:280268/148‑1 (MQ=255)
                          cACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGt                                                                                                                           >  3:212858/1‑149 (MQ=255)
                                    gggCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGc                                                                                                                 <  5:11664/149‑1 (MQ=255)
                                       cGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCgt                                                                                                              >  7:271522/1‑149 (MQ=255)
                                        gTACCTGGCCACCCGGCATCTGCTGGAGCTTGTCCATCGCGACGTCGCCTATAGCTGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCgtg                                                                                                             <  8:278698/149‑1 (MQ=255)
                                                      ggCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGAc                                                                                               >  2:164725/1‑149 (MQ=255)
                                                          tCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCa                                                                                           >  1:369562/1‑149 (MQ=255)
                                                          tCTGCTGGAGCTTGGCCATCGAGACGTCGCCTATATCGGCGGCCCAGCAGAAACAGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTc                                                                                            >  1:363894/1‑148 (MQ=255)
                                                                 gAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGcc                                                                                     <  1:164725/148‑1 (MQ=255)
                                                                  aGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCg                                                                                   >  7:249188/1‑149 (MQ=255)
                                                                  aGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCg                                                                                   >  2:363235/1‑149 (MQ=255)
                                                                     ttGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGt                                                                                >  6:246605/1‑149 (MQ=255)
                                                                       ggCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTg                                                                               <  3:176094/148‑1 (MQ=255)
                                                                        gCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGc                                                                             >  2:115240/1‑149 (MQ=255)
                                                                         ccATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGcc                                                                            <  4:315582/149‑1 (MQ=255)
                                                                                   gTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAgcggc                                                                  <  8:271522/149‑1 (MQ=255)
                                                                                      gCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAgcggcggc                                                               <  1:115240/149‑1 (MQ=255)
                                                                                       ccggtttCGGCGGGGCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAgcggcgg                                                                <  7:253252/141‑1 (MQ=255)
                                                                                       ccggtttCGGCGGGGCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAgcggcgg                                                                >  8:253252/7‑147 (MQ=255)
                                                                                       ccTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAgcggcggcg                                                              >  1:151277/1‑149 (MQ=255)
                                                                                        cTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAgcggcggcg                                                              >  2:111631/1‑148 (MQ=255)
                                                                                          atatCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAg                                                           >  5:198737/1‑149 (MQ=255)
                                                                                             tCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTg                                                        >  6:85840/1‑149 (MQ=255)
                                                                                                gcggcCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTg                                                     >  8:2156/1‑149 (MQ=255)
                                                                                                    ccccGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACGTCACCAGCCCGGGTGGCCAAGCAGCGGCGGCGCCGGTGCGGGAAg                                                 >  7:312713/1‑149 (MQ=255)
                                                                                                     cccGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAgg                                                >  3:99641/1‑149 (MQ=255)
                                                                                                     cccGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCACCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAgg                                                >  1:71073/1‑149 (MQ=255)
                                                                                                      ccgccgAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGaa                                                  >  1:187717/1‑146 (MQ=255)
                                                                                                      ccgccgAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGaa                                                  <  2:187717/146‑1 (MQ=255)
                                                                                                      ccgccgAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAgg                                                >  7:248226/1‑148 (MQ=255)
                                                                                                       cgccgAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGc                                               >  6:345586/1‑148 (MQ=255)
                                                                                                        gccgAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGc                                               >  2:56085/1‑147 (MQ=255)
                                                                                                         ccgAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGc                                                       <  4:187421/138‑1 (MQ=255)
                                                                                                         ccgAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGc                                                       >  3:187421/1‑138 (MQ=255)
                                                                                                                cGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGACGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCt                                     >  1:76016/1‑149 (MQ=255)
                                                                                                                  gggTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCg                                   >  4:229018/1‑149 (MQ=255)
                                                                                                                  gggTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCg                                   <  5:246605/149‑1 (MQ=255)
                                                                                                                   ggTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCgg                                  >  8:319091/1‑149 (MQ=255)
                                                                                                                         ccAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCACGCGGCCCTCGGCGAttt                            >  4:33581/1‑149 (MQ=255)
                                                                                                                         ccAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGc                                           >  5:88803/1‑134 (MQ=255)
                                                                                                                         ccAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGc                                           <  6:88803/134‑1 (MQ=255)
                                                                                                                         ccAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCgg                                  <  6:345732/143‑1 (MQ=255)
                                                                                                                         ccAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCgg                                  >  5:345732/1‑143 (MQ=255)
                                                                                                                         ccAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGAttt                            >  5:90398/1‑149 (MQ=255)
                                                                                                                                     caccGGCTTCGGCCATGGCACGGCGGAAGCCGCTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGaa                                                  >  6:126644/4‑115 (MQ=37)
                                                                                                                                     caccGGCTTCGGCCATGGCACGGCGGAAGCCGCTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGaa                                                  <  5:126644/112‑1 (MQ=37)
                                                                                                                                      accGGCTTCGGCCATGGCACGGCGGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGATTTTTGCCGGCAACGa               >  3:166814/3‑149 (MQ=255)
                                                                                                                                           cTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAg                                                 >  3:277587/1‑110 (MQ=255)
                                                                                                                                           cTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAg                                                 <  4:277587/110‑1 (MQ=255)
                                                                                                                                           cTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGATTTTTGCCGGCAACGACATGa          >  5:194455/1‑149 (MQ=255)
                                                                                                                                            ttCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGATTTTTGCCgg                     <  3:32929/137‑1 (MQ=255)
                                                                                                                                            ttCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGATTTTTGCCgg                     >  4:32929/1‑137 (MQ=255)
                                                                                                                                            ttCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGATTTTTGCCGGCAACGACATGAt         >  5:314577/1‑149 (MQ=255)
                                                                                                                                                  cgTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGATTTTTGCCGGCAACGACATGATCGGCtt   <  6:314577/149‑1 (MQ=255)
                                                                                                                                                  cgTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGATTTTTGCCGGCAACGACATGATCGGCtt   <  4:212858/149‑1 (MQ=255)
                                                                                                                                                   gTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGATTTTTGCCGGCAACGACATGATCGGCttt  <  5:85840/149‑1 (MQ=255)
                                                                                                                                                   gTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGATTTTTGCCGGCAACGACATGATCGGCttt  >  1:17847/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TCGATGCCGACCTGGTGCGCATTGACCACGAGCAGGGGGCGTACCTGGCCACCCGGCATCTGCTGGAGCTTGGCCATCGCGACGTCGCCTATATCGGCGGCCCCGCCGAAACCGGGGTTACCCAGCTGCGCCTGAGCGGCTTCCGCCGTGCCATGGCCGAAGCCGGTGCGCCCGTGCCGGGCAGCCGTGTACTGCACTGCGACTTCACCAGCCCGGGTGGCCATGCAGCGGCGGCGCAGGTGCTGGAAGGCAAGCGGCCCTCGGCGATTTTTGCCGGCAACGACATGATCGGCTTT  >  NC_002947/2804639‑2804934

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: ‑