Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,991,239 G→T 7.4% intergenic (‑108/+121) 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,2390GT 7.4% 201.4 / 4.0 63intergenic (‑108/+121)PP_4399/bkdRglutamine synthetase/Bkd operon transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (29/29);  new base T (3/2);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCAGACGACGGCCTTCGAGGAAGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCC  >  NC_002947/4991094‑4991387
                                                                                                                                                 |                                                                                                                                                     
tGCAGACGACGGCCTTCGAGGAAGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGac                                                                                                                                                     >  4:4739/1‑148 (MQ=255)
tGCAGACGACGGCCTTCGAGGAAGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGac                                                                                                                                                     >  8:83430/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                                                                                                                                >  5:147522/1‑149 (MQ=255)
                    gAAGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGgtgt                                                                                                                                >  1:220269/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                                                                                   <  5:145782/125‑1 (MQ=255)
                                                                                        cgtcgATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTc                                                                                   >  6:151709/1‑125 (MQ=255)
                                                                                        cgtcgATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTc                                                                                   <  5:151709/125‑1 (MQ=255)
                                                                                        cgtcgATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTc                                                                                   >  6:145782/1‑125 (MQ=255)
                                                                                        cgtcgATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAgg                                                            <  4:78610/148‑1 (MQ=255)
                                                                                         gtcgATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTAc                                                               <  2:31185/144‑1 (MQ=255)
                                                                                         gtcgATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTAc                                                               >  1:31185/1‑144 (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                              <  4:119554/142‑1 (MQ=255)
                                                                                                                            ttGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGcccc                              >  3:119554/1‑142 (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                    >  8:272867/1‑147 (MQ=255)
                                                                                                                                 cTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGc                    >  4:286911/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)
                                                                                                                                                 gACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGcc  >  6:16569/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TGCAGACGACGGCCTTCGAGGAAGCCTTCGGCGGTCATCACCTTGCCGCGCGGAACGCCGTTGAGGTCGGGGGTGACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCC  >  NC_002947/4991094‑4991387

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: