Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,586,720 G→A 6.0% A229V (GCG→GTG)  PP_1391 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,586,7200GA 6.0% 199.8 / 2.1 66A229V (GCG→GTG) PP_1391LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (29/33);  new base A (2/2);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

ATGTGTTTGCGCGACGTATACACCGCATTCAGGTTCAGTTCACGCGGTTTGGCCTCAGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACAAGCTGTGTCCGAAGTAGGCATGGGTCAGGCAGTTGTGCAGGCTCAGTTCCTCAACCC  >  NC_002947/1586579‑1586867
                                                                                                                                             |                                                                                                                                                    
aTGTGTTTGCGCGACGTATACACCGCATTCAGGTTCAGTTCACGCGGTTTGGCCTCAGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCt                                                                                                                                               >  2:102744/1‑149 (MQ=255)
 tgtgTTTGCGCGACGTATACACCGCATTCAGGTTCAGTTCACGCGGTTTGGCCTCAGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCtt                                                                                                                                              <  7:127874/149‑1 (MQ=255)
     tttGCGCGACGTATACACCGCATTCAGGTTCAGTTCACGCGGTTTGGCCTCAGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCACAGCCCTTTg                                                                                                                                            >  3:145608/1‑147 (MQ=255)
     tttGCGCGACGTATACACCGCATTCAGGTTCAGTTCACGCGGTTTGGCCTCAGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCACAGCCCTTTg                                                                                                                                            >  3:151198/1‑147 (MQ=255)
        gcgcGACGTATACACCGCATTCAGGTTCAGTTCACGCGGTTTGGCCTCAGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAgg                                                                                                                                        <  2:109821/148‑1 (MQ=255)
         cgcgACGTATACACCGCATTCAGGTTCAGTTCACGCGGTTTGGCCTCAGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAggg                                                                                                                                       <  8:12195/148‑1 (MQ=255)
         cgcgACGTATACACCGCATTCAGGTTCAGTTCACGCGGTTTGGCCTCAGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGt                                                                                                                                      <  4:145069/149‑1 (MQ=255)
          gcgACGTATACACCGCATTCAGGTTCAGTTCACGCGGTTTGGCCTCAGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGtc                                                                                                                                     <  8:231150/149‑1 (MQ=255)
          gcgACGTATACACCGCATTCAGGTTCAGTTCACGCGGTTTGGCCTCAGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGt                                                                                                                                      >  4:367532/1‑148 (MQ=255)
                  taCACCGCATTCAGGTTCAGTTCACGCGGTTTGGCCTCAGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCGCCTc                                                                                                                             >  3:142783/1‑149 (MQ=255)
                  taCACCGCATTCAGGTTCAGTTCACGCGGTTTGGCCTCAGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTc                                                                                                                             >  6:71382/1‑149 (MQ=255)
                    cacCGCATTCAGGTTCAGTTCACGCGGTTTGGCCTCAGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGt                                                                                                                           >  8:272320/1‑149 (MQ=255)
                     acCGCATTCAGGTTCAGTTCACGCGGTTTGGCCTCAGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGtt                                                                                                                          <  5:71382/149‑1 (MQ=255)
                     acCGCATTCAGGTTCAGTTCACGCGGTTTGGCCTCAGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGtt                                                                                                                          >  6:146201/1‑149 (MQ=255)
                         cATTCAGGTTCAGTTCACGCGGTTTGGCCTCAGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCa                                                                                                                      <  5:326289/149‑1 (MQ=255)
                            tCAGGTTCAGTTCACGCGGTTTGGCCTCAGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCTGGCTGGCGCAggtatc                                                                                                                                              >  1:201283/1‑117 (MQ=255)
                            tCAGGTTCAGTTCACGCGGTTTGGCCTCAGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCTGGCTGGCGCAggtatc                                                                                                                                              <  2:201283/122‑6 (MQ=255)
                                         aCGCGGTTTGGCCTCAGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATa                                                                                                     >  8:163191/1‑149 (MQ=255)
                                           gcgGTTTGGCCTCAGGCAGCAGACGCTCCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCGTTTGCAGGGTCTTCACCTCGTTGGGTCTGATGT‑TTCCGGATACg                                                                                                   >  6:54032/1‑149 (MQ=255)
                                                   gCCTCAGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACgg                                                                                            <  7:353642/148‑1 (MQ=255)
                                                     cTCAGGCAGCAGAAGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGa                                                                                         >  2:133592/1‑149 (MQ=255)
                                                        aGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGAc                                                                                        <  7:272320/147‑1 (MQ=255)
                                                        aGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGAc                                                                                        <  5:54032/147‑1 (MQ=255)
                                                         ggCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGGTGGCACTGATGT‑TTACGGATACGGGCACGGCGACctc                                                                                     >  3:232756/1‑148 (MQ=255)
                                                              cagACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGc                                                                                <  7:273911/149‑1 (MQ=255)
                                                                gACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCg                                                                              <  5:146201/149‑1 (MQ=255)
                                                                     aCCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTAc                                                                         >  6:29335/1‑149 (MQ=255)
                                                                     aCCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTAc                                                                         >  4:318989/1‑149 (MQ=255)
                                                                     aCCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTAc                                                                         >  5:115783/1‑149 (MQ=255)
                                                                        aGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACcgc                                                                      <  2:300143/149‑1 (MQ=255)
                                                                               cccGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGc                                                               >  4:46547/1‑149 (MQ=255)
                                                                               cccGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGGCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGc                                                               <  8:131519/149‑1 (MQ=255)
                                                                                ccGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGcc                                                              >  4:263857/1‑149 (MQ=255)
                                                                                 cGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCTGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCca                                                             >  3:118710/1‑149 (MQ=255)
                                                                                 cGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCca                                                             >  8:32475/1‑149 (MQ=255)
                                                                                  ggCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCcac                                                            <  5:314720/149‑1 (MQ=255)
                                                                                    cGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACaa                                                          <  7:32475/149‑1 (MQ=255)
                                                                                     gTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCACAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACAAg                                                         <  4:151198/149‑1 (MQ=255)
                                                                                     gTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCACAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACAAg                                                         <  4:145608/149‑1 (MQ=255)
                                                                                       aTGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCAAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGAGACCTGCTGGCCGGCTACCGCGAAATGCCACAAGCt                                                       >  6:251095/1‑149 (MQ=255)
                                                                                        tGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCTGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACAAGCt                                                       <  4:118710/148‑1 (MQ=255)
                                                                                        tGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACAAGCt                                                       <  1:88152/148‑1 (MQ=255)
                                                                                         ggcggcTGCGGCCTGGTAGGTTGGCATGATGGTGATAACTGCGCCATCCATGGCATCCCTTTGCATGGTCATCACCTCGTTGGCACTGATGTCTTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACAAGCt                                                       <  1:102744/148‑1 (MQ=255)
                                                                                            gTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACAAGCTGTGTc                                                  <  3:258806/149‑1 (MQ=255)
                                                                                                 cGGCCTGGTATGTTGGCAGCAGTGTGATACGTGCGCCAGCCAGCGCAGCCCTTTGCAGGGGCATCACCTCGTTGGCACGGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACAAGCTGTGTCCg                                                <  4:232756/146‑1 (MQ=255)
                                                                                                  ggCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCGCTGATGT‑TTCCGGATACGGGCACGGCGACctgctg                                                                                  >  4:173152/1‑111 (MQ=255)
                                                                                                  ggCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCGCTGATGT‑TTCCGGATACGGGCACGGCGACctgctg                                                                                  <  3:173152/111‑1 (MQ=255)
                                                                                                  ggCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCg                                                                              >  5:319408/1‑115 (MQ=255)
                                                                                                  ggCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCg                                                                              <  6:319408/115‑1 (MQ=255)
                                                                                                   gCCTGGTAGGTTGGCAGCATGGTTATACCTGCGCCAGCCAGCGCAGCCATAAGCAGGGTCAGCACCTCGTTGGCACTGATGT‑GTCCGGATACGGGCACGGCGACATGCAGGCCGGCTACCGCGAAAACCCAAAAGATGTGTCCGAAg                                             >  7:57814/1‑147 (MQ=255)
                                                                                                   gCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGa                                                                                         <  8:338514/103‑1 (MQ=255)
                                                                                                   gCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGa                                                                                         >  7:338514/1‑103 (MQ=255)
                                                                                                    ccTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACAAGCTGTGTCCGAAGTAg                                          >  6:138714/1‑149 (MQ=255)
                                                                                                    ccTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACAAGCTGTGTCCGAAGTAg                                          >  4:313559/1‑149 (MQ=255)
                                                                                                        gTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACAAGCTGTGTCCGAAGTAGGCAt                                      <  5:29335/149‑1 (MQ=255)
                                                                                                           ggTTGGCAGCATGGTGATACGTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGTTACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACAAGCTGTGTCCGAAGTAGGCATgg                                    <  2:207694/148‑1 (MQ=255)
                                                                                                               ggCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACAAGCTGTGTCCGAAGTAGGCATGGGTCa                                <  3:46547/148‑1 (MQ=255)
                                                                                                                   gcaTGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACAAGCTGTGTCCGAAGTAGGCATGGGTCAGGCAg                           <  8:57814/149‑1 (MQ=255)
                                                                                                                   gcaTGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACAAGCTGTGTCCGAAGTAGGCATGGGTCAGGCAg                           <  3:367532/149‑1 (MQ=255)
                                                                                                                     aTGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACAAGCTGTGTCCGAAGTAGGCATGGGTCAGGCAGtt                         <  6:115783/149‑1 (MQ=255)
                                                                                                                       ggTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACAAGCTGTGTCCGAAGTAGGCATGGGTCAGGCAGTtgt                       >  8:16986/1‑149 (MQ=255)
                                                                                                                             aCCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGaaa                                                                  <  5:282903/100‑1 (MQ=255)
                                                                                                                             aCCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACAAGCTGTGTCCGAAGTAGGCATGGGTCAGGCAGTTGTGCAGGc                 >  5:220129/1‑149 (MQ=255)
                                                                                                                             aCCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑ATCCGGATACGGGCACGGCGACCTGCTGGCCCGCTACAGCGaaa                                                                  >  6:282903/1‑100 (MQ=255)
                                                                                                                                   gccagccaGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACAAGCTGTGTCCGAAGTAGGCATGGGTCAGGCAGTTGTGCAGGCTCAGtt           >  6:79214/1‑149 (MQ=255)
                                                                                                                                       gccagcGCAGCCCTTTGCAGGGTCATCGCCTCGTCGGCACTGATGT‑TTCCGGATACGTGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACAAGCTGTGTCCGAAGTAGGCATGGGTCAGGCAGTTGTGCAGGCTCAGTTCCTc       <  4:142783/149‑1 (MQ=255)
                                                                                                                                             gcAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACAAGCTGTGTCCGAAGTAGGCATGGGTCAGGCAGTTGTGCAGGCTCAGTTCCTCAAccc  <  7:163191/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
ATGTGTTTGCGCGACGTATACACCGCATTCAGGTTCAGTTCACGCGGTTTGGCCTCAGGCAGCAGACGCACCAGTTCACCCCGGCGTATGGCGTCTGCGGCCTGGTAGGTTGGCAGCATGGTGATACCTGCGCCAGCCAGCGCAGCCCTTTGCAGGGTCATCACCTCGTTGGCACTGATGT‑TTCCGGATACGGGCACGGCGACCTGCTGGCCGGCTACCGCGAAATGCCACAAGCTGTGTCCGAAGTAGGCATGGGTCAGGCAGTTGTGCAGGCTCAGTTCCTCAACCC  >  NC_002947/1586579‑1586867

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