Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,213,7650AC29.7% 142.3 / 10.2 75V120G (GTG→GGG) yjaBpeptidyl‑lysine N‑acetyltransferase YjaB
Reads supporting (aligned to +/- strand):  ref base A (36/16);  new base C (1/21);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 2.66e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATCCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGT  >  NC_000913/4213518‑4214015
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                           
cGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATCCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCAACATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGTTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCAcc                                                                                                                                                                                                                                                           <  2:505108/250‑1 (MQ=255)
          aCGTTAATGAAATATCGCCAGTTCCACATCCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCCCCTTCTTATAGaa                                                                                                                                                                                                                                               <  1:389668/251‑1 (MQ=255)
              ttacgAAATATCGCCAGTGCCACATCCATGCGCAATCAGCGGTACTCAGTGATCGTGCGGTCGTTGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGTGTCCTAACCCTGAAGTTCAGCCCACCGGGATGAGAAACAATCGTCTACGCCCCCACATACGCCAGATTCAGCACTGGCTACGGTTTCCCCAAATCGTCCCCCTCAGGGCGTCCCGTAACCTTCAACCCCACCTTCGTATAGAAccc                                                                                                                                                                                                                                            <  2:186899/246‑1 (MQ=255)
                aTGAAATATCGCAAGTTCCACATCCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAAc                                                                                                                                                                                                                                         >  2:743783/1‑251 (MQ=255)
                  gAAATATCGCCAGTTCCACATCCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAAcc                                                                                                                                                                                                                                        <  1:437084/250‑1 (MQ=255)
                       atCGCCAGTTCCACATCCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTg                                                                                                                                                                                                                                   <  1:356715/250‑1 (MQ=255)
                               ttCCACACCCATGCGCAATCGGGGGTCCCCAGTGATAGTGGGGCCATGGCAATGCTTAAGCAGAAATAACCGTGTCCCATTTGGGGGGTACTAACCCTGAAGTCCAGCCCACCGGGATGAGAAAAAACCGCCTCCGCCCCCACATACCCCTGATGCAGCAACGAATACGGTTTCCCCAAATGCCCCCCCCCACGCCGTCCCTTAACCTTAAAACCCTCCTACTTATAGAACCCAACCGCCTGCTCAtttt                                                                                                                                                                                                                           <  2:291914/250‑1 (MQ=255)
                                tCCACTTCCATGATCAATCAGCGGTACTCTGTGATAGTGCGTTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATGGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGTATACGGTTTCCCCAAATCTTCCACCTCAGAGCGCCCCGTACCCTTAAAACCCCCCCTCTTATAGACCCCAACCGCCTGCTCATTTTg                                                                                                                                                                                                                          <  1:789238/250‑1 (MQ=255)
                                         aTGCGCAATCCGAGGTACTCAGTGTTAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCACACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCATCAACGGATACGGTATCCCCAAATCGTCCACCTCACAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGACCGCTACTTTTGCTCATTAACg                                                                                                                                                                                                                >  2:611198/1‑251 (MQ=255)
                                                aTCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTACCgttggttg                                                                                                                                                                                                         >  2:677043/1‑251 (MQ=255)
                                                           cacgttaTAGTGCGGTCATGGAAATGCTTAAGAAGAAATAAGCGTGTCCCCATTGGTGGGTACTAAACCTGAAGTTCACCCCACCGGGAGGAAAAAAAATCGCCTTCGCCCCCCCTACCGCCAGGTTCAGCAACGGTTACGGTTTCCCAAAATCGTCCCCCGCAGGGCGTCCCGTAACCTTAAACCCCCCCTTCTTATAGACCCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGt                                                                                                                                                                                              <  2:517520/245‑1 (MQ=255)
                                                                 aTAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATc                                                                                                                                                                                        >  1:190827/1‑251 (MQ=255)
                                                                   aGTGCGGTCATGGAAATGCTTAAGCAGAAATAACCGTGTCCCCATTGGTGGGTACTAACCCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTCCGCCCCCACATCCGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCCCAGAGCGTCCCGTAACCTTAAAACCCCCCTTCTTAGAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCg                                                                                                                                                                                       <  1:731890/250‑1 (MQ=255)
                                                                     tGCGTTCATGGCAATGATTAATCAGAAAAAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGAATGAGAAAAAACCGCCTCCGCCCCCCGATCCGCCAGATTCAGCAACGAATACGTTTTCCCCAAACCGCCCACCTCAGAGCGTCCCGTACCCTTAAACCCCCCCTTCTTATAGACCCCATCCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGTTGCCATCgag                                                                                                                                                                                     <  1:469395/250‑1 (MQ=255)
                                                                       cGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAgc                                                                                                                                                                                  <  1:553455/251‑1 (MQ=255)
                                                                        ggtaaggacAATGTTTAAGAAGAAATAACCGTGTCACCATGGGTGGGTACTAACCCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGCCCCGTAACCTTAAACCCCCCCTTCTTATAGACCCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAgcg                                                                                                                                                                                 <  2:306175/243‑1 (MQ=255)
                                                                        ggTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAgcg                                                                                                                                                                                 >  2:125981/1‑251 (MQ=255)
                                                                        ggTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAgcg                                                                                                                                                                                 >  1:286326/1‑251 (MQ=255)
                                                                            aTGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATg                                                                                                                                                                             >  2:43193/1‑251 (MQ=255)
                                                                                       caaGCAGAAATAATCGTGTCACCATTGTTGGTTACTAACCCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATCCGGTTTCCCCAAATCGTCCACCTCCGAGCGTCCCGTAACCTTAAAACCCCCCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGc                                                                                                                                                                   <  2:192788/249‑1 (MQ=255)
                                                                                         agaaGAAATAATCGTGTCACCATCGGTGGGTACTACACCTGAAGACCAGACCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCCAATCGCCCACCTCAGAGCGTCCCGTCACCTTAAAACCCCCCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTcaccacgcatc                                                                                                                                                                >  2:679540/4‑241 (MQ=255)
                                                                                         aGCAAAAATAACCGTGCCACCATTGGTGGGTACTAAACCGGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTCCCCCCCCACATACGCCAGATTCAGCACCGTATACGGTTTCCCAAAATCGCCCCCCTCAAAGCGCCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCAcc                                                                                                                                                                <  1:707484/251‑1 (MQ=255)
                                                                                            aGAAATAATCGTGTCACCATTGGTGGGTACTCAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCAACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTATTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGATCGCATGCTCCACCAGCACCCg                                                                                                                                                              >  2:82937/1‑250 (MQ=255)
                                                                                            aGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGAGCCATCGAGCGCGCATGCTCCACCAGCACCCg                                                                                                                                                              >  2:857019/1‑250 (MQ=255)
                                                                                             gAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGAc                                                                                                                                                            >  2:631745/1‑251 (MQ=255)
                                                                                                       tgtCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCTAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTAcgccgc                                                                                                                                                   >  1:199129/1‑250 (MQ=255)
                                                                                                          cACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGcc                                                                                                                                               >  1:256168/1‑251 (MQ=255)
                                                                                                            ccATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCcgc                                                                                                                                             >  1:278865/1‑251 (MQ=255)
                                                                                                            ccATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCcgc                                                                                                                                             >  1:56631/1‑251 (MQ=255)
                                                                                                             cagtggtggGTACTAAACCGGAAGTTCAGCCCCCCGGGAGGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTCCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCcgcg                                                                                                                                            <  2:709626/248‑1 (MQ=255)
                                                                                                               ttGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGAATGCTCCACCAGCACCCGACCTACGCCGCAGCCgcgc                                                                                                                                           >  1:304192/1‑250 (MQ=255)
                                                                                                                  ggtggGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGcac                                                                                                                                         >  2:296420/1‑250 (MQ=255)
                                                                                                                  ggggggTACTAACCCTGAAGTTCAGCCCACCGGGAAGAAAAAAAACCGCCTACCCCCCCCCACCCGCCAGATTCAGCAACGGATCCGTTTTCCCCAAATCGCCCCCCCCAGAGCTTCCCGAACCCTTAAACCCCCCCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGcaca                                                                                                                                        <  1:696525/248‑1 (MQ=255)
                                                                                                                    tggGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCCCGCACATc                                                                                                                                      >  2:64871/1‑251 (MQ=255)
                                                                                                                            aaaCCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAACCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATc                                                                                                                                <  1:378951/249‑1 (MQ=255)
                                                                                                                             aaCCTGAAGTCCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCCCATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAACCCCACCTACTTATAGACCCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGAt                                                                                                                             <  1:51006/251‑1 (MQ=255)
                                                                                                                                      ttCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGCTAACCAgcgc                                                                                                                    >  2:669298/1‑251 (MQ=255)
                                                                                                                                           cccACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGTTTTCCCCAAATCGTCCCCCTCAGAGCGCCCCGTACCCTTAAACCCCCCCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCa                                                                                                               <  1:43193/251‑1 (MQ=255)
                                                                                                                                                 gggATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCCCCTTCTTATAGACCCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCcgc                                                                                                                                             <  1:389120/215‑1 (MQ=255)
                                                                                                                                                 gggATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCcgc                                                                                                                                             >  2:389120/1‑215 (MQ=255)
                                                                                                                                                  ggATGAGAAAAAATCGCCTACGCCCCCACATCCGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGCCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATAtgctg                                                                                                        <  1:865451/251‑1 (MQ=255)
                                                                                                                                                   gaggagaAAAAATCGCTGCCGCCCCCACATACCCCATATGCGGCACCGTAACCGTTTTCCCCAAATGGCCCCCCTCGGGGCGTCCGGTACCCTAAAACCCCACCTTCTTATAGACCCCAACCGCCTGCTCATTTGGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCAGGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATAtgctgc                                                                                                       <  2:331378/248‑1 (MQ=255)
                                                                                                                                                   gATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGCGAGCGCCTGCTCCACCAGCACCCGACCTCCGCCGCCGCCGCGCACATCCGGATCGATACAAAGCGCATCCATAttctgc                                                                                                       >  2:216481/1‑251 (MQ=255)
                                                                                                                                                    aTGAGAAAAAATCGACTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGcc                                                                                                      >  1:56528/1‑251 (MQ=255)
                                                                                                                                                    aTGAAAAAAAATCCCCTACGCCCCCACAACCGCAAGAGCCAGCAACGGATACGGTTTCCCCAAACCGTCCACCCCAGAGCCTCCCGTAACCTTAAACCCCCCCTTCTTATAGACCCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGcc                                                                                                      <  2:12047/251‑1 (MQ=255)
                                                                                                                                                       agaaacaagctccgactcccccACATACGCCAAATTCACCAACGGATCCGGTTACCCCAACTCTTCCACCACAGAGGGTCCGGTACCTTTAAACCCCCCATTCTGAAAGACACCAACCGCCTGCTCATTTTGCTCATTATCGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATAGATAAACAGCGCATCCATATTCTGACCAc                                                                                                   <  2:173781/234‑1 (MQ=255)
                                                                                                                                                          aaaaaCCGCCTTCGCCCCACCAAACCCCATATGCAGCAACGGATACGGTTCCCCCAAACCGTCCACCTCAGAGCCCCCCGTAACCTTAAAACCCCCCTTCTTATAAAACCCAACCGCCTGCTCATTTTGCTCATTAACGGTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTa                                                                                                <  2:28495/251‑1 (MQ=255)
                                                                                                                                                             aaTCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGc                                                                                              >  1:346108/1‑250 (MQ=255)
                                                                                                                                                               ccGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGCCCCGTAACCTTAAAACCCCCCTTCTTATAGACCCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAAc                                                                                           <  1:743783/250‑1 (MQ=255)
                                                                                                                                                                cGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGTTTTCCCCAAACCGCCCCCCTCAGAGGGTCCCGTAACCTTAAACCCCCCCTTCTTATAGAACCCACCCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCg                                                                                                                                                                                                <  1:213669/149‑1 (MQ=255)
                                                                                                                                                                cGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCg                                                                                                                                                                                                >  2:213669/1‑149 (MQ=255)
                                                                                                                                                                     aCGCCCCCACATCCGCCAGATTCAGCAACGGATGCGGTTTCCCAAAATCTTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCCCCTTCTTAAAGACCCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATaa                                                                                       <  1:820764/249‑1 (MQ=255)
                                                                                                                                                                           ccACATCCGCCAGATCCAGCAACGGATACGGTTTCCCCAATTCGTCCTCCTCAGAGCTTCCCGTAACCTTAAACCCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAAcc                                                                               <  1:776306/251‑1 (MQ=255)
                                                                                                                                                                               aTACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCt                                                                           >  1:176091/1‑251 (MQ=255)
                                                                                                                                                                                         atccAGCAACGGATACGGTTCCCCCAAATCGTCCACCTCAGAGCGCCCCGTATCCTTAAACCCCACCTTCTTATAGACCCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATaaa                                                                                                                          <  2:250571/192‑1 (MQ=255)
                                                                                                                                                                                         aTTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATaaa                                                                                                                          >  1:250571/1‑195 (MQ=255)
                                                                                                                                                                                            cAGCAACGGATACGTTTCCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCCCCTTCTTATAGACCCAAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCAt                                                               <  1:857019/251‑1 (MQ=255)
                                                                                                                                                                                              gCAACGGATACGTTTTCCCCAAACCGTCCACCTCAGAGCGCCCCGTACCCTTAAACCCCCCCTTCTTATAGACCCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTa                                                             <  1:597190/251‑1 (MQ=255)
                                                                                                                                                                                                  cGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCCTTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCCGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGCCTGATCCCGCTCATTAACcg                                                         >  2:236059/1‑251 (MQ=255)
                                                                                                                                                                                                    gaatcgcttacccaaaaTAGTCCACCTCAGAGCGTCCCGTACCCTTAAACCCCCCCTTCTTATAGACCCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACcgcg                                                       <  2:47482/237‑1 (MQ=255)
                                                                                                                                                                                                    gATACTGTTTCCCCAAATAGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCAGACCTAAGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATACATCCAACCGGCTGAACCAGCTCATTAACcgcg                                                       >  2:708171/1‑251 (MQ=255)
                                                                                                                                                                                                    gATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACcgcg                                                       >  1:381132/1‑251 (MQ=255)
                                                                                                                                                                                                         ggTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGAccc                                                   >  1:459706/1‑250 (MQ=255)
                                                                                                                                                                                                         ggTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGAccc                                                   >  2:211059/1‑250 (MQ=255)
                                                                                                                                                                                                                    aTCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCtt                                       >  1:574093/1‑251 (MQ=255)
                                                                                                                                                                                                                    aTCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCtt                                       >  1:134359/1‑251 (MQ=255)
                                                                                                                                                                                                                       gTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCg                                    <  1:296420/251‑1 (MQ=255)
                                                                                                                                                                                                                        ccccccTCAGAGCGTCCCGTAACCTAAAACCCCCCCTCCTTAAAGACCCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCgg                                   <  1:425673/247‑1 (MQ=255)
                                                                                                                                                                                                                         ccACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATcc                                                                      >  2:738749/1‑215 (MQ=255)
                                                                                                                                                                                                                              tCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCCGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCagga                              >  2:473691/1‑250 (MQ=255)
                                                                                                                                                                                                                               cAGAGCGTCCCGTAACCTTAAACCCCACCTTCTTATAGACCCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCagga                              <  1:64871/249‑1 (MQ=255)
                                                                                                                                                                                                                                  agCGTCCCGTAACCTTAAAACCCCCCTTCTTATAGACCCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAggaaggaa                         <  1:125981/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cccGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGgaacgaac                    >  1:535960/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gtacccTTAAACCCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAg                 <  1:677043/247‑1 (MQ=255)
                                                                                                                                                                                                                                                  aaaaCCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGtcctcc          >  1:548254/1‑250 (MQ=255)
                                                                                                                                                                                                                                                         aCCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCggg  >  2:433303/1‑250 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                           
CGCTGGGGCAACGTTAATGAAATATCGCCAGTTCCACATCCATGCGCAATCAGCGGTACTCAGTGATAGTGCGGTCATGGCAATGCTTAAGCAGAAATAATCGTGTCACCATTGGTGGGTACTAAACCTGAAGTTCAGCCCACCGGGATGAGAAAAAATCGCCTACGCCCCCACATACGCCAGATTCAGCAACGGATACGGTTTCCCCAAATCGTCCACCTCAGAGCGTCCCGTAACCTTAAAACCCACCTTCTTATAGAACCCAACCGCCTGCTCATTTTGCTCATTAACGTTGGTTGTCAGTTCCGGTGCCATCGAGAGCGCATGCTCCACCAGCACCCGACCTACGCCGCAGCCGCGCACATCAGGATCGATAAACAGCGCATCCATATGCTGCCCACTTAGCAACATAAATCCAACCGGCTGATCCCGCTCATTAACCGCGACCCACAACGGCGCTTCCGGCAGGAAGGAACGAACTAGGTCCTCCAGCTCGGT  >  NC_000913/4213518‑4214015

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