Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,467,957 T→A 8.1% intergenic (+16/+16) queF → / ← PP_2161 NADPH‑dependent 7‑cyano‑7‑deazaguanine reductase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,467,9570TA 8.1% 171.1 / 3.6 59intergenic (+16/+16)queF/PP_2161NADPH‑dependent 7‑cyano‑7‑deazaguanine reductase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/25);  new base A (3/2);  total (32/27)
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

TTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTGCCAGGTTCACGCCGTCCACCAGGTTGTTGTCTGGCGTGGCTTCCT  >  NC_002947/2467810‑2468103
                                                                                                                                                    |                                                                                                                                                   
ttGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCAATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCACCCCGGATAACAAACGGCTGGTCCGGAAGTAATACAGAGCGCCCCATTg                                                                                                                                                    >  4:280253/1‑149 (MQ=255)
   aagaaTCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCg                                                                                                                                                  <  3:280253/148‑1 (MQ=255)
          tgctgcAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGc                                                                                                                                          >  7:147431/1‑149 (MQ=255)
            ctgcAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCAAAGCGGccccc                                                                                                                                            <  5:334887/145‑1 (MQ=255)
             tgcAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTAt                                                                                                                                       >  8:130272/1‑149 (MQ=255)
                aGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCTCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTc                                                                                                                                    <  1:41387/149‑1 (MQ=255)
                aGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCAAAGCGGccccc                                                                                                                                            <  4:326373/141‑1 (MQ=255)
                aGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCAAAGCGGccccc                                                                                                                                            >  3:326373/1‑141 (MQ=255)
                     gAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATc                                                                                                                               >  4:116532/1‑149 (MQ=255)
                            tGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGc                                                                                                                        >  1:228071/1‑149 (MQ=255)
                                 gTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCtgtg                                                                                                                    >  8:68900/1‑148 (MQ=255)
                                  tCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCtgtg                                                                                                                    >  2:310016/1‑147 (MQ=255)
                                  tCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCtgtg                                                                                                                    >  3:150745/1‑147 (MQ=255)
                                     aTGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGc                                                                                                               >  3:289216/1‑148 (MQ=255)
                                      tGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCTCGTTCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATCCAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCg                                                                                                              >  1:102735/1‑149 (MQ=255)
                                      tGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCg                                                                                                              <  3:116532/148‑1 (MQ=255)
                                        cgcgcTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTAttt                                                                                                                                     >  2:170924/1‑124 (MQ=255)
                                        cgcgcTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTAttt                                                                                                                                     <  1:170924/124‑1 (MQ=255)
                                          cgcTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTg                                                                                                         <  7:130272/149‑1 (MQ=255)
                                                       gcgGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCa                                                                                            <  8:224595/149‑1 (MQ=255)
                                                          gTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAggg                                                                                         <  2:228071/149‑1 (MQ=255)
                                                             ggCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTgg                                                                                      <  8:147431/149‑1 (MQ=255)
                                                             ggCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCACCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGTGTgg                                                                                      >  2:365363/1‑149 (MQ=255)
                                                             ggCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGGCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTgg                                                                                      <  1:130993/149‑1 (MQ=255)
                                                                 ggACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCa                                                                                                                  >  5:358902/1‑118 (MQ=255)
                                                                 ggACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCa                                                                                                                  <  6:358902/118‑1 (MQ=255)
                                                                   aCATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGt                                                                                  >  3:78689/1‑147 (MQ=255)
                                                                       cAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGt                                                                                        >  5:63220/1‑137 (MQ=255)
                                                                        aaCCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCg                                                                           <  8:204135/149‑1 (MQ=255)
                                                                        aaCCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCg                                                                           >  7:155944/1‑149 (MQ=255)
                                                                                   cagcaCCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGc                                                                <  1:310016/149‑1 (MQ=255)
                                                                                       tccGGGCTGATCGGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTc                                                             >  1:189253/2‑148 (MQ=255)
                                                                                        ccGGGCTGATCGGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCa                                                                                            <  1:78301/116‑1 (MQ=255)
                                                                                        ccGGGCTGATCGGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCa                                                                                            >  2:78301/1‑116 (MQ=255)
                                                                                                gATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCATCCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCt                                                     <  1:293497/149‑1 (MQ=255)
                                                                                                     gCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTg                                               >  7:121000/1‑149 (MQ=255)
                                                                                                                aaaCGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTGCCAGGTTCACg                                    >  2:290835/1‑149 (MQ=255)
                                                                                                                 aaCGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTGCCAGGTTCACGc                                   >  2:235840/1‑149 (MQ=255)
                                                                                                                  aCGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCa                                                                                    <  1:180296/99‑1 (MQ=255)
                                                                                                                  aCGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCa                                                                                    >  2:180296/1‑99 (MQ=255)
                                                                                                                     gCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTGCCAGGTTCACGCCGTc                               <  5:258001/149‑1 (MQ=255)
                                                                                                                              gCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTGCCAGGTTCACGCCGTCCACCAGgtt                      >  1:185217/1‑149 (MQ=255)
                                                                                                                                    ctACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTGCCAGGTTCACGCCGTCCACCAGGTTGTTGTc                <  4:289216/148‑1 (MQ=255)
                                                                                                                                    atagcgggggCCGCAAAGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTGCCAGGTTCACGCCGTCCACCAGGTTGTTGTc                >  8:99328/6‑149 (MQ=255)
                                                                                                                                      aCAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTGCCAGGTTCACGCCATCCACCAGGTTGTTGTCTg              >  1:299730/1‑149 (MQ=255)
                                                                                                                                         gggggCCGCAAAGCGGCCCCCGCTATTTCAGATCCCAATGCTGTGCAGCGAGTTGGCAATGCAGCGCAGGGTCGCAGTGAGGTCGGGGTGGTATGCATCGAAGCGCTCGATTGCCAGGTTCACGCCGTCCACCATGTTGTTGTCTggag           >  2:295261/1‑147 (MQ=255)
                                                                                                                                          ggggCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTGCCAGGTTCACGCCGTc                               <  5:268343/128‑1 (MQ=255)
                                                                                                                                          ggggCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTGCCAGGTTCACGCCGTc                               >  6:268343/1‑128 (MQ=255)
                                                                                                                                          ggggCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTGCCAGGTTCACGCCGTCCACCAGgttgtt                   >  2:220261/1‑140 (MQ=255)
                                                                                                                                          ggggCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTGCCAGGTTCACGCCGTCCACCAGgttgtt                   <  1:220261/140‑1 (MQ=255)
                                                                                                                                          ggggCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTGCCAGGTTCACGCCGTCCACCAGGTTGTTGTCTGGCGt          >  4:65224/1‑149 (MQ=255)
                                                                                                                                          ggggCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTGCCAGGTTCACGCCGTCCACCAGGTTGTTGTCTGGCGt          <  4:150745/149‑1 (MQ=255)
                                                                                                                                          ggggCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTGCCAGGTTCACGCCGTCCACCAGGTTGTTGTCTGGCGt          <  3:65224/149‑1 (MQ=255)
                                                                                                                                           gggCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTGCCAGGTTCACGCCGTCCACCAGGTTGTTGTCTGGCGTg         <  8:121000/149‑1 (MQ=255)
                                                                                                                                           gggCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTGCCAGGTTCACGCCGTCCACCAGGTTGTTGTCTGGCGTg         <  1:290835/149‑1 (MQ=255)
                                                                                                                                              ccGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTGCCAGGTTCACGCCGTCCACCAGGTTGTTGTCTGGCGTGGCt      <  4:78689/149‑1 (MQ=255)
                                                                                                                                                gCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTGCCAGGTTCACGCCGTCCCCCAGGTTGTTGTCTGGCGTGGCTTc    >  3:233128/1‑149 (MQ=255)
                                                                                                                                                gCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGGCGGCCTCGAAGCGCTCGATTGCCAGGTTCACGCCGTCCACCAGGTTGTTGTCTGGCGTGGCTTc    <  8:155944/149‑1 (MQ=255)
                                                                                                                                                 cTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTGCCAGGTTCACGCCGTCCACCAGGTTGTTGTCTGGCGTGGCTTcc   >  3:78800/1‑149 (MQ=255)
                                                                                                                                                    tGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTGCCAGGTTCACGCCGTCCACCAGGTTGTTGTCTGGCGTGGCTTCCt  <  4:168891/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TTGAAGAATCTGCTGCAGCCGGAGCACCTGACGGTCTATGCGCGCTATGTGCGCCGCGGTGGGCTGGACATCAACCCGTATCGCAGCACCGGGCCGATCAGCCCGGATAACAAACGGCTGGTCCGGCAGTAATACAGGGGGCCGCTTTGCGGCCCCCGCTATTTCAGATCCCCATGCTGTGCAGCGAGTTGGCAATGCTGCGCAGGGTGGCAGTGAGGTCGGGGTGGTCGGCCTCGAAGCGCTCGATTGCCAGGTTCACGCCGTCCACCAGGTTGTTGTCTGGCGTGGCTTCCT  >  NC_002947/2467810‑2468103

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: