Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,991,313 G→T 8.0% intergenic (‑182/+47) 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,3130GT 8.0% 194.2 / 5.4 62intergenic (‑182/+47)PP_4399/bkdRglutamine synthetase/Bkd operon transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (33/24);  new base T (2/3);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

ACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTAGCGCACCTGCTTCAGGGCAAAACTCGAGCGAATATTGGCTACCCCTGG  >  NC_002947/4991169‑4991458
                                                                                                                                                |                                                                                                                                                 
aCGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAAc                                                                                                                                                <  2:281001/148‑1 (MQ=255)
  gCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCcg                                                                                                                                             >  1:58361/1‑149 (MQ=255)
   cATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCcgc                                                                                                                                            >  4:131729/1‑149 (MQ=255)
   cATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCcgc                                                                                                                                            <  6:310244/149‑1 (MQ=255)
   cATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCcgc                                                                                                                                            <  8:224307/149‑1 (MQ=255)
    aTTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTTCACCcgc                                                                                                                                            >  2:117536/1‑148 (MQ=255)
    aTTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCcgcg                                                                                                                                           >  6:149020/1‑149 (MQ=255)
     ttCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTTCACCcgcga                                                                                                                                          <  8:22881/149‑2 (MQ=255)
     ttCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCgcgc                                                                                                                                          >  4:255277/1‑149 (MQ=255)
        gATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTTCACCcgcgaaga                                                                                                                                       <  1:117536/149‑5 (MQ=255)
         aTTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTAc                                                                                                                                      <  3:131729/149‑1 (MQ=255)
          tttCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTTCACCcgcgaa                                                                                                                                         >  3:132823/1‑143 (MQ=255)
          tttCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTTCACCcgcgaa                                                                                                                                         <  4:132823/145‑3 (MQ=255)
                          gAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCtt                                                                                                                     <  3:255277/149‑1 (MQ=255)
                            gCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTtgt                                                                                                                   <  3:166521/149‑1 (MQ=255)
                             cTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTtgtg                                                                                                                  >  4:347310/1‑149 (MQ=255)
                                               gCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACt                                                                                                >  4:324421/1‑149 (MQ=255)
                                                     ctatgttgttgtCAGCCTACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGc                                                                                           >  5:138875/4‑148 (MQ=39)
                                                      tatgttgttgtCAGCCTACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGc                                                                                           >  7:40194/3‑147 (MQ=255)
                                                      cTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGc                                                                                           >  3:274147/1‑147 (MQ=255)
                                                           tgttgtAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCa                                                                                    >  6:226155/1‑149 (MQ=255)
                                                                         aacaacaTAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGgcgc                                                                                                           <  3:357516/112‑1 (MQ=255)
                                                                         aacaacaTAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGgcgc                                                                                                           >  4:357516/1‑112 (MQ=255)
                                                                         aacaacaTAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCcgg                                                                       >  1:8815/1‑148 (MQ=255)
                                                                         aacaacaTAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCCCAAGGTCATGCCATTGGCcgg                                                                       >  3:260688/1‑148 (MQ=255)
                                                                           caacaTAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCAt                                                                                   <  1:283825/134‑1 (MQ=255)
                                                                           caacaTAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCAt                                                                                   >  2:283825/1‑134 (MQ=255)
                                                                            aacaTAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCcggcaac                                                                   >  2:162790/1‑149 (MQ=255)
                                                                              caTAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGt                                                                                                                                  <  3:325185/84‑1 (MQ=255)
                                                                              caTAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGt                                                                                                                                  >  4:325185/1‑84 (MQ=255)
                                                                               aTAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCcggcaac                                                                   >  7:293788/1‑146 (MQ=255)
                                                                               aTAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCcggcaac                                                                   <  8:293788/146‑1 (MQ=255)
                                                                                   gCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGgcaacggcaacg                                                            >  3:103402/1‑149 (MQ=255)
                                                                                        gggTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTc                                                       >  5:361894/1‑149 (MQ=255)
                                                                                           tgtTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTg                                                    <  4:274147/149‑1 (MQ=255)
                                                                                             ttCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTa                                                  <  5:149020/149‑1 (MQ=255)
                                                                                             ttCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTa                                                  <  5:226155/149‑1 (MQ=255)
                                                                                               cAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTAgc                                                >  5:270330/1‑149 (MQ=255)
                                                                                                aaaaTATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCcggcaacg                                                                  <  4:6106/130‑1 (MQ=255)
                                                                                                aaaaTATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCcggcaacg                                                                  >  3:6106/1‑130 (MQ=255)
                                                                                                   atatCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCTCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTAGCGCAc                                            <  1:162790/149‑1 (MQ=255)
                                                                                                        aaGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTAGCGCACCTGc                                        <  6:138875/148‑1 (MQ=255)
                                                                                                         agcagcCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTAGCGCACCTGCtt                                      >  6:244519/1‑149 (MQ=255)
                                                                                                             gcCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTAGCGCACCTGCTTCAgg                                  <  2:8815/149‑1 (MQ=255)
                                                                                                              ccTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTAGCGCACCTGCTTCAggg                                 >  4:142481/1‑149 (MQ=255)
                                                                                                                 aTTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCACCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTAGCGCACCTGCTTCAGGGCa                               >  4:16008/1‑148 (MQ=255)
                                                                                                                   tGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCAGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATTCCATTGGCCGGCAACGGCAACGCCGTCTTGTAGCGCACCTGCTTCAGGGCaa                              >  3:20443/1‑147 (MQ=255)
                                                                                                                    gggTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTAGCGCACCTGCTTCAGGGCaaaa                            <  4:103402/148‑1 (MQ=255)
                                                                                                                       tGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCtt                                                     <  4:28914/120‑1 (MQ=255)
                                                                                                                       tGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCtt                                                     >  3:28914/1‑120 (MQ=255)
                                                                                                                            gtgCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTAGCGCACCTGCTTCAGGGCAAAACTCGAGCg                    >  3:213171/1‑148 (MQ=255)
                                                                                                                               ccGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTAGCGCACCTGCTTCAGGGCAAAACTCGa                       <  5:319113/142‑1 (MQ=255)
                                                                                                                               ccGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTAGCGCACCTGCTTCAGGGCAAAACTCGa                       >  6:319113/1‑142 (MQ=255)
                                                                                                                                cGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTAGCGCACCTGCTTCAGGGCAAAACTCGAGCGAatat               >  1:154954/1‑149 (MQ=255)
                                                                                                                                 ggCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTAGCGCACCTGCTTCAGGGCAAAACTCGAGCGAATAtt              <  5:244519/149‑1 (MQ=255)
                                                                                                                                    ctctTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTAGCGCACCTGCTTCAGGGCAAAACTCGAGCGAATATTGGc           <  8:40194/149‑1 (MQ=255)
                                                                                                                                         cgcgGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCACGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTAGCGCACCTGCTTCAGGGCAAAACTCGAGCGAATATTGGCTAc        <  4:213171/147‑1 (MQ=255)
                                                                                                                                         cgcgGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGCGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTAGCGCACCTGCTTCAGGGCAAAACTCGAGCGAATATTGGCTAccc      >  4:275865/1‑149 (MQ=255)
                                                                                                                                         cgcgGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGCGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTAGCGCACCTGCTTCAGGGCAAAACTCGAGCGAATATTGGCTAccc      >  4:270855/1‑149 (MQ=255)
                                                                                                                                             ggTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTAGCGCACCTGCTTCAGGGCAAAACTCGAGCGAATATTGGCTACCCCTgg  <  6:84156/149‑1 (MQ=255)
                                                                                                                                             ggTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTAGCGCACCTGCTTCAGGGCAAAACTCGAGCGAATATTGGCTACCCCTgg  >  4:190383/1‑149 (MQ=255)
                                                                                                                                              gTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTAGCGCACCTGCTTCAGGGCAAAACTCGAGCGAATATTGGCTACCCCTgg  >  6:57175/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
ACGCATTCGATTTCGTCGATGCCCTGGAGCTGAGCAGTGCTCATGACGCTTGTCCTTGTTGTTGTAGGCTGACAACAACATAGGCTGCGGGTGTTCAAAATATCAAGCAGCCTATTGGGTGGCTGTGCCGGCCTCTTCGCGGGTGAACCCGCGCCTACAGGTGGATGCGATCCTTGTGCCGGCGCCGCCCCCTACTCCCGCAAGGTCATGCCATTGGCCGGCAACGGCAACGCCGTCTTGTAGCGCACCTGCTTCAGGGCAAAACTCGAGCGAATATTGGCTACCCCTGG  >  NC_002947/4991169‑4991458

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: