Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947955,1490AG30.9% 110.5 / 65.2 85intergenic (+103/‑210)cyoE/alaCprotoheme IX farnesyltransferase/aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (27/31);  new base G (12/14);  total (39/46)
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.51e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATT  >  NC_002947/955002‑955294
                                                                                                                                                           |                                                                                                                                                     
gACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACaa                                                                                                                                                       <  1:70711/148‑1 (MQ=255)
gACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAg                                                                                                                                                      >  6:115434/1‑149 (MQ=255)
gACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAg                                                                                                                                                      >  8:110954/1‑149 (MQ=255)
    gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                  <  1:178054/149‑1 (MQ=255)
    gCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTc                                                                                                                                                  <  6:59229/149‑1 (MQ=255)
      cAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                                <  2:130788/149‑1 (MQ=255)
         aCCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTg                                                                                                                                               >  4:200383/1‑147 (MQ=255)
         aCCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTa                                                                                                                                             <  7:110954/149‑1 (MQ=255)
               gCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGc                                                                                                                                       <  3:91131/149‑1 (MQ=37)
                cGGATGTGCTGATGACCTACGCACGCTGATACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCgg                                                                                                                                    <  7:369384/149‑1 (MQ=255)
                cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                    <  6:390430/149‑1 (MQ=37)
                cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                    <  4:193851/149‑1 (MQ=38)
                cGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                    <  7:68262/149‑1 (MQ=38)
                  gATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGt                                                                                                                                  <  3:252166/149‑1 (MQ=37)
                  gATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGt                                                                                                                                  <  3:200383/149‑1 (MQ=255)
                   aTGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTTCAAGGTTCCTGTAGGAGCGGGtt                                                                                                                                 >  5:278431/1‑149 (MQ=255)
                      tgCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTATTTTTATTTTTCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCgg                                                                                                                                    <  3:157075/148‑1 (MQ=18)
                      tgCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTAc                                                                                                                             <  1:292744/149‑1 (MQ=255)
                            tgaCCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGa                                                                                                                       <  8:206521/149‑1 (MQ=255)
                                cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAagag                                                                                                                   <  3:103106/149‑1 (MQ=255)
                                cTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTACTGTGCCGGCCTCTTCGCGGGTAAAACCCCGCCAAAAAGGTTGCAGTAGGTGCGGGTTTACCCGGGAagag                                                                                                                   >  8:41795/1‑149 (MQ=255)
                                                    gcggttttCACCAACACCCCGCCCCCTGTGCCGGGGTTTTTTTGGGGGGTTT‑‑‑‑‑TTTTCCTGGGCCGCCTTCTTCGGGGTTAAACCCGCTCATACAAGGTTCCTGTCGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaataaaaa                                                                                              <  8:214040/144‑1 (MQ=25)
                                                    cccTTTTTCACAAAACCCCCGGCCCGTGGGGGGGGTTTTTTTTTGGGGGTTTT‑‑‑‑‑TTTTCCTGGGCCGCCCTCTCCGCGTCAAACCCCGCTCAAACATGTTCCTTGTCGCCGGGGTTTTACCCGCGAAGAGGCCAGCACAGaaaaataaaa                                                                                               <  8:331184/148‑1 (MQ=255)
                                                     ccTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaatcagg                                                                                               >  5:171222/1‑144 (MQ=255)
                                                     ccTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGaaaaataaaaa                                                                                              >  5:393561/1‑149 (MQ=255)
                                                       tGTTTCACGATAACGCCGTCCACTGTGCCGGGGTTTTTTTTTGGTGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAATGTTCCt                                                                                                                                                <  8:185795/100‑1 (MQ=37)
                                                       tGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCt                                                                                                                                                >  7:185795/1‑100 (MQ=255)
                                                             aCGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTTAACCCGCAGCTAACAGGTTGCGGTAGGAGCGGGGTTTACCGCGAAGAGGGCCGCCCAGAAAAATAAAAATACCAGAt                                                                                     >  8:18783/1‑149 (MQ=255)
                                                               gAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTAccc                                                                                                                           >  7:42004/1‑110 (MQ=21)
                                                               gAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTAccc                                                                                                                           <  8:42004/110‑1 (MQ=21)
                                                                   aCCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGGCCTCTTCGCGGGTAAACCCGCTCCGACAAGGTTCCGGTAGGAGCGGGTTTTCCCGCGAAGAGGCCAGGACAGAAAAGTAAAAATAACAGAttttg                                                                                 >  1:94684/1‑147 (MQ=37)
                                                                     cccGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTAGCGGGTTAAACCGCAACTAAAAAGGTACGGTTTGAGCCGGGTTTACCGCGAAGAGGCCCCCACAGAAAAAAAAAAAAAAATGAATTTaaaa                                                                              >  2:212701/1‑148 (MQ=255)
                                                                            aCTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGt                                                                                                                                  >  1:183175/1‑91 (MQ=17)
                                                                            aCTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGt                                                                                                                                  <  2:183175/91‑1 (MQ=17)
                                                                                   cGGGTTTTTTTTTTGGTGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTATTCGCGTGTAAACCCGCTCCTACAAGGTTCCTGTAGTAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGaaaa                                                               <  5:115434/149‑1 (MQ=37)
                                                                                   cGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGTTAACCCCGCTCCTACAAGGTTCCTGTTGGTGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAGTAAAAATACCTGATTTTCAAAAAATAGATCTGaaaa                                                               <  2:94684/149‑1 (MQ=38)
                                                                                    gggtTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTTTGCCGGCCTCTTCGCGGTTAAACCCGCTCCAAAAAGTTTCCTTTAGTAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCCTGATTTTCAAAAAATAGATCTGaaaa                                                               <  1:212701/149‑1 (MQ=37)
                                                                                         tttttttttGGGGTTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGTTTAAACCCGATCCTAAAAGTTTCATTTAGTAGCGTTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTc                                                           <  1:350494/148‑1 (MQ=37)
                                                                                          ttttttttGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCCCGCCGACAAGGTTCGTGTAGGAGCGGGGTTTACCCGGGAGAGGGCAGCACAGAAAAATAAAAATAACTGATTTTAAAAAAAAAGATGAGAAAAtttgt                                                          >  8:164609/1‑146 (MQ=25)
                                                                                           tttttttGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTCGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCta                                                         <  6:393561/148‑1 (MQ=255)
                                                                                              ttttGGTGTTTT‑‑‑‑‑TTTTCCTGTCCGGGCTGTTTCGGGGTAAACCCGCCTCTAACATGTTTCTGTTTGGTCGGGTTTAACCCGCGAAGAGGTCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCtaata                                                      <  2:334230/148‑1 (MQ=255)
                                                                                                  gggggTTT‑‑‑‑‑TTTTCCGTTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAATGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTc                                                           <  7:362125/139‑1 (MQ=255)
                                                                                                   ggtgTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTTTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACa                                                 <  8:362436/145‑1 (MQ=255)
                                                                                                   gggtgaaT‑‑‑‑‑TTTTCCTTTGCCGGCCTCTTCGCTGTAAAACCCTATCCAAAAATTTTCATTTAGTAGCGTTTTAACCCGCTAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAAc                                                  <  5:239903/142‑1 (MQ=255)
                                                                                                   gggtTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTTTCCTGTAGTAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACa                                                 <  3:130161/148‑1 (MQ=255)
                                                                                                    gggTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGTTTCCTTTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATaaa                                                    <  3:315176/144‑1 (MQ=255)
                                                                                                    gggTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAAACCGCTACTAAAAAGTTGCAGTACGAGCGGGTTTACCCCCCAAGAGGCCACCACAACAAAAAAAAAAAACCAGATTTTCAAAAAATATATAAGAAAAGTTATAaaaaaaaaa                                                >  4:315176/1‑145 (MQ=19)
                                                                                                                     cTGTGCCGGCCTCTTCGCGGGTTAACCCCCTCCTACAAAGTTACGGTTGGAGCCGGTTTTACCGCGAAGAGGCCCGCACAGAAAAATAAAAATAACTGATTTTGAAAAAATAGATATGAAAAGTTATAATAAAAAAGGAAATTATACtt                                    >  4:318948/1‑149 (MQ=255)
                                                                                                                     cTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAAAAGGTTCCTGTAGGATCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCt                                     <  7:18783/148‑1 (MQ=255)
                                                                                                                          ccGGCCTCTTCGCGGTTAAACCCGCTCATACAATTTTCATTTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAg                               <  5:328298/149‑1 (MQ=38)
                                                                                                                             ggCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAgg                              >  6:302673/1‑148 (MQ=255)
                                                                                                                                     cgcgGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCt                        <  5:302673/146‑1 (MQ=255)
                                                                                                                                      gcgGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGttt                    >  2:219305/1‑149 (MQ=255)
                                                                                                                                      gcgGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGttt                    >  7:317770/1‑149 (MQ=255)
                                                                                                                                        gggTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATc                          <  5:75339/141‑1 (MQ=255)
                                                                                                                                        gggTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATc                          >  6:75339/1‑141 (MQ=255)
                                                                                                                                        gggTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCg                  <  1:219305/149‑1 (MQ=255)
                                                                                                                                           tAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCtaat                                                       <  6:49300/109‑1 (MQ=21)
                                                                                                                                           tAAACCCGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCtaat                                                       >  5:49300/1‑109 (MQ=21)
                                                                                                                                           tAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATc                                       >  2:211502/1‑125 (MQ=255)
                                                                                                                                           tAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATc                                       <  1:211502/125‑1 (MQ=255)
                                                                                                                                               cccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAt           >  3:60302/1‑149 (MQ=37)
                                                                                                                                                ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTGCCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGAATTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc          >  8:294081/1‑149 (MQ=37)
                                                                                                                                                ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATa                                                                         >  6:381874/1‑86 (MQ=17)
                                                                                                                                                ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATa                                                                         <  5:381874/86‑1 (MQ=17)
                                                                                                                                                ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCtaa                                                        >  8:308346/1‑103 (MQ=21)
                                                                                                                                                ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCtaa                                                        <  7:308346/103‑1 (MQ=21)
                                                                                                                                                ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGc                             <  8:188615/130‑1 (MQ=37)
                                                                                                                                                ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGc                             >  7:188615/1‑130 (MQ=37)
                                                                                                                                                ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc          >  2:86423/1‑149 (MQ=38)
                                                                                                                                                ccGCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtc          >  8:208071/1‑149 (MQ=38)
                                                                                                                                                ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCATCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCAc              <  5:155396/145‑1 (MQ=255)
                                                                                                                                                ccGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCAc              >  6:155396/1‑145 (MQ=255)
                                                                                                                                                 cGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTAtct         >  1:189160/1‑149 (MQ=255)
                                                                                                                                                  gCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGc                             >  3:192966/1‑128 (MQ=37)
                                                                                                                                                  gCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGc                             <  4:192966/128‑1 (MQ=37)
                                                                                                                                                  gCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGc                >  3:261141/1‑141 (MQ=37)
                                                                                                                                                  gCTCCTACAGGAACCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGc                <  4:261141/141‑1 (MQ=37)
                                                                                                                                                      cTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGAATCCTGTTTCGGCAGTATCTTCTGa    >  4:14162/1‑149 (MQ=255)
                                                                                                                                                       tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt   >  2:280930/1‑149 (MQ=255)
                                                                                                                                                       tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt   >  1:234888/1‑149 (MQ=255)
                                                                                                                                                       tACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAt   >  1:229313/1‑149 (MQ=255)
                                                                                                                                                        aCAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGc                >  6:260468/1‑135 (MQ=255)
                                                                                                                                                        aCAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGc                <  5:260468/135‑1 (MQ=255)
                                                                                                                                                        aCAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCAc              >  6:106693/1‑137 (MQ=255)
                                                                                                                                                        aCAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCAc              <  5:106693/137‑1 (MQ=255)
                                                                                                                                                          aaGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGAtt  >  1:227449/1‑147 (MQ=255)
                                                                                                                                                           |                                                                                                                                                     
GACAGCCAGACCGCTGCGGATGTGCTGATGACCTACGCACGCTGAAACTGCTGCCTGTTTCACGAAAACCCCGGCCACTGTGCCGGGGTTTTTTTTTGGGGGTAT‑‑‑‑‑TTTTCCTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGTTCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGAAAAATAAAAATACCTGATTTTCAAAAAATAGATCTGAAAAGTTCTAATAAAACAGGAAATTATCCTTTACAGGCATCCTGTTTCGGCACTATCTTCTGATT  >  NC_002947/955002‑955294

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: