Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,667,1260AC31.8% 97.9 / 10.8 66Y241D (TAT→GAT) hcaTputative 3‑phenylpropionate transporter
Reads supporting (aligned to +/- strand):  ref base A (37/8);  new base C (1/20);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 1.34e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGC  >  NC_000913/2666877‑2667364
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              
tGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATAc                                                                                                                                                                                                                                               >  1:445174/1‑251 (MQ=255)
   cTTCGCTACCCTGGCAAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATAcccc                                                                                                                                                                                                                                            >  1:434744/1‑251 (MQ=255)
      cGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCAcc                                                                                                                                                                                                                                         <  2:162221/251‑1 (MQ=255)
         tACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGCGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGCCTTCCGCCACCACGCCCAGCGACCACAAATACCCCAccgcc                                                                                                                                                                                                                                      >  1:235629/1‑251 (MQ=255)
         tACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCAccgcc                                                                                                                                                                                                                                      >  1:540704/1‑251 (MQ=255)
                    gcagcaATATAACACATGGCGGCCAGGTGGCAGACCGTGACGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCAGCACATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGAATTCCGCCCCCACGCCCACCGCCCACCAATACCCCCCCccccaggccgcga                                                                                                                                                                                                                            >  2:34260/1‑247 (MQ=255)
                     aaacaataaaactgaggcgcgcCAGGTGGCAAACCGTTAAGGGGCCGCAATGCAGAATTTGCACCACTATCACCCACGGTAACTCCGTAGTTGCCCCCATAATGCCCCAGGCCACTACGCCGAAACTCGCCGAGACCAACGGCAAATCGCGTTCCCTCCAACGGCGGAAAATTTTATTACCCAGCGCAAAGATAATGCCGCCCGCCACCACGCCCTGCGCCCACAAATACCCCACCGCCGAGGCCGAGTAg                                                                                                                                                                                                                          <  2:518612/231‑1 (MQ=255)
                         aaTATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGcc                                                                                                                                                                                                                        >  1:108988/1‑249 (MQ=255)
                                cGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCCAAGATAATTACTTCCCCCACCACGCCCAGCGACCACAAATACCCCACCcccgaggccgcgtcacccgccccccg                                                                                                                                                                                                               >  2:435322/1‑234 (MQ=255)
                                  cATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGAAATGCAGAATTTGCACCCCTACCAACCACGGCAACGCCGTAGTTGCCCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATACCGCGTGCACTACAACGGCGGAAAAGTTTATTCCTCAGCGCAAAGATAATGACTTCCGCCACCCCGCCCAGCGCCCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCtgcctgcc                                                                                                                                                                                                             <  1:257475/251‑1 (MQ=255)
                                   aTGGCGGCCAGGGGGACGACCGGGAGGGCGCCGCAATGCAGAATTGGACCCACTACCACCCACGGCAACGCCGAAGTTGCTCCCATAAGGCCCCACCGCCCCCCGCAGCAAACCGCCGAGAGCAACAGCATAGCGCTTGCCCTACACCGGCGGAAAAGTTTATTCCCAAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCa                                                                                                                                                                                                            <  2:261529/251‑1 (MQ=255)
                                         gcCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAAt                                                                                                                                                                                                      >  1:135128/1‑251 (MQ=255)
                                           cAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGGTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCCCCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATgg                                                                                                                                                                                                    >  2:430888/1‑251 (MQ=255)
                                            aGGTGGCAAACCGTGAAGGTGCCGCAATGCAGAATTTGCACCCCTATCACCCACGCCAACGCCGTAGTTGCCCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAATTTTATTACTCAGCGCAAAGATAATGACTCCCGCCACCACGCCCAGCGACCACAAATCCCCCTCCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGgc                                                                                                                                                                                                   <  1:456094/251‑1 (MQ=255)
                                                ggCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAACCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATATGGCCTTCCGCCACCACGCCCAGCGACCAAAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTa                                                                                                                                                                                               <  1:422292/251‑1 (MQ=255)
                                                      ccggtaaGGTGCCGCACTCCCGAATTACCACCACTACCAACCCCGTCAACGCTGTAGTTCTCCCCATAAGCCCCCAGCCCCCTCCGCCGAAAACCGCGAAAACCAACACCATACCTCGTGCACAAAAACGGCGAAAAAGTTAATTACTCAGCGCAAAGAAAATGTCTCCCTCCCCCACGCCCCGCTACAACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTaaatctg                                                                                                                                                                                         <  2:380578/246‑5 (MQ=255)
                                                       cgtaaaGGTGCCCCAATGCAGTATTTGCCCCAATATAAACCCCGGCAACTCCGTAGTTCCCCCCATAATGCCCCAGCGCATTCCGCCGCAAATCGCCGATATCAGCAGCAGACCGCGTGCACTACACCGGCGCAAAAGTTTATTACCCAGAGCAAAGAAAATGACTCCCGCCCCCACGCCTAGCGACCAAAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGt                                                                                                                                                                                        <  2:509866/247‑1 (MQ=255)
                                                        gTGAAGGTGCCGCAATGCAGCATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCCGCGACCACAAAAACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCCGTAAATTGCGCTAAAACCGt                                                                                                                                                                                        >  2:122331/1‑250 (MQ=255)
                                                         tgcaGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAAAGCCGCAGTTGATCCCATAATGCCCCAGCGCACTGAGCCGCAAATCGACGAGATCAACAGCATATCGCGTGCACTACACCGGCGGAAAAGTTTATTACTCAGAGCCAAGGTAATGACTTCCGCAACCACGCCCCGCGACCAAAAATACCCCACCGCCGAGGCCGAGTAGCAAGCTGCCTGCCAGTAAATGGCGCTaaaccccta                                                                                                                                                                                       >  1:565373/4‑244 (MQ=255)
                                                         tGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCCCGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTa                                                                                                                                                                                       >  1:565376/1‑250 (MQ=255)
                                                         tGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGCGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAACACCgca                                                                                                                                                                                       >  1:303423/1‑248 (MQ=255)
                                                         tGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTCAAACCGTa                                                                                                                                                                                       >  1:550821/1‑250 (MQ=255)
                                                           aaGGTGCCGAAATGCAGAATTTGCACCACTACCAACCTCGGCAACCCCGTAGTTGCTCCCAAAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATACCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATa                                                                                                                                                                                    <  1:122331/251‑1 (MQ=255)
                                                             ggTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTCCAACGGCGGAAAAATTTATTACTCAGCACCAACATAATGACTTCCGCCACAACGCCCAGCGACCCCAAATCCCCCCCCGCCGAGGCCGAGTACCCCGCCACCAGCCAGAAAATGGAGCTCAAACCGTAATatg                                                                                                                                                                                  >  2:147492/1‑249 (MQ=255)
                                                              gTGCCGCAATGCAAAATTTGCACCACAATCACCCGCGGCAACGCCGTAGGTGCTCCCAAAAGGCCCCAGCGCATTACGCCGCAAACCCCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTCCTCATCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGCCCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGc                                                                                                                                                                                 <  1:267723/251‑1 (MQ=255)
                                                               ggccccaatgaaaattttccaccacTATCACCCCCGGCAACGCCTTGTTTGCCCCAATAATGCCCCAGCGCACTCCCCCGCAAATCGCCGAGATAAACAGCATACCGCGTGCCCTACAACGGCGGAAAAGTTTATTACCCAGCGAAAAGATAATGCCTCCCGCCACCACGCCCAGCGACCACAAATCCCCCCCCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGcc                                                                                                                                                                                <  2:392373/232‑1 (MQ=255)
                                                                    cAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCAtg                                                                                                                                                                           >  1:495976/1‑251 (MQ=255)
                                                                        gCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACCCCCCGCGACCCCAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAAATGCGCGAAAACCGTAATAGGCCCCATGGGcc                                                                                                                                                                       >  2:116668/1‑251 (MQ=255)
                                                                        gCAGAATTTGCACCACTACCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAACCGCCGAGATCAACAGCATACCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCCCCACGCCCAGCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGcc                                                                                                                                                                       <  1:98028/251‑1 (MQ=255)
                                                                                  caccaccctCAACCACGGCAACGCCGTAGTTGCTCCAATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATACCGCGGGCACTACAACGGCGGAAAAGTTTATTCCTCGGCGCAAAGATAATGACTTCCGCCCCCACGCCCAGCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATa                                                                                                                                                             <  1:225551/243‑1 (MQ=255)
                                                                                  caccacTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCCCCCCGCCCAGCGACCACAAATACCCCCCCGCCCAGGCCGAGTAGCCAGCTGCCTGCCAGTAAACGGCGCCAAAACCGAAAAAAGCCGCATGGTCCCCCCGCAAt                                                                                                                                                              >  2:498034/1‑250 (MQ=255)
                                                                                             aCCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAAAAGGCCGCATGTGCCCCCTGCAATAAACAAACGc                                                                                                                                                    >  1:21561/1‑249 (MQ=255)
                                                                                                  ggCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTAATCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGcc                                                                                                                                               >  1:66683/1‑249 (MQ=255)
                                                                                                        gCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGCAGCCAGCTGCCCGCCAGTAAATGGCCCTAAAACCCAAATAAGCCGCATGTGCCCCCGGCAATAACCAAACCCCGGCCAGAAAgcg                                                                                                                                       >  2:9427/1‑251 (MQ=255)
                                                                                                                 gcacccATAATGCCCCAGCGCCCTACGCCGCAAAGCGCCGAGAACAACAGCATATCGCGTGCACTAAAACGGCGGAAAATTTTATTACTCAGCACAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAACACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCt                                                                                                                               <  1:283117/247‑1 (MQ=255)
                                                                                                                  cccccAAAATGCCCCAGCGCCCTCCGCCCCAATCGCCGAGAACAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGc                                                                                                                            <  1:435322/249‑1 (MQ=255)
                                                                                                                   tCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGc                                                                                                                            >  1:12505/1‑251 (MQ=255)
                                                                                                                         tAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAAcc                                                                                                                       >  1:79753/1‑251 (MQ=255)
                                                                                                                                    cgcACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCCTGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAAccc                                                                                                            >  1:89618/1‑250 (MQ=255)
                                                                                                                                                              cAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCctgct                                                                                  >  1:446202/1‑251 (MQ=255)
                                                                                                                                                                aCAGCAAAGCGCGTGCACTACAACGGCGAAAAAGTTTATTACTCAGCGAAAAGATATTGCCTTCCGCCACCCCGCCCAGCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTgg                                                                                <  1:283172/251‑1 (MQ=255)
                                                                                                                                                                    cATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGAc                                                                            >  2:52183/1‑251 (MQ=255)
                                                                                                                                                                          gtggcacctaaaccggcgaAAAATTTTATTACCAAGCGAAAAGATAAGTACTCCCGCCACCGCGCCCAGCGACCAAAAAACCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGccc                                                                      <  2:550821/238‑1 (MQ=255)
                                                                                                                                                                          gcgTGCACGACAACGGCGGAAAAGTTTATTCCCCAGCGAAAAGATAATGACTTCCGCCCCCACGCCCAGCGACCCCAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGccc                                                                      <  1:171615/251‑1 (MQ=255)
                                                                                                                                                                           ggTGCACTACAACGGCGAAAAATTTTATTACTCAGCCAAAAGATAAGGCCTCCCCCCCCCCCGCCCAGCACCCACAAACCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGcccc                                                                     <  2:12505/250‑1 (MQ=255)
                                                                                                                                                                               cACTACAACGGCGGAAAATTTTATTACTCAGCGAAAAGATAATGCCTTCCGCCACCCCGCCCAGCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTtgt                                                                 <  1:429339/251‑1 (MQ=255)
                                                                                                                                                                               cACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTtgt                                                                 >  1:117146/1‑251 (MQ=255)
                                                                                                                                                                                      aCGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGCGTAGCCAGCTGCCTGCCAGTAAATTGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACCCAGACCAACCGGTGCTCTCCTGCTGGCGACCTTCCCCTTGTGGCTGa                                                           >  2:148020/1‑250 (MQ=255)
                                                                                                                                                                                       cggcggAAAAGTTTATTACCCAGCGCAAAGATAATGCCTCCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAAt                                                         <  1:452332/251‑1 (MQ=255)
                                                                                                                                                                                             aaaaGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTcgg                                                   >  1:465150/1‑251 (MQ=255)
                                                                                                                                                                                                aGTTTATTACCCAGCGCAAAAAAAATGCCTCCCGCCCCCACGCCCAGCGACCAAAATCACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTcggacg                                                <  2:433087/251‑1 (MQ=255)
                                                                                                                                                                                                    tattaCTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTCATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATAGTCGGACGGATg                                            >  1:285450/1‑251 (MQ=255)
                                                                                                                                                                                                     tttccTCAGCGCAAAAAAAATGCCTCCCGCCCCCACGCCCAGCGCCAACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATg                                            <  2:355597/246‑1 (MQ=255)
                                                                                                                                                                                                          tCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGaa                                        >  2:70697/1‑249 (MQ=255)
                                                                                                                                                                                                                aaaaaaaaaggcttccccccaccaCCCCCGGCGACAACAAAACCCCCCCCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGagc                                <  2:21561/233‑1 (MQ=255)
                                                                                                                                                                                                                    aataataACTCCCGCCACCACCCCCAGCGACCACAACTCCCCCGCCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAgcagca                            <  1:274783/250‑1 (MQ=255)
                                                                                                                                                                                                                      tAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATg                          >  1:180848/1‑251 (MQ=255)
                                                                                                                                                                                                                              cccGCCCCCACGCCCAGCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCa                   <  1:150516/249‑1 (MQ=255)
                                                                                                                                                                                                                               ccGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACg                 >  1:220307/1‑251 (MQ=255)
                                                                                                                                                                                                                               ccGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACg                 >  1:477194/1‑251 (MQ=255)
                                                                                                                                                                                                                                      cccccccagcccccACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACg          <  2:540704/239‑1 (MQ=255)
                                                                                                                                                                                                                                      cacGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCCGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAACCGCGCAGCCGTGATGCCCCGCCCAAcc          >  2:191116/1‑250 (MQ=255)
                                                                                                                                                                                                                                       acGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTATATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCTAGTTCTGGCAAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGt         >  1:180504/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGCAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCCGCATGGATGCCACGCCCAACGTc        >  1:325910/1‑250 (MQ=255)
                                                                                                                                                                                                                                              ccacccaaaaaTCCCCCCCCGCCGAGGCCGAGTAGCCAGCTGCCGGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAgc  <  2:252696/243‑1 (MQ=255)
                                                                                                                                                                                                                                                       aaTACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGc                                                                               >  1:326562/1‑165 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              
TGACTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGC  >  NC_000913/2666877‑2667364

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