Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,375,5280CA30.8% 54.5 / 12.0 40D220E (GAC→GAAdcuDputative transporter DcuD
Reads supporting (aligned to +/- strand):  ref base C (22/5);  new base A (0/12);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 1.58e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATGGTCACCTTATTTCCGACGCTGGTGAGTCTGGGAGTAAGTCGTCTCTCTGCGGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACG  >  NC_000913/3375291‑3375769
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                 
gATGGTCCCCTGATTTCCGACGCTGGGGAGTCTGGGAGTGATTCGTCTCTCTGCGGTGGCAGTTATCGCAACCACGATGGCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTGCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGGCATTATCTCGGTGGCGATCTCCCACCTTTTCGTGCAACGCGCTTTTGACAAAAAAGATaa                                                                                                                                                                                                                                        <  2:366297/249‑1 (MQ=255)
   ggTCACCTTATTTCCGACGCTGGTGAGTCTGGGAGTAAGTCGTCTCTCTGCGGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAat                                                                                                                                                                                                                                    >  2:364389/1‑250 (MQ=255)
   ggTCACCTTATTTCCGACGCGTGTGAGTCTGGGAGGAAGTCGTTGCTCTTCGGTTGCAGTTATCGCAACCACGATGTCCATTGATGGGGGGATGCTGGAAACGAACTTCATTTGTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATTCTTCTTCCACTGCCAGCTTCCTGTCGCCTCTTGCGTCATTATCGCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAat                                                                                                                                                                                                                                    <  2:154243/250‑1 (MQ=255)
            aTTTCCGACGCTGGTGAGTCTGGGAGTAAGTCGTCTCGCTGCGGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCGTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGGGCAACGCGCTTTTGACAAAAAATATAAAAATATCAATCAc                                                                                                                                                                                                                           <  2:265958/250‑1 (MQ=255)
               cccGACGCTGGTGAGTCTGGGAGGTAGTCGTCTCTCTGCTGGTGCAGTTAACGCAACCACGATATTCATTGAGTGGGGGATTCTGGAAACGAACTCCAGTTGTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACGTCTTCCGCTACCAGCTTCCGGTCGCCTCTGGCGTCATTAGCACGGTGGCGAGCTCCCACTTTTTCGTGCAACTCGCTTTGGACAAAAAAGAGAAAAATATCAATCACGAAc                                                                                                                                                                                                                       <  2:314682/250‑1 (MQ=255)
                        ggTGAGTCTGGGAGTAAGTCGTCTCTCTGCGGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAAGATCAATCACGAACAGGCagag                                                                                                                                                                                                               <  1:413626/250‑1 (MQ=255)
                                ggggAGTAAGTCGTCTCTCTGCGGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGTGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTGGCATGGTGAGCTCGGGGTCGATCGCCCTCTTGTTCGTGCAACGCGCTTTTGAGAAAATAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGct                                                                                                                                                                                                      <  1:212550/250‑1 (MQ=255)
                                       aaGTCGTCTCTCTGCGGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGAt                                                                                                                                                                                                >  1:101643/1‑250 (MQ=255)
                                                       ggTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGctct                                                                                                                                                                                >  1:53263/1‑251 (MQ=255)
                                                                     cAACCACGATATCCATTGAGTGGGGGATTCTGGAAACGTACTCCATTTATGATGCCCAGGTAACGGGAATAACAATTGCCACATACTTCTTCCACTACCCGCTTCCGGTCGCCTCTTGCGTCATCATCTCGGTGGCGATCTCCCCCTTTTTCGTGCAACGCGCTTTTGACAAAACAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCAtttt                                                                                                                                                                   >  1:296042/1‑250 (MQ=255)
                                                                           cGATGTCCATTGAGTGGGGGATTCTGGAAACGCCCTCCCTTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATACAAATATCAATCACGAACAGGCAGTGCAAAACGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCtgtg                                                                                                                                                            >  2:220565/1‑251 (MQ=255)
                                                                            gATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCtgtg                                                                                                                                                            >  2:329161/1‑250 (MQ=255)
                                                                            gATGTCCATTGAGTGGGGGATACTGGAAACGAACTCCATTTTTGCTGCTCAGGTAGGGGAAATGAAAATTGCCACATCCTTCTCCCACTCCCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTGCGTGCACCGCGCTTTTGAAAAAAAAGATAAAAATACCAACCACGAACAGGCAGAGCAAAAAGCTCGCTATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGa                                                                                                                                                           <  1:281363/251‑1 (MQ=255)
                                                                               gTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGc                                                                                                                                                        >  1:146304/1‑251 (MQ=255)
                                                                                                    gaaaaCCAACTCCTTTTTCGCTGCCCAGGTAGGGGGAATGAAAATTGCCACATACTTCCCCCCCCCCCAGCTCCCGGTCGCCTCTTGCGTCATTATCGCGGTGGCGATCCCCCACTTTTTCGTGCAACCCGTTTTTGAAAAAAAAGAAAAAAATACCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATgctcggc                                                                                                                                    <  2:91531/248‑1 (MQ=255)
                                                                                                      aaaCGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATgctcggctcg                                                                                                                                 >  1:408796/1‑251 (MQ=255)
                                                                                                           aaCTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTACCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGtt                                                                                                                            >  1:165997/1‑251 (MQ=255)
                                                                                                            accccttttttGCTGCCCAGGTAGCGGGAATGAAAATGGCCACATACTTCTCCCACTCCCAGCTTCCGGTCCCCTCTGGCGTCATTATCGCGGTGGCGATCCCCCCCTTTTTCGTGCAACGCGCTTTTGAAAAAAAAAATAAAAATATCAACCACAAACAGGCAGAGAAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTc                                                                                                                           <  1:348650/245‑1 (MQ=255)
                                                                                                                cATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCATCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTCAATGTGATGCCGTTACTCCTGATGCCCGGCTCGATGTTCCTCg                                                                                                                       >  2:233942/1‑251 (MQ=255)
                                                                                                                         ggccacggaaGCGGAAATGAAAATTGCCCCATTCTTCTTCCCCTCCCGGCTTCCGGTCGCCTCTTGCGTCTTTATCTGGTTGGCTATCCCCCCCTTTTTGTTGCACCCCCTTTTTAAAAAAAAAAATAAAAAAACCACTCAGGACCGGACAGAGAAAAAAGCTCTCGATAATCTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCg                                                                                                              <  1:80285/242‑1 (MQ=255)
                                                                                                                           cccgggtggCGGAAATGAAAATCGCCACATTCTTCTTCCCCCCCCAGCTCCGGTTCCCCTCTTGCCCCATTACCTGGGTGGCGACCCCCCACTTTTTCTTGCGCCGCGCGCTTGCAAAAAAAAAAAAAAAAATCATCCACGAGCAGGCAAAGAAAAAGGCTCTCAAGAATTTCCCTCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCggg                                                                                                            <  1:35660/243‑1 (MQ=255)
                                                                                                                            ccAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGc                                                                                                           >  1:62361/1‑251 (MQ=255)
                                                                                                                                          tGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCCGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCTCCCACGTCGGGCTGATGCAGTccga                                                                                              >  1:253908/1‑247 (MQ=255)
                                                                                                                                           gAAAATTGCCACATACTTTTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCCACGCGCTTTTGACAAAAAAGCTAACAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTACCGCCGCTCTATTAAGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTAGGCTCGCTGTTCATCCCCTACGTCGGGCAGATGCAGTccgcac                                                                                            >  2:299850/1‑246 (MQ=255)
                                                                                                                                            aaaaTTGCCACTTCCTTCTCCCACTCCCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCCCCCCCTTTTTCGTGCAACGCGCTTTTAAAAAAAAAGAAAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACt                                                                                           <  1:364389/251‑1 (MQ=255)
                                                                                                                                            aaaaTTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACt                                                                                           >  2:347854/1‑251 (MQ=255)
                                                                                                                                             aaaTTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCATCTTGCGTCATTATCTCGGAGGCGCTCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGATCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCCGAACTg                                                                                          >  1:277188/1‑251 (MQ=255)
                                                                                                                                                      caTACTTCTTCCACTCCCAGTTTCCGGTCCCCTCTGGCGTCATTAGCTCGGTGGCGATCTCCCATTTTTTCGTGCAACGCGTTTTTGAAAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGAAAAAAGCTCTGGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCtggtg                                                                                 <  1:329084/251‑1 (MQ=255)
                                                                                                                                                          cttcttCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATATCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTg                                                                                                                              >  1:77439/1‑202 (MQ=255)
                                                                                                                                                                   cTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTCGCGATCTCCAACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATACTGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTat                                                                    >  1:222075/1‑250 (MQ=255)
                                                                                                                                                                              cGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACt                                                         >  1:371739/1‑251 (MQ=255)
                                                                                                                                                                               ggTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACtg                                                        >  1:234384/1‑251 (MQ=255)
                                                                                                                                                                                             gTCATTATCTCGGTGGCGATCTCCTACTTTTTCGTGCAACGCGCTTTTGCCAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTc                                                                                                               >  2:264099/1‑183 (MQ=255)
                                                                                                                                                                                             gTCATTATCGCGGTGGCGACCCCCCACTTTTTCTGGCAACGCGTTTTTGAAAAAAAAAAAAAAAAAATCAACCACGACAAGGCAGAGAAAAAAGCTCTCGAAAATTCCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTc                                                                                                               <  1:264099/183‑1 (MQ=255)
                                                                                                                                                                                                tttaccggggTGGCGATCTCCCTCTTTTCCGTGCAACGTGCTTTAGAAAAAAAAGATAAAAATATAAACCACGAACGGGCAAAGAAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTttgttg                                         <  2:62361/242‑1 (MQ=255)
                                                                                                                                                                                                 ttATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGCACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTttgttg                                         >  1:28508/1‑249 (MQ=255)
                                                                                                                                                                                                  taccgcGGTGGCGATCTCCCCTTTTTGCGTCAAACGCGCTTTGGAAAAAAAAAATAAAAATATCAATCCCAAACAGGCAGAGAAAAAGGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGTCTGTGACGATGTTTGTTGAg                                       <  1:176274/245‑1 (MQ=255)
                                                                                                                                                                                                       cGGTGGCGATCTCCCACTTTGTCGTGCAACGCGCTTTTGAAAAAAAAGAAAAAAATATCAACCCCAAACAGGCAGAGAAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGttctt                                  <  1:329161/250‑1 (MQ=255)
                                                                                                                                                                                                                    cctcTTTTTCGTGCAACGCGCTTTTGAAAAAAAAAAAAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTAGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACt                    <  1:110697/248‑1 (MQ=255)
                                                                                                                                                                                                                                        gcTTTTGACAAAAAAGATAAACATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGCTGTTTGTTGCGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACg  >  2:188380/1‑250 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                 
GATGGTCACCTTATTTCCGACGCTGGTGAGTCTGGGAGTAAGTCGTCTCTCTGCGGTGGCAGTTATCGCAACCACGATGTCCATTGAGTGGGGGATTCTGGAAACGAACTCCATTTTTGCTGCCCAGGTAGCGGGAATGAAAATTGCCACATACTTCTTCCACTACCAGCTTCCGGTCGCCTCTTGCGTCATTATCTCGGTGGCGATCTCCCACTTTTTCGTGCAACGCGCTTTTGACAAAAAAGATAAAAATATCAATCACGAACAGGCAGAGCAAAAAGCTCTCGATAATGTCCCGCCGCTCTATTACGCCATTTTACCTGTGATGCCGTTAATCCTGATGCTCGGCTCGCTGTTCCTCGCCCACGTCGGGCTGATGCAGTCAGAACTGCATCTGGTGGTGGTGATGTTACTGAGTTTGACTGTGACGATGTTTGTTGAGTTCTTCCGCAAGCATAACTTGCGCGAAACAATGGACG  >  NC_000913/3375291‑3375769

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