Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,991,235 G→A 7.7% intergenic (‑104/+125) PP_4399 ← / ← bkdR glutamine synthetase/Bkd operon transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,991,2350GA 7.7% 186.3 / 4.1 65intergenic (‑104/+125)PP_4399/bkdRglutamine synthetase/Bkd operon transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base A (3/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

CATCTGCAGACGACGGCCTTCGAGGAAGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGC  >  NC_002947/4991090‑4991379
                                                                                                                                                 |                                                                                                                                                 
cATCTGCAGACGACGGCCTTCGAGGAAGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCt                                                                                                                                                <  5:122566/149‑1 (MQ=255)
 aTCTGCAGACGACGGCCTTCGAGGAAGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCt                                                                                                                                                >  3:409/1‑148 (MQ=255)
    tGCAGACGACGGCCTTCGAGGAAGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGac                                                                                                                                             >  8:83430/1‑148 (MQ=255)
    tGCAGACGACGGCCTTCGAGGAAGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGac                                                                                                                                             >  4:4739/1‑148 (MQ=255)
      cAGACGACGGCCTTCGAGGAAGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGac                                                                                                                                             <  3:231293/146‑1 (MQ=255)
         acgacgGCCTTCGAGGAAGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGacaacaa                                                                                                                                        >  3:78610/1‑148 (MQ=255)
           gacgGCCTTCGAGGAAGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTCAGCCTacaacaacaa                                                                                                                                     <  3:129471/149‑2 (MQ=255)
           gacgGCCTTCGAGGAAGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGGCAGCCTacaacaacaa                                                                                                                                     <  4:258869/149‑2 (MQ=255)
            acgGCCTTCGAGGAAGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATa                                                                                                                                    <  7:83430/149‑1 (MQ=255)
                    cGAGGAAGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGc                                                                                                                              <  1:324407/147‑1 (MQ=255)
                     gAGGAAGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCTATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCgg                                                                                                                            <  4:235556/148‑1 (MQ=255)
                     gAGGAAGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGTCTGctg                                                                                                                            <  8:69111/148‑3 (MQ=255)
                     gAGGAAGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCgg                                                                                                                            <  6:293040/148‑1 (MQ=255)
                        gAAGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGgtgt                                                                                                                        >  1:220269/1‑149 (MQ=255)
                        gAAGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGgtgt                                                                                                                        >  5:147522/1‑149 (MQ=255)
                          aGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTc                                                                                                                      <  2:317923/149‑1 (MQ=255)
                             cTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCaa                                                                                                                    >  5:294635/1‑148 (MQ=255)
                             cTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGATCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCaa                                                                                                                    >  5:299728/1‑148 (MQ=255)
                                 ggcggTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAata                                                                                                                >  8:233566/1‑148 (MQ=255)
                                      tcatcaCCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAag                                                                                                           >  1:269139/1‑148 (MQ=255)
                                        atcaCCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAg                                                                                                                                   >  2:45180/1‑122 (MQ=255)
                                        atcaCCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAg                                                                                                                                   <  1:45180/122‑1 (MQ=255)
                                                  cgcgcgGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGt                                                                                             >  4:283894/1‑149 (MQ=255)
                                                      cgGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCt                                                                                         <  5:44371/149‑1 (MQ=255)
                                                              cGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGc                                                                                 >  4:22805/1‑149 (MQ=255)
                                                                    ggTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTc                                                                           <  6:147522/149‑1 (MQ=255)
                                                                      tCGGTGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTTCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTcgc                                                                         <  6:299728/149‑1 (MQ=255)
                                                                      tCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTcgc                                                                         <  6:294635/149‑1 (MQ=255)
                                                                          gggTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGt                                                                     <  3:4739/149‑1 (MQ=255)
                                                                          gggTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGt                                                                     <  3:283894/149‑1 (MQ=255)
                                                                          gggTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGt                                                                     <  3:245709/149‑1 (MQ=255)
                                                                              gccGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAAc                                                                 >  2:94955/3‑149 (MQ=255)
                                                                                 gCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCcg                                                              >  2:254097/1‑149 (MQ=255)
                                                                                       gATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGAGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGcc                                                          <  1:254097/148‑1 (MQ=255)
                                                                                           tcgtcgATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAgg                                                    >  4:124936/1‑149 (MQ=255)
                                                                                            cgtcgATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTc                                                                           >  6:145782/1‑125 (MQ=255)
                                                                                            cgtcgATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTc                                                                           >  6:151709/1‑125 (MQ=255)
                                                                                            cgtcgATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTc                                                                           <  5:151709/125‑1 (MQ=255)
                                                                                            cgtcgATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTc                                                                           <  5:145782/125‑1 (MQ=255)
                                                                                            cgtcgATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAgg                                                    <  4:78610/148‑1 (MQ=255)
                                                                                             gtcgATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTAc                                                       >  1:31185/1‑144 (MQ=255)
                                                                                             gtcgATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTAc                                                       <  2:31185/144‑1 (MQ=255)
                                                                                             gtcgATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTg                                                  >  6:167022/1‑149 (MQ=255)
                                                                                             gtcgATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTg                                                  >  5:196981/1‑149 (MQ=255)
                                                                                                      cTGGAGCTGAGCAGTTCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATc                                         <  5:200278/149‑1 (MQ=255)
                                                                                                      cTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATc                                         <  5:205323/149‑1 (MQ=255)
                                                                                                               aGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGcc                                 <  2:220269/148‑1 (MQ=255)
                                                                                                                  aGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAAAATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGATGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGgcg                             >  5:348407/1‑149 (MQ=255)
                                                                                                                    tGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGcgcc                           >  5:34890/1‑149 (MQ=255)
                                                                                                                          tGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGt                                                                     <  4:59319/101‑1 (MQ=255)
                                                                                                                          tGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGt                                                                     >  3:59319/1‑101 (MQ=255)
                                                                                                                                ttGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGcccc                      >  3:119554/1‑142 (MQ=255)
                                                                                                                                ttGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGcccc                      <  4:119554/142‑1 (MQ=255)
                                                                                                                                  gTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAAACCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTccc              <  2:269139/148‑1 (MQ=255)
                                                                                                                                    ctatgttgttgtCAGCCTACAACAACATATGCTGCGGGTGTTAAAAATATCAAGCAGCATATTGGGTGGATGTGCAGGCCTCTTCGAGGTTGAACCAGCTCCTACAGGTGGATGCGATCCTTGTGACGGCGCAGACCCCTAATCACGc            >  1:207904/4‑148 (MQ=255)
                                                                                                                                    ctatgttgttgtCAGCCTACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGc            >  1:258709/4‑148 (MQ=39)
                                                                                                                                    ctatgttgttgtCAGCCTACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGc            >  1:241903/4‑148 (MQ=39)
                                                                                                                                     cTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATc                                         >  6:105080/1‑118 (MQ=255)
                                                                                                                                     cTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATc                                         <  5:105080/118‑1 (MQ=255)
                                                                                                                                     cTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGc            >  4:286911/1‑147 (MQ=255)
                                                                                                                                     cTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGc            >  8:272867/1‑147 (MQ=255)
                                                                                                                                            ttgtAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCAt    <  7:233566/148‑1 (MQ=255)
                                                                                                                                            ttgtAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATg   >  5:258426/1‑149 (MQ=255)
                                                                                                                                              gtATGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGCCATGc  <  8:296730/148‑1 (MQ=255)
                                                                                                                                              gtAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGc  >  3:153939/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CATCTGCAGACGACGGCCTTCGAGGAAGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGC  >  NC_002947/4991090‑4991379

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: