Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,348,2260AG73.3% 97.7 / 26.9 60D173D (GAT→GAC) dcuBanaerobic C4‑dicarboxylate transporter DcuB
Reads supporting (aligned to +/- strand):  ref base A (10/6);  new base G (21/23);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.23e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GACTACAGCGATTGCCCCGAGGAAAATCCACATTGCCAGCCAGTTGCTTTTCGGCAGTTTTTTATCCAGCAGCGTCGCGGTATCACCGTAAACATACTCACGGTTTTCCGGTACGGAGATGAATTTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGATCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATCGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTCCGGACGGATGTTGTTCTTAATGGCGACGTCGTAGATGATCGGCAGAATGGTGTAAACCACATGACCCGTACCGCAAAGAATGGTCAGTGTACAGGTCACAAACGGCGCGAC  >  NC_000913/4347998‑4348471
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
gACTACAGCGATTGCCCCGAGGAAAATCCACATTGCCAGCCAGTTGCTTTTCGGCAGTTTTTTATCGAGCAGCGTCGCGGTATCACCGTAAACATACTCACGGTTTTCCGGTACGGAGATGAATTTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGAc                                                                                                                                                                                                                                 <  1:300940/251‑1 (MQ=255)
  cTACAGTGATTGCCCCGAGGAAAATCCACATTGCCAGCCAGTTGCTTTTCGGCAGTTTTTTATCGAGCAGCGTCGCGGTATCACCGTAAACATACTCACGGTTTTCCGGTACGGAGATGAATTTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAAGGTCGTTGGAATGGTGATTGCCAGCAGGTGGAGGAACTCAAGATGGCGAcc                                                                                                                                                                                                                                >  1:110442/1‑250 (MQ=255)
    aCAGCGATTGCCCCGAGGAAAATCCACTGTTCCAGCCAGTTGCTGTTCTGCAGTTTTTTATCGAGCAGTGTCGCGGTATCACCGGAAACATACTCACGGTTTTCCGGTACGGAGATGAATGTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGGTTAGCGTCGATGGAATGGAGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCAt                                                                                                                                                                                                                              <  1:165797/250‑1 (MQ=255)
      ctctgccgatctgaagAAAATCCACATTGCCATCCAGTTGCTTTTCGGCAGTTTTGTATCGAGCAGCGTCGCGGTATCACCGTAAACATACGCACGGTTTGCCGGTACGGAGATGAATTTCTGGAACTCTTCGTCGTTATCCAGAGCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAATGTCGATGGAATGGTGATTGCCAGCAGGGCGAGGAACTCAAGATGGCTACCAta                                                                                                                                                                                                                             <  2:375597/234‑2 (MQ=255)
        cgcttGCCCCGGGGAAACTCCACATTGCCAGCCAGTTGCTTTTCGGCCGTTTTTTATCGAGCAGCGTCGCGGTATCACCGTAAACATCCTCACGGTTTTCCGGTACGGAGATGAATTTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGCTCGCCAGGATACCGATTAACGTCGATGGACTGGGGATTGCCAGCAGGTCGAGGAACTCaa                                                                                                                                                                                                                                          >  2:198028/4‑234 (MQ=255)
        cGATTGCCCCGAGGAAAATCCACATTGCCAGCCAGTTGCTTTTCGGCAGTTTTTTATCGAGCAGCGTCGCGGTATCACCGTAAACATACTCACGGTTTTCCGGTACGGAGATGAATTTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCaa                                                                                                                                                                                                                                          <  1:198028/234‑1 (MQ=255)
         gATTTCCCCGAGGAAAATCAACATTGCCAGACAGTTGGTATTCGGCAGTATTTTATCAAAGAGCGTCGCGGTATCACCGAAAACATACTCACGGTTTTACGGTATGGAGATGAATTTCTGGAACTATTCGTCTTGATCCAGATATTGACCGTGGAACCAGCTGAAGATACCGATCGCCAGGATAACGATTAACGTCGATGGTGTGGTGATTGGCAGCAGAGCGACGAACTCAAGaaggagggcatcaaa                                                                                                                                                                                                                          >  2:346591/1‑235 (MQ=255)
           gtGCCCCGATGAAAATCCACATGTCCAGCCGGTTTCTGTTCGTCAGTTTTTTATCCAGCTGCGTCTCGGTATCTCCGGAAACATACTCACGGTTTTCCGGTACGGAGATGTATTTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGATTCATGTTAAGATGCCGATCTCCAGGATCCCTATTGACGTCGATGGAATGGTGATTGTCATCAGATCTAGGAATTCAAGATTGCTACCATCAACGGTg                                                                                                                                                                                                                      <  2:345307/250‑1 (MQ=255)
             gCCCCGAGGAAAATCCACATTGCCAGCCAGTTGCTTTTCGGCAGTTTTTTATCGAGCAGCGTCGCGGTATCACCGTAAACATACTCACGGTTTTCCGGTACGGAGATGAATTTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGAGACCGATTAACGTCGATGGAATGGTGATTGTCAGCAGGTCGAGGAAGTCAAGATGGCGACCATCAAAGGTGa                                                                                                                                                                                                                     <  2:205198/250‑1 (MQ=255)
                 cGAGGAAAATCCACATTGCCAGCCAGTTGCTTTTCGGCAGTTTTTTATCGAGCAGCGTCGCGGTATCACCGTAAACATACTCACGGTTTTCCGGTACGGAGATGAATTTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACAtt                                                                                                                                                                                                                 >  1:386805/1‑250 (MQ=255)
                 cGAGGAAAATCCACATTGCCAGCCAGTTGCTTTTCGGCAGTTTTTTATCCAGCAGCGTCGCGGTATCACCGTAAACATACTCACGGTGTTCCGGTACGGAGATGAATTTCTGGAACTCTGCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGATCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTa                                                                                                                                                                                                                <  1:235403/251‑1 (MQ=255)
                  gAGGAAAATCCACATTGCCAGCCAGTTGCTTTTCGGCAGTTTTTTATCGAGCAGCGTCGCGGTATCACCGTAAACATACTCACGGTTTTCCGGTACGGAGATGAATTTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATCCCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCAta                                                                                                                                                                                                                             >  1:375597/1‑236 (MQ=255)
                   aGGAAAATCCACATTGCCAGCCAGTTGCTTTTCGGCAGTTTTTTATCGAGCAGCGTCGCGGTATCACCGTAAACATACTCACGGTTTTCCGGTACGGAGATGAATTTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTAcc                                                                                                                                                                                                              <  1:145033/251‑1 (MQ=255)
                                      gccagTTGCTTTTCGGCAGTTTTTTATCGAGCAGCGTCGCGGTATCACCGTAAACATACTCACGGTTTTCCGGTACGGAGATGAATTTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGac                                                                                                                                                                                           <  1:301861/251‑1 (MQ=255)
                                        cagTTGCTTTTCGGCAGTTTTTTATCCAGCAGCGTCGCGGTATCACCGTAAACATACTCACGGTTTTCCGGTACGGAGATGAATTTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGATCGCGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATCGCAACCAGAGacac                                                                                                                                                                                         >  2:154607/1‑251 (MQ=255)
                                        cagTTGCTTTTCGGCAGTTTTTTATCCAGCAGCGTCGCGGTATCACCGTAAACATACTCACGGTTTTCCGGTACGGAGATGAATTTCTGGAACTCCTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGATCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATCGCAACCCGAGacac                                                                                                                                                                                         >  2:250292/1‑251 (MQ=255)
                                         agTTGCTTTTCGGCAGTTTTTTATCGAGCAGCGTCGCGGTATCACCGTAAACATACTCACGGTTTTCCGGTACGGAGATGAATTTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACg                                                                                                                                                                                        >  1:405106/1‑251 (MQ=255)
                                                 ttCGGCAGTTTTTTATCCAGCCGCGTCGCGGTATCACCGTAAACATACTCACGGTTTTCCGGTACGGAGATGAATTTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCCGCAGATCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCGGCATCGCAACCGGAGacacacccgcagc                                                                                                                                                                                >  2:171397/1‑242 (MQ=255)
                                                                 cgagcagcGTCGCGGTATCACCGTACACATACTCACGGTTTTCCGGTACGGAGATGAATTTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGGCGCGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACGCACGCAACCGACACCGGACTGGCg                                                                                                                                                                >  2:214043/3‑251 (MQ=255)
                                                                     cggcGTCGCGGTACAACCGTAAACATACTCGCGGTTTTCCGGTACGGAGATGAATTTCTGGAACTCTTCGTCTTTATCCAGACCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGATCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATCGCAACCAGAGACACAACCGCAACCGTCACCGGACTGTCGATaa                                                                                                                                                            <  2:204695/249‑1 (MQ=255)
                                                                                tATCAACGTAAACATACTCACGGTTTTCCGGTACGGAGATGAATTTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGAAACCAGAGACACGACAGCAACCGACACCGGACTGGCGATAATCCCCATCtgt                                                                                                                                                 >  2:234464/1‑251 (MQ=255)
                                                                                    cccGTAAAGCTACTCACGGTTTTCCAGTACGGAGATGAATTTCTGGAACTCTTCGCCTTTATCCAGATCTTTACCGCGGAACCAGCAGAAGATACAGATAGCACGGATACCGATTAACGTCGATGGAATGATGATTGACAGCAGGTCGAGGAACTCAAGATGACGCCCATCAAAGGTGGCAATACCAAGCATTGCAACCAGAGACACAACCGCATCCGACACAGGACTGGCCATAATCCCCATCtgtgaac                                                                                                                                             >  2:389915/2‑248 (MQ=255)
                                                                                                   aCGGTTTTCCGGTACGGAGATGAATTTCTGGAACTCTTCGTCTTTATCCAGCTCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGACTGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACGACCGCAACAGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTgccg                                                                                                                              >  1:164562/1‑251 (MQ=255)
                                                                                                        tttCCGGTACGGAGATGAATTTCTGGAACTCTTCGTCTTTATCCAGATCGTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAAGGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCAt                                                                                                                          <  2:386805/250‑1 (MQ=255)
                                                                                                                aCGGAGATGAATTTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGCTGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTc                                                                                                                 >  1:33633/1‑251 (MQ=255)
                                                                                                                aCGGAGATGAATTTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGAGAATCCCCGGCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTc                                                                                                                 <  1:239011/251‑1 (MQ=255)
                                                                                                                 cGGAGATGAATTTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTcc                                                                                                                <  1:144032/251‑1 (MQ=255)
                                                                                                                   gagaTGAATTTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGATCGAGGAACTCAAGATGGAGACCATCAAAGGTGACATTACCCAGCATCGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTCcgg                                                                                                              >  1:365680/1‑251 (MQ=255)
                                                                                                                    agatggttttCTGGAACTCTTCGTCTTTATCCAGCTCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGATCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATCGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCtgtg                                                                                                                                                >  1:110202/8‑216 (MQ=255)
                                                                                                                    agaTGAATTTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGATCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATCGCAACCAGAGACACGACCGCAACCGACACCGAACTGGCGATAATCCCCATCtgtg                                                                                                                                                <  2:110202/216‑1 (MQ=255)
                                                                                                                            ttCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTCCGGACGGAtgtt                                                                                                     >  1:148261/1‑251 (MQ=255)
                                                                                                                               tGGAACTCTTCGTCTTTATCCAGATCTTTACAGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGATCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATCGCAACCAGAGACACGACAGCAACCGACCCCGGACTGGCGATAATCACCATATGTGCACCGATAGAACTTGCCGCCATCGGCAGTTCCGGACGGATgttgtt                                                                                                  >  1:207521/1‑251 (MQ=255)
                                                                                                                                ggAACTCTTCGTCTTTATCCAGATCTTTACCGACGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTGACGTCGATGGAATGGCGCTTGCCAGCAGGTCGAGGCACTCAAGATGGCGACCATCAAAGGAGACATTACCCAGCATTGCAACCAGAGACACGACCGCCACCCACACCGAACTGGCGATAATCCCCATCTGTGCACCGATACAACTTGCCGCCCTCGGACGTTCCCGACGGAtgttgcac                                                                                                 >  1:138128/1‑248 (MQ=255)
                                                                                                                                        tCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCCTCGGCCGTTCCGGACGGATGTTGTTATTAATggg                                                                                         >  2:96713/1‑250 (MQ=255)
                                                                                                                                             tttATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGTTGACATTACCCAGCATTGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTCcgg                                                                                                              <  1:178969/225‑1 (MQ=255)
                                                                                                                                             tttATCAAGATCTTTACCGCGGAACCAGCTGCAGATACCGATCGCCAGGATACCGCTTAACGTCGATGGCATGGAGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGAAACGACAGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTCcgg                                                                                                              >  2:178969/1‑225 (MQ=255)
                                                                                                                                                   cAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGCGATTGCCAGCAGGTCGAGGAACTCACGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACGCCCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCAt                                                                                                                          >  1:337367/1‑207 (MQ=255)
                                                                                                                                                   cAGATCTTTACCGCGGAACCAGATGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCATCATTGCAACCAGAGACACGTCCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACCTGCCGCCAt                                                                                                                          <  2:337367/207‑1 (MQ=255)
                                                                                                                                                      aTCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTCCGGACGGATGTTGTTCTTAATGGCGAcgtc                                                                                   <  1:215272/243‑1 (MQ=255)
                                                                                                                                                           tACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTCCGGACGGCTGTTGTTCTTAATGGCGACGTCGTAGATGATCGGc                                                                      >  1:21006/1‑251 (MQ=255)
                                                                                                                                                           tACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACAACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTCCGGACGGATGTTGTTCTTAATGGCGACGTCGTAGATGATCgg                                                                       <  1:284776/250‑1 (MQ=255)
                                                                                                                                                                    aCCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTCCGGACGGATGTTGTTCTTAATGGCGACGTCGTAGATGATCGGCAGAATGgtg                                                             <  1:96713/251‑1 (MQ=255)
                                                                                                                                                                            aaGATACCGATCGCCAGGATTCCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTCCGGACGGATGTTGTTCTTAATGGCGACGTCGTAGATGATCGGCAGAATGGTGTAAACcac                                                     <  1:302359/251‑1 (MQ=255)
                                                                                                                                                                              gATACCGATCGCCAGGATACCGATTAACGTCGATGGCATGGTGATTGCCAGCAGGTAGCGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGCCACGACCGCAACAGACACCGGACTGGCGCTAATCCCCATCTGTGCACCGATAGAACTTGCCGCAATCGGACGTTCAGGCCGGATGTTGTTCTTAATGGAGACTCCGTAGATGAACGGCAGAATGGTGTAAACcaaat                                                   >  2:132112/1‑248 (MQ=255)
                                                                                                                                                                                     aTCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTCCGGACGGATGTTGTTCTTAATGGCGACGTCGTAGATGATCGGCAGAATGGTGTAAACCACATGACCCGt                                            <  1:234464/251‑1 (MQ=255)
                                                                                                                                                                                          cAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGATCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATCGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTCCGGACGGATGTTGTTCTTAATGGCGACGTCGTAGATGATCGGCAGAATGGTGTAAACCACATGACCCGTACCGc                                       >  1:65026/1‑251 (MQ=255)
                                                                                                                                                                                             gATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTCCGGACGGATGTTGTTCTTAATGGCGAc                                                                                      <  1:196239/201‑1 (MQ=255)
                                                                                                                                                                                             gATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGGGACATTACCCAGCATTGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTCCGGACGGATGTTGTTCTTAATGGCGAc                                                                                      >  2:196239/1‑201 (MQ=255)
                                                                                                                                                                                             gATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGATCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATCGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTCCGGACGGATGTTGTTCTTAATGGCGACGTCGTAGATGATCGGCAGAATGGTGTAAACCACATGACCCGTACCGCaa                                     >  1:159049/1‑250 (MQ=255)
                                                                                                                                                                                                   gATTAACGTCGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTCCGGACGGATGTTGTTCTTAATGGCGACGTCGTAGATGATCGGCAGAATGGTGTAAACCACATGACCCGTACCGCAAAGAATgg                              <  1:214043/251‑1 (MQ=255)
                                                                                                                                                                                                            ggATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTCCGGACGGATGTTGTTCTTAATGGCGACGTCGTAGATGATCGGCAGAATGGTGTAAACCACATGACCCGTACCGCTAAGAATGGTCAGCGTAc                     <  2:405106/250‑1 (MQ=255)
                                                                                                                                                                                                            cGATGGAATGGTGATTGCCAGCAGGTCGCGGCACTCAAGCTGGCGACCATCAAAGGAGACATTACCCAGCATTGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCACATCtgtg                                                                                                                                                >  2:36856/1‑128 (MQ=255)
                                                                                                                                                                                                            cGATGGAATGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCtgtg                                                                                                                                                <  1:36856/128‑1 (MQ=255)
                                                                                                                                                                                                                  aaTGGTGATTGCCAGCAGGTCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACGACCGCACCCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTCCGGACGGATGTTGTTCTTAATGGCGACGTCGTAGATGATCGGCAGAATGGTGTAAACCACATGACCCGTACCGCAAAGAATGGTCAGCGTACAGGTc                >  1:392704/1‑250 (MQ=255)
                                                                                                                                                                                                                   aTGGTGATTGCCAGCAGGTCGATGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCACCCAGAGACGCGACCGCAACCGACATCGGACTGGCGATAATTCCCATCTTTGGACCGAGAGAACTGGCCGCCATCGGACGTTGCGGACGGATGTTGTTCTTAATTGCGGCGTCGTAGATGATCGGCAGAGTGGGTTAAACCACATGACGCGTACCGCGAAGAATGGTCAGCGTACATGTcac              <  2:21006/251‑1 (MQ=255)
                                                                                                                                                                                                                   aTGGTGATTGCCAGCAGGGCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACGGCCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTCCGGACGTATGTTGTTCTTAATGGCGACGGCGTAGATTATCGGCAGAATGGTGTAAACCACATGACCCGTACCGCAAAGAATGGTCAGCGTACAGGTcac              <  2:392704/251‑1 (MQ=255)
                                                                                                                                                                                                                    tGGTTATTGCCAGCAGATCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATCGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTCCGGACGGATGTTGTTCTTAATGGCGACGTCGTAGATGATCGGCAGAATGGTGTAAACCCAATAAACCGAACCGCAAAGAATGGTCAATGTACAggaaac              >  2:240050/1‑246 (MQ=255)
                                                                                                                                                                                                                       tGATTGCCAGCAGGTCGAGGCACTCACGATGGCGACCATCAAAGGTGACATTACCCAGCATTGCAACCAGAGACACGACAGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCCCCGATAGAACTTGCCGCCATCGGACGTTCCGGACGGATGTTGTTCTTACTGGCGACGTCGTAGATGATCGGCAGAATGGTGTAAACCACATGCCCCGTACCGCAAAGAATGGTCAGCGTACAGGTCACAAAc          >  2:258823/1‑251 (MQ=255)
                                                                                                                                                                                                                             ccAGCAGATCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATCGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTCCGGACGGATGTTGTTCTTAATGGCGACGTCGTAGATGATCGGCAGAATGGTGTAAACCACATGACCCGTTCCGCAAAGAATGGTCAGTGTACAGGTCACAAACGGcgcg    <  1:171397/251‑1 (MQ=255)
                                                                                                                                                                                                                               agcagATCGAGGAACTCAAGATGGCGACCAGCAAAGGTTACATGACCCAGCATCGCAACCAGAGACCCGACCGCAACCGGTACCGGTCTGTCGATAATCCCCATCTGTGCACCGATAGAACTTTCCGCCAGCGGACGTTCCGGACTGATGTTGTTCTTAATGGCGACGTCGTATATGATCGGCAGAATGGTGTAAACCACAGGACCCGTACCGCAAAGAATGGTCAGTGTACAGGTCACAAACGGCGCGAc  <  2:106958/251‑1 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
GACTACAGCGATTGCCCCGAGGAAAATCCACATTGCCAGCCAGTTGCTTTTCGGCAGTTTTTTATCCAGCAGCGTCGCGGTATCACCGTAAACATACTCACGGTTTTCCGGTACGGAGATGAATTTCTGGAACTCTTCGTCTTTATCCAGATCTTTACCGCGGAACCAGCTGAAGATACCGATCGCCAGGATACCGATTAACGTCGATGGAATGGTGATTGCCAGCAGATCGAGGAACTCAAGATGGCGACCATCAAAGGTGACATTACCCAGCATCGCAACCAGAGACACGACCGCAACCGACACCGGACTGGCGATAATCCCCATCTGTGCACCGATAGAACTTGCCGCCATCGGACGTTCCGGACGGATGTTGTTCTTAATGGCGACGTCGTAGATGATCGGCAGAATGGTGTAAACCACATGACCCGTACCGCAAAGAATGGTCAGTGTACAGGTCACAAACGGCGCGAC  >  NC_000913/4347998‑4348471

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