Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,667,1260AC36.4% 63.6 / 12.6 57Y241D (TAT→GAT) hcaTputative 3‑phenylpropionate transporter
Reads supporting (aligned to +/- strand):  ref base A (30/5);  new base C (0/20);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.05e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCAC  >  NC_000913/2666880‑2667375
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                            
cTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAACGCAAAGATAATGACTTCCGCCACCACGCCCCGCGACCACAAATAcccc                                                                                                                                                                                                                                                          >  2:311524/1‑251 (MQ=255)
     cTACCCTGGCGAGCAGCATTATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTCTCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGAAACTCGCCGAGATCAACAGCATATCGCGTGCACTACACCGGCGGAAAAAGTTCTTCCTAACCGCAAAGCAAATTACTTCCGCCACCACGCCCACCGACCACAAATACCCCAcccc                                                                                                                                                                                                                                                    >  2:382335/1‑249 (MQ=255)
      cccccagccgggcagaaagaaaaCGCATGGCGGCCGGGTGGCAGACCTGGAAGGGGCCGCAATGCAGAATTTGCCCCCCTATCAGCCACGGCAACGCCGTAGTTGCTCCCGTAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGACCAACAGCATATCGCGTTCACGACATCGGCGGAAAAGGTTAGTACGCAGCGCAAAGATAACGACTGCCGCCACCAGGCCCAGCGACCACAAATGCACCAccgcc                                                                                                                                                                                                                                                   <  2:167438/230‑1 (MQ=255)
         ccTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGGTTATTACTAAGCGCAAAGATAATGACTTCCTCCACCACGCCCAGCGACCACAAATACCCCAccgccg                                                                                                                                                                                                                                                  >  1:445783/1‑249 (MQ=255)
                          aaccgcgtgcgggcCAGGTGGCAACCGGGGAAGTGCCCGCAAGGCGGAATTTGCACCCTTACCAACCAGGGCAACGCCTAGTTTGCCCCAATAATGCCCCACCGCCTACCGCCGCAAATCGCGGAGATCAACACAAGCTCGCTTGCCCTACACCGGCGGAAAAGTTTATTCCTCAGCGCAAAGATAATGCCGCCCGCCACCCCGCCCAGCTCCCACAAAACCCCCACCGCCGAGGCCGAGTAGCCAGctgc                                                                                                                                                                                                                              <  2:110945/240‑1 (MQ=255)
                              gcagggcggcCAGGTGGCAGCCCGTAAAGGTGCCGCAATTCAAAATTTGCCCCCCTACCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTCCTCAGCGCAAAGATAATGACTTCCGCCACCCCGCCCAGCGACCCCAAATCCCCCACCGCCGAGGCCGAGTAGCCAGctgcctgc                                                                                                                                                                                                                          <  2:86977/247‑1 (MQ=255)
                                aTGGCGTCAAGGTGGCAGACCGTGAAGGCGCCGCAATGCAGAATTTGCACCCCTATCAACCACGGCAAAGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGAAAATCGCCGAGATCAACAGCATATCGCGGGCACTACCACGGCGGAAAAGGTTATTAATCAGCGCAAAGGATATGACTTTCGCCGCCACGCCCCGCGACCACAAATACCCCACCGCCCAGGCCCAGCAGCCacacgcccgccc                                                                                                                                                                                                                        >  2:395352/1‑240 (MQ=255)
                                 gggcgccCAGGTGGCAGACCGGGAAGGTGCCGAAATGCAGAATTTGCACCACAACCAACCACGGCAACGCCGTATTTGCTCCAATAATGCCCCAGCGCACTACGCCGCAAAGCGCCGAGATCAACAGTATTTCGCGTGACCTACAACGGCGCAAAAGTTTATTACCCAGCGAAAAGATATTGCCTTCCGCCCCCACGCCCAGCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAg                                                                                                                                                                                                                       <  1:5383/250‑1 (MQ=255)
                                    cggccgggggcagaaccttGAAGGTGCCGCAACGCATAATTGCCCCCCCCACCACCCACGCCAACCCCGTAGTTGCCCCCATAAGCCCCCAGCCCCCTCCGCCGCAAATCCCCGAGACCAACACCATACCGCGTGCCCCACAACGGCGAAAAATTTTATTACCCAGCGCAAAGAAAATGACTCCCGCCACCCCGCCCAGCCACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTaa                                                                                                                                                                                                                    <  2:143443/233‑1 (MQ=255)
                                        cAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCCGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCCCCCAGGCCGAGTAACCAGCTGCCTGCCAGTAAATgg                                                                                                                                                                                                                >  2:90947/1‑251 (MQ=255)
                                                           gTGCCGCAATGCAGAATTTGCACCACTATAACCCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCAATACAACGGCGGAAAAGTTTATTACTCAGCGCCAAGCTAATGACTTCCGCCACCACGCCCAGCGACCACAAATACACCCCCGCCGCGGCCGAGCAGCCCCCTGCATGCCAGTAAATGGCGCTAAAACCGTaactggc                                                                                                                                                                                             >  2:412507/1‑246 (MQ=255)
                                                               cGCAATGCAGAATTTGCCCCCCTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATACCGCGTGCACTACAACGGCGAAAAAGTTTATTCCTCAGCGCAAAGATAAGGACTTCCGCCACCCCGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCa                                                                                                                                                                                         <  1:75345/251‑1 (MQ=255)
                                                                                       aTCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCCCCGCCGAGGCCGAGTAGCCAGCTGCCTGCCCGTAAATGGCGCTAAAAACGTAATAGGCCGCATGTGCCCCCTGCacaaaaccaac                                                                                                                                                                  >  2:258430/1‑241 (MQ=255)
                                                                                       aTCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCtgcctgcc                                                                                                                                                                                                                         >  1:87055/1‑196 (MQ=255)
                                                                                        ccAACCACGGCAACCCCGTAGTTGCTCCAATAATTCCCCAGCGCTCTCCGCCGCAAACCGCCGATATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACCCAGCGCAAAGATAATGCCTCCCGCCCCCACGCCCAGCGACCACAAATCCCCCCCCGCCTAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATaaccaaac                                                                                                                                                                  <  1:144315/249‑1 (MQ=255)
                                                                                          acccACGCCAACGTCGTAGTTGCCCCAATACTCCCCCAGCTCACTCCTCCGCAAATCGCCGAGAACAACAGCATATCGCGTGCACTACAACTGCGGAAAAGTTTATTACTCAGCGCAAATATAATGACTTCCGCCACCCCGCCCAGCGACCACAAATCCCCCACCGCCGTGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATaaacaaactca                                                                                                                                                               <  1:335668/249‑4 (MQ=255)
                                                                                          aaCCACGGCAACGCCGTAGTTGCTCACATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACCCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCCGTAAATGGCGCTAAAACCGTCATAGGCCGCATGTTCCCCCTGCAATAAACAAACGc                                                                                                                                                                >  1:366756/1‑250 (MQ=255)
                                                                                                 gCAACGCCGTAGTTGCCCCCATAATGCCCCAGCGCACTACGCCGCAATCCCCCGAGATAAACAGCATATCGCGTGCACTACACCGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGCCTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGa                                                                                                                                                       <  1:463709/251‑1 (MQ=255)
                                                                                                   aaCGCCGTAGTTCTTCCCATAATCCCCCACCCCACTACCCCGCAAACCGCCAAGATCAACAGCATATCGCGTGCACTACACCGGCGGAAAAGTTTATTACTCAGCGCAAATATAATGCCTTCCTCCTCCTCGCCCAGGGACCACAAATACCCCTCCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGaaa                                                                                                                                                     <  1:466686/251‑1 (MQ=255)
                                                                                                    ccGCCGAAGTTGCTCCCATAATGCCCCGGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCCCAACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAg                                                                                                                                                    <  1:238456/250‑1 (MQ=255)
                                                                                                         gTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATTTGCCCCCTGCAATAAACAAACGCAGGCCCGAAAGCGcc                                                                                                                                                >  1:414381/1‑250 (MQ=255)
                                                                                                            ggtcgcccccAAAATGCCCCGGCGCCCGCCGCCGAAAATCGCCGAAAACAAAGGCTGATCGCGTGCACAACACCGGCGGAAAATTTTATTACCCAGGGCAAATATAATGCCTCCCGCCCCCACCCCCAGCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAg                                                                                                                                              <  2:401354/243‑1 (MQ=255)
                                                                                                                gCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTg                                                                                                                                         >  1:248850/1‑251 (MQ=255)
                                                                                                                gCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGCTTATTACTCCGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGCAGCCACCTGACTGCCAGAAAATGGCGCCACAACCGTAACAGCCCGCATGTGCCCCCTGCCCTACCCAAACCCAGGCCAGCAAGCCCCCGTTCTg                                                                                                                                         >  2:144520/1‑251 (MQ=255)
                                                                                                                gCTCCAATAATGCCCCAGCCCACGACCCCCAAAATCGCCGAGACCACCGGAATATCGCGTGCACTACAACGGCGCAAAAGTTTAGTCCCCGACGCAAAGATACTGACTTCCGCCACCACGCCCAGCGACCACAAATCCCCCACCGCCGAGGCCGCGTAGCCAGCTGCCTGCCAGTAAATGGAGCTAAAACAGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCATTTCTg                                                                                                                                         <  2:131792/251‑1 (MQ=255)
                                                                                                                gCTCACATTATGCCACAGCTCACTACGCCGCAAATCGTCGAGTTCAACAGCATACCGCTTTCACTACAACGTCGGAAAATTTTATTCCTCAGCGAAAAGATAACGATTTCCGCCACCACGCCCGCCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCGGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCt                                                                                                                                          <  1:309918/250‑1 (MQ=255)
                                                                                                                                  cgcACTACGCCGCAAATCGACGAGATCAAAAGTATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAACATAATGACCTCCGCCACCACGCCCAGACACCACAAATACCCCACCCCCGAGCCCGCCTACGCCCATGCCCGCCAGTAAATCGCGGTAAAACCGCAATACGACGCACGCGCCCACCCCAATAAACAAACACAGGCCCACAAACGCCAGCACTGGCGCACCAGCGccccac                                                                                                                        >  2:160617/1‑246 (MQ=255)
                                                                                                                                   gcACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCAc                                                                                                                      >  1:83446/1‑251 (MQ=255)
                                                                                                                                     aCTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGGTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCCCGCCCAGCGACCACAAATACCCCACCCCCGAGGCAGAGTACACCGCTGCCAGCCATTAAAAGGCGCCCAAACCGTAAGAAGACGCATATGCCCCCCGCCATAAACAAAACCAGGCAAAAAACCGCCCGGTCCGGCGAACCAGCGACAACCAccc                                                                                                                    >  2:448674/1‑249 (MQ=255)
                                                                                                                                        aCGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAACGGTTATTACTCAGCGCAAAAATAAAGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCAACGCCGAGGCCCGGTAGCCCACTGCCCGCCAGTAAATGGCGCTACAACCGCAATAGGCCGCATGTTCACCCCGCCATAAACAAACCCAGGCCACAAAGCACNCGGTCTGGCGGACCAGGCCCAACCACGCAgc                                                                                                                 >  2:212492/1‑250 (MQ=255)
                                                                                                                                          ggccgcAAACCGCCTAGATCACAAGCATATCGCGTGCATTCCAACGGCTGAAAAGTTTATTACTCATCGCAAAGTTAATGCCTTCCCCCACCCCGCCCAGCGACCCCAAATTCCCCACCGCCGATGCCGAGTATCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGAc                                                                                                                <  1:453507/250‑1 (MQ=255)
                                                                                                                                             cgcAAACCGCCGCGACAAACATCATATCGCGTGCCCTCCACCGGGGAAAAATTTTATACCCCAGCGCAAAGAAACGGCCTCCCGCCCCCACGCCCACCAACCACAAATCCCCCCCCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCaa                                                                                                             <  2:186677/251‑1 (MQ=255)
                                                                                                                                                  aTCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCCGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCCCCACCCCCGCAAACCACCCggg                                                                                                        >  2:27816/1‑250 (MQ=255)
                                                                                                                                                    cGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGc                                                                                                      >  1:167583/1‑251 (MQ=255)
                                                                                                                                                       cGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGctct                                                                                                   >  1:462943/1‑251 (MQ=255)
                                                                                                                                                          gaTCAACAGCATATCGCTTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTCCCGCCACCACGCCCAGCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCct                                                                                                <  1:233729/251‑1 (MQ=255)
                                                                                                                                                          gaTCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGTCCAGCGACAACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAAAGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCCCCCACGCAGACCAACCGGTGCTCTCct                                                                                                >  1:114342/1‑251 (MQ=255)
                                                                                                                                                              aaCAGCATATCGCGTGCCCTACAACGGCGGAAAATTTTATTCCACAGCGCAAAAATAATGACTTCCGCCCCCACGCCCAGCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCctgctg                                                                                            <  1:90947/251‑1 (MQ=255)
                                                                                                                                                                     tatCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCAAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAAAAGGCCGCATGTGCCCCCTGCAATAAACAAAAGCAGGCCAGAAAGCGCCAGGTCTGGCGAACCAGCGCCAACCACGCAGAACAACCGGGGCTCTCCTGCTGGCggcct                                                                                     >  1:431058/1‑247 (MQ=255)
                                                                                                                                                                            ggCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCGGCTGGCGACTTGCCCCtt                                                                              <  1:358329/250‑1 (MQ=255)
                                                                                                                                                                              cacaaCAACGGGGAAAAAGTTTATTCCCAAGCGCAAAGATAATGCCTCCCCCCCCCCCGCCCAGCGACCACAAAACCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTtgt                                                                            <  2:22354/247‑1 (MQ=255)
                                                                                                                                                                              cACTACAACGGCGAAAAAGTTTATTACCAAGCGAAAAGAAAATGCCTCCCGCCACCACGCCCAGCGACCAAAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAAGGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCGCTCCTGCTGGCGACTTGCCCCTtgt                                                                            <  2:414381/251‑1 (MQ=255)
                                                                                                                                                                                           gAAAAGTTTATTACCCACCGCAAAGAAAAGGACTTCCCCCACCCCCCCCAGCGACCACAAATCCCCCCCCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGTCTTGCCCCTTGTGGCTGAAtcgtcg                                                               <  2:299035/251‑1 (MQ=255)
                                                                                                                                                                                                      tcccggacacAAAGATAATGCCTTCCGCACCCGCGCCAAGCCACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATgag                                                     <  2:241157/241‑1 (MQ=255)
                                                                                                                                                                                                       aCTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGCGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATgaga                                                    >  1:376492/1‑250 (MQ=255)
                                                                                                                                                                                                       aCTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATgaga                                                    >  1:201660/1‑250 (MQ=255)
                                                                                                                                                                                                              caaaaGATAATCCCTCCCGCCCCCCCGCCCGGCGCCCACAAAGCCCCCCCCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGag                                            <  2:50366/249‑1 (MQ=255)
                                                                                                                                                                                                                       aaTGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATgg                                    >  1:118187/1‑251 (MQ=255)
                                                                                                                                                                                                                          gACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGAAGCATGGATg                                 >  1:351989/1‑251 (MQ=255)
                                                                                                                                                                                                                           aCTTCCGCCACCCCGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTCAAACCGAAATAGGCCGCATGTGCCCCCTGCAATCAACAAACGAAAGCCAGAAAGCGCCAGTTATGGCGAACCAGCACCAACCACGCAGCCCAACCGGTGCTCTCCTGCTGCCGACTTGCCCCTTGTCGCCGCATCGTCGGACGGCTGACAACGCCGagccacctggaagc                                >  2:239757/1‑240 (MQ=255)
                                                                                                                                                                                                                                 cccacccCGCCCAGCGACAACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGcc                          <  1:346305/250‑1 (MQ=255)
                                                                                                                                                                                                                                      cacGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGCTGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACg                     >  2:163399/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTc                   >  1:447320/1‑250 (MQ=255)
                                                                                                                                                                                                                                         ccccagggcccACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCa                  <  2:167583/242‑1 (MQ=255)
                                                                                                                                                                                                                                              gCGACCACAAATCCCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATaaa                                                                                                                                                                       <  2:367850/98‑1 (MQ=255)
                                                                                                                                                                                                                                              gCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATaaa                                                                                                                                                                       >  1:367850/1‑98 (MQ=255)
                                                                                                                                                                                                                                                       aaTACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACCGCGCGAGGATc    >  2:197422/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          aCCCCAACCGCCGTGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCAc  <  1:405367/251‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                            
CTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCCGCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTTGCTCCCATAATGCCCCAGCGCACTACGCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAGATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCACCGCCGAGGCCGAGTAGCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCCAGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCAC  >  NC_000913/2666880‑2667375

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