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.

CTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCC‑GCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTT‑GCTCCCATAATGCCCCAGCGCACTAC‑GCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAG‑ATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCAC  >  NC_000913/2666880‑2667375
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                            
cTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCC‑GCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTT‑GCTCCCATAATGCCCCAGCGCACTAC‑GCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAACGCAAAG‑ATAATGACTTCCGCCACCACGCCCCGCGACCACAAATAcccc                                                                                                                                                                                                                                                          >  2:311524/1‑251 (MQ=255)
     cTACCCTGGCGAGCAGCATTATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCC‑GCAATGCAGAATTTGCACCACTCTCAACCACGGCAACGCCGTAGTT‑GCTCCCATAATGCCCCAGCGCACTAC‑GCCGAAACTCGCCGAGATCAACAGCATATCGCGTGCACTACACCGGCGGAAAAAGTTCTTCCTAACCGCAAAG‑CAAATTACTTCCGCCACCACGCCCACCGACCACAAATACCCCA‑cccc                                                                                                                                                                                                                                                    >  2:382335/1‑249 (MQ=255)
      cccccagccgggcagaaagaaaaCGCATGGCGGCCGGGTGGCAGACCTGGAAGGGGCC‑GCAATGCAGAATTTGCCCCCCTATCAGCCACGGCAACGCCGTAGTT‑GCTCCCGTAATGCCCCAGCGCACTAC‑GCCGCAAATCGCCGAGACCAACAGCATATCGCGTTCACGACATCGGCGGAAAAGGTTAGTACGCAGCGCAAAG‑ATAACGACTGCCGCCACCAGGCCCAGCGACCACAAATGCACCA‑ccgcc                                                                                                                                                                                                                                                   <  2:167438/230‑1 (MQ=255)
         ccTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCC‑GCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTT‑GCTCCCATAATGCCCCAGCGCACTAC‑GCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGGTTATTACTAAGCGCAAAG‑ATAATGACTTCCTCCACCACGCCCAGCGACCACAAATACCCCA‑ccgccg                                                                                                                                                                                                                                                  >  1:445783/1‑249 (MQ=255)
                          aaccgcgtgcgggcCAGGTGGCA‑ACCGGGGAAGTGCCCGCAAGGCGGAATTTGCACCCTTACCAACCAGGGCAACGCC‑TAGTTTGCCCCAATAATGCCCCACCGC‑CTACCGCCGCAAATCGCGGAGATCAACACAAGCTCGCTTGCCCTACACCGGCGGAAAAGTTTATTCCTCAGCGCAAAG‑ATAATGCCGCCCGCCACCCCGCCCAGCTCCCACAAAACCCCCA‑CCGCCGAGGCCGAGTA‑GCCAGctgc                                                                                                                                                                                                                              <  2:110945/240‑1 (MQ=255)
                              gcagggcggcCAGGTGGCAGCCCGTAAAGGTGCC‑GCAATTCAAAATTTGCCCCCCTACCAACCACGGCAACGCCGTAGTT‑GCTCCCATAATGCCCCAGCGCACTAC‑GCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTCCTCAGCGCAAAG‑ATAATGACTTCCGCCACCCCGCCCAGCGACCCCAAATCCCCCA‑CCGCCGAGGCCGAGTA‑GCCAGctgcctgc                                                                                                                                                                                                                          <  2:86977/247‑1 (MQ=255)
                                aTGGCGTCAAGGTGGCAGACCGTGAAGGCGCC‑GCAATGCAGAATTTGCACCCCTATCAACCACGGCAAAGCCGTAGTT‑GCTCCCATAATGCCCCAGCGCACTAC‑GCCGAAAATCGCCGAGATCAACAGCATATCGCGGGCACTACCACGGCGGAAAAGGTTATTAATCAGCGCAAAGGATA‑TGACTTTCGCCGCCACGCCCCGCGACCACAAATACCCCA‑CCGCCCAGGCCCAGCA‑GCCacacgcccgccc                                                                                                                                                                                                                        >  2:395352/1‑240 (MQ=255)
                                 gggcgccCAGGTGGCAGACCGGGAAGGTGCC‑GAAATGCAGAATTTGCACCACAACCAACCACGGCAACGCCGTATTT‑GCTCCAATAATGCCCCAGCGCACTAC‑GCCGCAAAGCGCCGAGATCAACAGTATTTCGCGTGACCTACAACGGCGCAAAAGTTTATTACCCAGCGAAAAG‑ATATTGCCTTCCGCCCCCACGCCCAGCGACCACAAATCCCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAg                                                                                                                                                                                                                       <  1:5383/250‑1 (MQ=255)
                                    cggccgggggcagaaccttGAAGGTGCC‑GCAACGCATAATTGCCCCCCCCACCACCCACGCCAACCCCGTAGTT‑GCCCCCATAAGCCCCCAGCCCCCTCC‑GCCGCAAATCCCCGAGACCAACACCATACCGCGTGCCCCACAACGGCGAAAAATTTTATTACCCAGCGCAAAG‑AAAATGACTCCCGCCACCCCGCCCAGCCACCACAAATCCCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTaa                                                                                                                                                                                                                    <  2:143443/233‑1 (MQ=255)
                                        cAGGTGGCAGACCGTGAAGGTGCC‑GCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTT‑GCTCCCATAATGCCCCAGCGCACTAC‑GCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCCGCGCAAAG‑ATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCA‑CCCCCCAGGCCGAGTA‑ACCAGCTGCCTGCCAGTAAATgg                                                                                                                                                                                                                >  2:90947/1‑251 (MQ=255)
                                                           gTGCC‑GCAATGCAGAATTTGCACCACTATAACCCACGGCAACGCCGTAGTT‑GCTCCCATAATGCCCCAGCGCACTAC‑GCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCAATACAACGGCGGAAAAGTTTATTACTCAGCGCCAAG‑CTAATGACTTCCGCCACCACGCCCAGCGACCACAAATACACCC‑CCGCCGCGGCCGAGCA‑GCCCCCTGCATGCCAGTAAATGGCGCTAAAACCGTaactggc                                                                                                                                                                                             >  2:412507/1‑246 (MQ=255)
                                                               c‑GCAATGCAGAATTTGCCCCCCTATCAACCACGGCAACGCCGTAGTT‑GCTCCCATAATGCCCCAGCGCACTAC‑GCCGCAAATCGCCGAGATCAACAGCATACCGCGTGCACTACAACGGCGAAAAAGTTTATTCCTCAGCGCAAAG‑ATAAGGACTTCCGCCACCCCGCCCAGCGACCACAAATACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCa                                                                                                                                                                                         <  1:75345/251‑1 (MQ=255)
                                                                                       aTCAACCACGGCAACGCCGTAGTT‑GCTCCCATAATGCCCCAGCGCACTAC‑GCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAG‑ATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCC‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCCGTAAATGGCGCTAAAAACGTAATAGGCCGCATGTGCCCCCTGCacaaaaccaac                                                                                                                                                                  >  2:258430/1‑241 (MQ=255)
                                                                                       aTCAACCACGGCAACGCCGTAGTT‑GCTCCCATAATGCCCCAGCGCACTAC‑GCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAG‑ATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCtgcctgcc                                                                                                                                                                                                                         >  1:87055/1‑196 (MQ=255)
                                                                                        ccAACCACGGCAACCCCGTAGTT‑GCTCCAATAATTCCCCAGCGCTCTCC‑GCCGCAAACCGCCGATATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACCCAGCGCAAAG‑ATAATGCCTCCCGCCCCCACGCCCAGCGACCACAAATCCCCCC‑CCGCCTAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATaaccaaac                                                                                                                                                                  <  1:144315/249‑1 (MQ=255)
                                                                                          acccACGCCAACGTCGTAGTT‑GCCCCAATACTCCCCCAGCTCACTCC‑TCCGCAAATCGCCGAGAACAACAGCATATCGCGTGCACTACAACTGCGGAAAAGTTTATTACTCAGCGCAAAT‑ATAATGACTTCCGCCACCCCGCCCAGCGACCACAAATCCCCCA‑CCGCCGTGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATaaacaaactca                                                                                                                                                               <  1:335668/249‑4 (MQ=255)
                                                                                          aaCCACGGCAACGCCGTAGTT‑GCTCACATAATGCCCCAGCGCACTAC‑GCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAG‑ATAATGACTTCCGCCACCACCCCCAGCGACCACAAATACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCCGTAAATGGCGCTAAAACCGTCATAGGCCGCATGTTCCCCCTGCAATAAACAAACGc                                                                                                                                                                >  1:366756/1‑250 (MQ=255)
                                                                                                 gCAACGCCGTAGTT‑GCCCCCATAATGCCCCAGCGCACTAC‑GCCGCAATCCCCCGAGATAAACAGCATATCGCGTGCACTACACCGGCGGAAAAGTTTATTACTCAGCGCAAAG‑ATAATGCCTTCCGCCACCACGCCCAGCGACCACAAATACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGa                                                                                                                                                       <  1:463709/251‑1 (MQ=255)
                                                                                                   aaCGCCGTAGTT‑CTTCCCATAATCCCCCACCCCACTAC‑CCCGCAAACCGCCAAGATCAACAGCATATCGCGTGCACTACACCGGCGGAAAAGTTTATTACTCAGCGCAAAT‑ATAATGCCTTCCTCCTCCTCGCCCAGGGACCACAAATACCCCT‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGaaa                                                                                                                                                     <  1:466686/251‑1 (MQ=255)
                                                                                                    ccGCCGAAGTT‑GCTCCCATAATGCCCCGGCGCACTAC‑GCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCCCAACAACGGCGGAAAAGTTTATTACTCAGCGCAAAG‑ATAATGACTTCCGCCACCACGCCCAGCGACCACAAATCCCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAg                                                                                                                                                    <  1:238456/250‑1 (MQ=255)
                                                                                                         gTAGTT‑GCTCCCATAATGCCCCAGCGCACTAC‑GCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAG‑ATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATTTGCCCCCTGCAATAAACAAACGCAGGCC‑CGAAAGCGcc                                                                                                                                                >  1:414381/1‑250 (MQ=255)
                                                                                                            ggtcgcccccAAAATGCCCCGGCGCCCGCC‑GCCGAAAATCGCCGAAAACAAAGGCTGATCGCGTGCACAACACCGGCGGAAAATTTTATTACCCAGGGCAAAT‑ATAATGCCTCCCGCCCCCACCCCCAGCGACCACAAATCCCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAg                                                                                                                                              <  2:401354/243‑1 (MQ=255)
                                                                                                                gCTCCCATAATGCCCCAGCGCACTAC‑GCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAG‑ATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTg                                                                                                                                         >  1:248850/1‑251 (MQ=255)
                                                                                                                gCTCCCATAATGCCCCAGCGCACTAC‑GCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGCTTATTACTCCGCGCAAAG‑ATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCA‑CCGCCGAGGCCGAGCA‑GCCACCTGACTGCCAGAAAATGGCGCCACAACCGTAACAGCCCGCATGTGCCCCCTGCCCTACCCAAACCCAGGCC‑AGCAAGCCCCCGTTCTg                                                                                                                                         >  2:144520/1‑251 (MQ=255)
                                                                                                                gCTCCAATAATGCCCCAGCCCACGAC‑CCCCAAAATCGCCGAGACCACCGGAATATCGCGTGCACTACAACGGCGCAAAAGTTTAGTCCCCGACGCAAAG‑ATACTGACTTCCGCCACCACGCCCAGCGACCACAAATCCCCCA‑CCGCCGAGGCCGCGTA‑GCCAGCTGCCTGCCAGTAAATGGAGCTAAAACAGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCATTTCTg                                                                                                                                         <  2:131792/251‑1 (MQ=255)
                                                                                                                gCTCACATTATGCCACAGCTCACTAC‑GCCGCAAATCGTCGAGTTCAACAGCATACCGCTTTCACTACAACGTCGGAAAATTTTATTCCTCAGCGAAAAG‑ATAACGATTTCCGCCACCACGCCCGCCGACCACAAATCCCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCGGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCt                                                                                                                                          <  1:309918/250‑1 (MQ=255)
                                                                                                                                  cgcACTAC‑GCCGCAAATCGACGAGATCAAAAGTATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAC‑ATAATGACCTCCGCCACCACGCCCAGACACCACAAATACCCCA‑CCCCCGAGCCCGCCTACGCC‑CATGCCCGCCAGTAAATCGCGGTAAAACCGCAATACGACGCACGCGCCCACCCCAATAAACAAACACAGGCCCACAAA‑CGCCAGCACTGGCGCACCAGCGccccac                                                                                                                        >  2:160617/1‑246 (MQ=255)
                                                                                                                                   gcACTAC‑GCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAG‑ATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCAc                                                                                                                      >  1:83446/1‑251 (MQ=255)
                                                                                                                                     aCTAC‑GCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGGTTATTACTCAGCGCAAAG‑ATAATGACTTCCGCCACCCCGCCCAGCGACCACAAATACCCCA‑CCCCCGAGGCAGAGTA‑CACCGCTGCCAGCCATTAAAAGGCGCCCAAACCGTAAGAAGACGCATATGCCCCCCGCCATAAACAAAACCAGGCA‑AAAAACCGCCCGGTCCGGCGAACCAGCGACAACCAccc                                                                                                                    >  2:448674/1‑249 (MQ=255)
                                                                                                                                        aC‑GCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAACGGTTATTACTCAGCGCAAAA‑ATAAAGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCA‑ACGCCGAGGCCCGGTA‑GCCCACTGCCCGCCAGTAAATGGCGCTACAACCGCAATAGGCCGCATGTTCACCCCGCCATAAACAAACCCAGGCC‑ACAAAGCACNCGGTCTGGCGGACCAGGCCCAACCACGCAgc                                                                                                                 >  2:212492/1‑250 (MQ=255)
                                                                                                                                          ggccgcAAACCGCCTAGATCACAAGCATATCGCGTGCATTCCAACGGCTGAAAAGTTTATTACTCATCGCAAAG‑TTAATGCCTTCCCCCACCCCGCCCAGCGACCCCAAATTCCCCA‑CCGCCGATGCCGAGTA‑TCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGAc                                                                                                                <  1:453507/250‑1 (MQ=255)
                                                                                                                                             cgcAAACCGCCGCGACAAACATCATATCGCGTGCCCTCCACCGGGGAAAAATTTTATACCCCAGCGCAAAG‑AAACGGCCTCCCGCCCCCACGCCCACCAACCACAAATCCCCCC‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCaa                                                                                                             <  2:186677/251‑1 (MQ=255)
                                                                                                                                                  aTCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAG‑ATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCCGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCCCCACCCCCGCAAACCACCCggg                                                                                                        >  2:27816/1‑250 (MQ=255)
                                                                                                                                                    cGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAG‑ATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGc                                                                                                      >  1:167583/1‑251 (MQ=255)
                                                                                                                                                       cGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAG‑ATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGctct                                                                                                   >  1:462943/1‑251 (MQ=255)
                                                                                                                                                          gaTCAACAGCATATCGCTTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAG‑ATAATGACTCCCGCCACCACGCCCAGCGACCACAAATCCCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCct                                                                                                <  1:233729/251‑1 (MQ=255)
                                                                                                                                                          gaTCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAG‑ATAATGACTTCCGCCACCACGTCCAGCGACAACAAATACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAAAGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCCCCCACGCAGACCAACCGGTGCTCTCct                                                                                                >  1:114342/1‑251 (MQ=255)
                                                                                                                                                              aaCAGCATATCGCGTGCCCTACAACGGCGGAAAATTTTATTCCACAGCGCAAAA‑ATAATGACTTCCGCCCCCACGCCCAGCGACCACAAATCCCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCctgctg                                                                                            <  1:90947/251‑1 (MQ=255)
                                                                                                                                                                     tatCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAG‑ATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCA‑CCGCCGAGGCCAAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAAAAGGCCGCATGTGCCCCCTGCAATAAACAAAAGCAGGCC‑AGAAAGCGCCAGGTCTGGCGAACCAGCGCCAACCACGCAGAACAACCGGGGCTCTCCTGCTGGCggcct                                                                                     >  1:431058/1‑247 (MQ=255)
                                                                                                                                                                            ggCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAG‑ATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCGGCTGGCGACTTGCCCCtt                                                                              <  1:358329/250‑1 (MQ=255)
                                                                                                                                                                              cacaaCAACGGGGAAAAAGTTTATTCCCAAGCGCAAAG‑ATAATGCCTCCCCCCCCCCCGCCCAGCGACCACAAAACCCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTtgt                                                                            <  2:22354/247‑1 (MQ=255)
                                                                                                                                                                              cACTACAACGGCGAAAAAGTTTATTACCAAGCGAAAAG‑AAAATGCCTCCCGCCACCACGCCCAGCGACCAAAAATCCCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAAGGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCGCTCCTGCTGGCGACTTGCCCCTtgt                                                                            <  2:414381/251‑1 (MQ=255)
                                                                                                                                                                                           gAAAAGTTTATTACCCACCGCAAAG‑AAAAGGACTTCCCCCACCCCCCCCAGCGACCACAAATCCCCCC‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGTCTTGCCCCTTGTGGCTGAAtcgtcg                                                               <  2:299035/251‑1 (MQ=255)
                                                                                                                                                                                                      tcccggacacAAAG‑ATAATGCCTTCCGCACCCGCGCCAAGCCACCACAAATCCCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATgag                                                     <  2:241157/241‑1 (MQ=255)
                                                                                                                                                                                                       aCTCAGCGCAAAG‑ATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCA‑CCGCCGCGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATgaga                                                    >  1:376492/1‑250 (MQ=255)
                                                                                                                                                                                                       aCTCAGCGCAAAG‑ATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATgaga                                                    >  1:201660/1‑250 (MQ=255)
                                                                                                                                                                                                              caaaaG‑ATAATCCCTCCCGCCCCCCCGCCCGGCGCCCACAAAGCCCCCC‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGag                                            <  2:50366/249‑1 (MQ=255)
                                                                                                                                                                                                                       aaTGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATgg                                    >  1:118187/1‑251 (MQ=255)
                                                                                                                                                                                                                          gACTTCCGCCACCACGCCCAGCGACCACAAATACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGAAGCATGGATg                                 >  1:351989/1‑251 (MQ=255)
                                                                                                                                                                                                                           aCTTCCGCCACCCCGCCCAGCGACCACAAATACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTCAAACCGAAATAGGCCGCATGTGCCCCCTGCAATCAACAAACGAAAGCC‑AGAAAGCGCCAGTTATGGCGAACCAGCACCAACCACGCAGCCCAACCGGTGCTCTCCTGCTGCCGACTTGCCCCTTGTCGCCGCATCGTCGGACGGCTGACAACGCCGagccacctggaagc                                >  2:239757/1‑240 (MQ=255)
                                                                                                                                                                                                                                 cccacccCGCCCAGCGACAACAAATCCCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGcc                          <  1:346305/250‑1 (MQ=255)
                                                                                                                                                                                                                                      cacGCCCAGCGACCACAAATACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGCTGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACg                     >  2:163399/1‑251 (MQ=255)
                                                                                                                                                                                                                                         gCCCAGCGACCACAAATACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTc                   >  1:447320/1‑250 (MQ=255)
                                                                                                                                                                                                                                         ccccagggcccACAAATACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCa                  <  2:167583/242‑1 (MQ=255)
                                                                                                                                                                                                                                              gCGACCACAAATCCCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATaaa                                                                                                                                                                       <  2:367850/98‑1 (MQ=255)
                                                                                                                                                                                                                                              gCGACCACAAATACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATaaa                                                                                                                                                                       >  1:367850/1‑98 (MQ=255)
                                                                                                                                                                                                                                                       aaTACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACCGCGCGAGGATc    >  2:197422/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          aCCCCAACCGCCGTGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCAc  <  1:405367/251‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                            
CTTCGCTACCCTGGCGAGCAGCAATATAACGCATGGCGGCCAGGTGGCAGACCGTGAAGGTGCC‑GCAATGCAGAATTTGCACCACTATCAACCACGGCAACGCCGTAGTT‑GCTCCCATAATGCCCCAGCGCACTAC‑GCCGCAAATCGCCGAGATCAACAGCATATCGCGTGCACTACAACGGCGGAAAAGTTTATTACTCAGCGCAAAG‑ATAATGACTTCCGCCACCACGCCCAGCGACCACAAATACCCCA‑CCGCCGAGGCCGAGTA‑GCCAGCTGCCTGCCAGTAAATGGCGCTAAAACCGTAATAGGCCGCATGTGCCCCCTGCAATAAACAAACGCAGGCC‑AGAAAGCGCCAGTTCTGGCGAACCAGCGCCAACCACGCAGACCAACCGGTGCTCTCCTGCTGGCGACTTGCCCCTTGTGGCTGAATCGTCGGACGGATGAGAAAGCCGAGCAGCATGGATGCCACGCCCAACGTCAACAGCGCGAGGATCAC  >  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: ‑