Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,374,0840AG41.0% 36.5 / 67.5 64N288S (AAC→AGC) ycjOputative ABC transporter membrane subunit YcjO
Reads supporting (aligned to +/- strand):  ref base A (18/18);  new base G (5/20);  total (25/39)
Fisher's exact test for biased strand distribution p-value = 3.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGTTTCGTATCGTCACGCTGCCCGCAATTATGCCGGTCCTGGCGACGGTGGTGACACTGCGCACCATCTGGATGTTCTACATGTTCGCGGATGTTTATTTGCTGACGACCAAGGTCGATATTCTCGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCGTTTGTCGACTACTTCCGTAACAGTCTGGTGGTGT  >  NC_000913/1373841‑1374330
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
ggTTTCGTATCGTCACGCTGCCCGCAATTATGCCGGTCCTGGCGACGGTGGTGACACTGCGCACCATCTGGATGTTCTACATGTTCGCGGATGTTTATTTGCTGACGACCAAGGTCGATATTCTCGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGACCCTCAAt                                                                                                                                                                                                                                                 >  2:42040/1‑251 (MQ=255)
   ttcatagcttCACGCTTCCCGCAATTATGCCGGTCCGGGGAACGGTGGTGACCCTGCGCACACTCTTGATTTTCTACATGTTCGCGGCTGTTTATTTGCTGACGACAAAGTACGATATTCTCGGTGTTTATCTCTACAAACCCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATgg                                                                                                                                                                                                                                               <  2:474527/241‑1 (MQ=255)
                gCTGCCCGCAATTATGCCGGTCCTGGCGACGGTGGTGACACTGCGCACCATCTGGATGTTCTACATGTTCGCGGATGTTTATTTGCTGACGACCAAGGTCGATATTCTCGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGc                                                                                                                                                                                                                                 >  1:128807/1‑251 (MQ=255)
                         aaTTATGCCGGTCATGGCGCCGGTGGTGACACTGCGCACCATCTGGATGTTCTAAATGTTCGCGGATGTTTATTTGCTGACGACCAAGGTCGATATTCTCGGTGTCTATCTCTACCAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCCTTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAACCGCCCACTCAgg                                                                                                                                                                                                                        >  2:134935/1‑250 (MQ=255)
                             aTGCCGGTCCTGGCGACGGTGGTGACACTGCGCACCATCTGGATGTTCTACATGTTCGCGGATGTTTATTTGCTGACGACCAAGGTCGATATTCTCGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATAAACGCACACTCAGTCGCa                                                                                                                                                                                                                    <  1:185919/251‑1 (MQ=255)
                              ggCCGGTCCTGGCGACGGTGGTGACACTGCGCACCATCTGGATGTTCTACATGTTCGCGGATGTTTATTTGCTGACGACCAAGGTCGATATTCTCGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATAAACGCACACTCAGTCGCAt                                                                                                                                                                                                                   <  1:328248/250‑1 (MQ=255)
                                    tCCTGGCGACGGTGGTGACACTGCGCACCATCTGGATGTTCTACATGTTCGCGGATGTTTATTTGCTGACGACCAAGGTCGATATTCTCGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGttt                                                                                                                                                                                                             >  1:77297/1‑251 (MQ=255)
                                            gcGGTGGGGACACTGCGCACCAGCTGGATGTTCTACAGGTTCGCGGGGGGGTATGTGCTGCCGCCCAAGGTCGACATGCTCGGTGTATACCTCTACAAAACCTCCTTTGCCTATAATGATTTAGGAAAAGGGGCGGCGATCTCGGTGGTGCTCTTCAGCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAAGCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTGACTGCg                                                                                                                                                                                                     <  2:220661/250‑1 (MQ=255)
                                                    gACACTGCGCACCATCTGGATGTTCTACATGTTCGCGGATGTTTATTTGCTGACGCCCAAGGTCGATATTCTCGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGAACACTCAGTCGCATCGGTTTTTACTGCGGGCTGgcg                                                                                                                                                                                             >  2:304860/1‑251 (MQ=255)
                                                       acTGCGCACCCTCTGGATGTTCTACCTGTTCGCGGATGTTTATTTGCTGACGACCAAGGTCGATATTCTCGGTGTATATCTCTACCAAACCGCCTTTGCATTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCCTTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATCAACGCACACTCAGTCGCATCGGtt                                                                                                                                                                                                              >  2:435360/1‑231 (MQ=255)
                                                       acTGCGCACCATCTGGATGTTCTACATGTTCGCGGATGTTTATTGGCTGACGACCAAGGTCGATATTCTCGGTGTAGATCGCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCGCTTCATCATTATGTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCGCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGtt                                                                                                                                                                                                              <  1:435360/231‑1 (MQ=255)
                                                                     ggATGTTCTACATGTTCGCGGATGTTTATTTGCTGACGACCAAGGTCGATATTCTCGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCTCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCtcatcat                                                                                                                                                                               <  1:70765/248‑1 (MQ=255)
                                                                         gTTCTTCAGGTTCGCGGATGTTTATGTGCTTAAGACCAGGGTCGATATTCTCGGTGTATAGCTCTACAAAACCGCCTTTGCCTTTAAGGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTGCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATAATCGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTAGCGCTGTTTCTCATCATCACGCt                                                                                                                                                                         <  2:60694/250‑1 (MQ=255)
                                                                              aCATGTTCGCGGATGTTTATTTGCTGACGACCAAGGTCGATATTCTCGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAGCGGGTGATCCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTGTTGACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTc                                                                                                                                                                    <  1:48080/250‑1 (MQ=255)
                                                                              aCATGTTCGCGGATGTTTATTTGCTGACGACCAAGGTCGATATTCTCGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGtt                                                                                                                                                                      <  1:402288/248‑1 (MQ=255)
                                                                                  gTTCGCGGATGTTTATTTGCTGACGACCAAGGTCGATATTCTCGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCAtt                                                                                                                                                                <  1:42040/250‑1 (MQ=255)
                                                                                         gATGTTTATTTGCTGACGACCAAGGTCGATATTCTCGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGGAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCAAACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACCCTGTTTCCATTTTTtgtg                                                                                                                                                         >  2:462659/1‑250 (MQ=255)
                                                                                            gTTTATTTGCTGACGACCAAGGTCGATATTCTCGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATg                                                                                                                                                      <  1:89182/250‑1 (MQ=255)
                                                                                                    gCTGACGACCAAGGTCGATATTCTCGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCCTCACGCTGTTTCCCTTTTTTGTGATGCTGATGAcc                                                                                                                                             >  2:139794/1‑251 (MQ=255)
                                                                                                     cTGACGACCAAGGTCGATATTCTCGGTGTATATCTCTACAAAACAGCCTTTGCCTTTAATGATTTAGGAAACGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCt                                                                                                                                            >  1:286953/1‑251 (MQ=255)
                                                                                                           tccAAGGTCGATATTCTCGGTGTATATCTCTACAAAGCCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCGCGTCATCATGATTTTCGCTGTCATTCTGCTGGCCAGGAAACGGGTGAACCTCAATGGCAACAAGTAGACGCACACTCAGTCTCATGGGTGTTTACTGCGGGCTGGCGCTGTTTCCCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGtt                                                                                                                                        <  1:226261/248‑1 (MQ=255)
                                                                                                                ggTCGATATTCTCGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGgc                                                                                                                                 <  1:96199/251‑1 (MQ=255)
                                                                                                                 gTCGATATTCTCGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGTGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGgcg                                                                                                                                <  1:206575/251‑1 (MQ=255)
                                                                                                                    gATATTCTCGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGcgcg                                                                                                                              >  1:12116/1‑250 (MQ=255)
                                                                                                                    gATATTCTCGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATAATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCCTTTTTTGTGATGCTGATGACCTCGTTCAAGGGcgcg                                                                                                                              >  1:221366/1‑250 (MQ=255)
                                                                                                                     atatTCTCGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGAGGCTGATGACCTCGTTCAAGGGCGCGaa                                                                                                                            <  1:114813/251‑1 (MQ=255)
                                                                                                                       atTCTCGCGGTATCTCTCTGAAAAACCGCCTTTCCCATTGATGATATCGGACAAGCGGCGGCGATCTCGGTGGTGCTGTTCATCATTATTTTAGCTGTCATTCTGCTGCCCAGGAAACGCGCGAACCTCAATGCCCACAAATAAAAGCACACTCAGTCGCATCGGTTTTTACTGCCGGCTGGCGCTGTTTCTCATAATCACCCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGaccg                                                                                                                          >  2:280452/1‑248 (MQ=255)
                                                                                                                         tctcGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAaga                                                                                                                         <  1:66914/250‑1 (MQ=255)
                                                                                                                            cGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCg                                                                                                                     >  2:200112/1‑251 (MQ=255)
                                                                                                                            cGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCg                                                                                                                     >  2:256394/1‑251 (MQ=255)
                                                                                                                                        ctACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCCCGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCAt                                                                                                         >  2:239947/1‑251 (MQ=255)
                                                                                                                                        ctACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGCGGTGCTCTTCATCATTAGTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCGCGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGc                                                                                                           <  1:134935/249‑1 (MQ=255)
                                                                                                                                          aCAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATcc                                                                                                       <  1:418097/251‑1 (MQ=255)
                                                                                                                                          aCAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATcc                                                                                                       <  1:304860/251‑1 (MQ=255)
                                                                                                                                             aaaCCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGGTTCGCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGGCCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTAc                                                                                                    <  1:323643/251‑1 (MQ=255)
                                                                                                                                              aaCCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGCGGTGCTCTTCCTCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTACCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTagg                                                                                                   >  2:46475/1‑249 (MQ=255)
                                                                                                                                                ccGCCTGTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGGTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACgct                                                                                                 <  1:239947/251‑1 (MQ=255)
                                                                                                                                                  gcctttgcctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATAAACGCCCACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACgctgc                                                                                               >  1:282793/1‑251 (MQ=255)
                                                                                                                                                         cctttAATGATTTAGGAAAAGCGGCGGCGATCTCGGGGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCgca                                                                                        <  2:136899/251‑1 (MQ=255)
                                                                                                                                                              aaTGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGt                                                                                   <  1:61902/251‑1 (MQ=255)
                                                                                                                                                                             gcggcggcgATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACGGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCAtt                                                                    <  1:325176/251‑1 (MQ=255)
                                                                                                                                                                                 cggcgATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGGTCAAGGGCGCTAAAGAGGCGATCTCACGGCATCCTACGCTGCTGCCGCAGCAGTGGACTCTGGGGCATTACGt                                                                <  2:424371/251‑1 (MQ=255)
                                                                                                                                                                                  ggcgATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCTCAGCAGTGGACTCTGTAGCATTACGTc                                                               <  1:366908/251‑1 (MQ=255)
                                                                                                                                                                                          cGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACAttt                                                        <  2:77297/250‑1 (MQ=255)
                                                                                                                                                                                          cGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACAtt                                                         <  2:401510/249‑1 (MQ=255)
                                                                                                                                                                                            gtggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt                                                                                                                                                          <  2:389056/150‑1 (MQ=255)
                                                                                                                                                                                            gtggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTtgt                                                                                                                                                          >  1:389056/1‑150 (MQ=255)
                                                                                                                                                                                             tggtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGCGCATTACGTCGACATTTTTa                                                     >  2:21732/1‑250 (MQ=255)
                                                                                                                                                                                               gtgCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTTGAGCATTACGTCGACATTTTTAAcc                                                  <  1:392975/251‑1 (MQ=255)
                                                                                                                                                                                                    ctTCATCATTATTTTCGCGGTCATTCTGCTGACCAGGAAACGGGGTAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCga                                               <  1:283932/249‑1 (MQ=255)
                                                                                                                                                                                                            ttattTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGa                                                                         <  2:390544/215‑1 (MQ=255)
                                                                                                                                                                                                            ttattTTCGATGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGTGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGa                                                                         >  1:390544/1‑215 (MQ=255)
                                                                                                                                                                                                              attTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCGt                                   >  1:73309/1‑251 (MQ=255)
                                                                                                                                                                                                              attTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATAAACGCACACGCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTGCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCg                                    <  1:227465/250‑1 (MQ=255)
                                                                                                                                                                                                              attTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCg                                    <  2:337631/250‑1 (MQ=255)
                                                                                                                                                                                                                 ttCGCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCGttt                                 <  1:13094/250‑1 (MQ=255)
                                                                                                                                                                                                                    gCTGTCATTCTGCTGACCAGGAAACGGGTGAGCCGCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCGACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCGTTTGTCg                             <  1:231882/251‑1 (MQ=255)
                                                                                                                                                                                                                                ctgACCAGGAAACGGGTGACCCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCGTTTGTCGACTACTTCCGt                  >  1:110154/1‑250 (MQ=255)
                                                                                                                                                                                                                                      aGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCCTTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCGTTTGTCGACTACTTCCGTAACAGTc           >  1:219494/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTccc                                     <  1:51313/225‑1 (MQ=255)
                                                                                                                                                                                                                                      aGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTccc                                     >  2:51313/1‑225 (MQ=255)
                                                                                                                                                                                                                                         aaaCGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCGTTTGTCGACTACTTCCGTAACAGTCtgg        <  1:21732/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aCGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCGTTTGTCGACTACTTCCGTAACAGTCtggtg      <  1:228476/251‑1 (MQ=255)
                                                                                                                                                                                                                                               gTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCGTTTGTCGACTACTTCCGTAACAGTCtggtggggt  >  2:154371/1‑248 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      
GGTTTCGTATCGTCACGCTGCCCGCAATTATGCCGGTCCTGGCGACGGTGGTGACACTGCGCACCATCTGGATGTTCTACATGTTCGCGGATGTTTATTTGCTGACGACCAAGGTCGATATTCTCGGTGTATATCTCTACAAAACCGCCTTTGCCTTTAATGATTTAGGAAAAGCGGCGGCGATCTCGGTGGTGCTCTTCATCATTATTTTCGCTGTCATTCTGCTGACCAGGAAACGGGTGAACCTCAATGGCAACAAATAAACGCACACTCAGTCGCATCGGTTTTTACTGCGGGCTGGCGCTGTTTCTCATCATCACGCTGTTTCCATTTTTTGTGATGCTGATGACCTCGTTCAAGGGCGCGAAAGAGGCGATCTCACTGCATCCTACGCTGCTGCCGCAGCAGTGGACGCTGGAGCATTACGTCGACATTTTTAACCCGATGATTTTCCCGTTTGTCGACTACTTCCGTAACAGTCTGGTGGTGT  >  NC_000913/1373841‑1374330

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