Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,579,7570CA28.2% 52.9 / 11.8 39C345F (TGC→TTC) ydeMputative anaerobic sulfatase maturation enzyme YdeM
Reads supporting (aligned to +/- strand):  ref base C (19/9);  new base A (0/11);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 1.45e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTT‑GCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTA‑CTTTAGTAAT‑ACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCGTTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGC  >  NC_000913/1579520‑1580003
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        
aGGATTTTCGCAAGGGAAGATGATGATGCCGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGGCATAAGAAGGGGGCTAAGACTGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTT‑GCAACATCCATAATGTGGTCAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTTCGCAAAATTAAGAAAc                                                                                                                                                                                                                                             >  1:265243/1‑251 (MQ=255)
         gCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTT‑GCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCTCAtt                                                                                                                                                                                                                                    >  1:387042/1‑251 (MQ=255)
           aaGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTT‑GCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGTTTGAAAAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCTCATTGt                                                                                                                                                                                                                                  <  1:405728/251‑1 (MQ=255)
                  gatgatgatGACGAAATACCCGACAGGAATTNATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTT‑GCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTA‑Cttt                                                                                                                                                                                                                          <  2:320484/251‑1 (MQ=255)
                  gatgatgatGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTT‑GCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTA‑Cttt                                                                                                                                                                                                                          <  1:324696/251‑1 (MQ=255)
                   atgatgatgACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTT‑GCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTA‑CTTTa                                                                                                                                                                                                                         <  2:164088/251‑1 (MQ=255)
                             cGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTT‑GCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTA‑CTTTAGTAAT‑ACGAt                                                                                                                                                                                                              <  2:712419/251‑1 (MQ=255)
                              gAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTT‑GCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTA‑CTTTAGTAAT‑ACGATg                                                                                                                                                                                                             >  1:461276/1‑251 (MQ=255)
                                 acacCCGACAGGAATTTATGTCTGTTGCCTTGTTAAGTTCATAAGAAGGAGTCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTT‑GCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTA‑CTTTAGTAAT‑ACTATGCtt                                                                                                                                                                                                          <  1:358519/251‑1 (MQ=255)
                                         cAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTT‑GCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGGTTGACAAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTA‑CTTTAGTAAT‑AAGATGCTTAggcaaac                                                                                                                                                                                                  >  2:117238/1‑246 (MQ=255)
                                              aTTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGGGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTT‑GCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTA‑CTTTAGTAAT‑ACGATGCTTAGGACAACCGcc                                                                                                                                                                                              >  2:656234/1‑250 (MQ=255)
                                                  aTGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTT‑GCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTA‑CTTTAGAAAT‑ACGATGCTTAGGACAACCGCCATTg                                                                                                                                                                                          >  2:12749/1‑250 (MQ=255)
                                                                              aagGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTT‑GCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTA‑CTTTAGTAAT‑ACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGAc                                                                                                                                                              >  2:647501/1‑250 (MQ=255)
                                                                               aggaggCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTATCCATTTGTTTT‑GCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTATACCATGGCGT‑TCATATAA‑GGTACCATGGTTGTATAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTA‑CTTTAGTAAT‑ACGATGCTTAGAACAACCGCCATTGCAGATAGTTTTATATGCACATTGCTGACAt                                                                                                                                                            >  2:264700/1‑251 (MQ=255)
                                                                                ggaggaTAAGAATGGAGTTAAAAAAGGGCGTTATTAATTTTTCTCCATTTGTTTT‑GAAACATCCTAAAGGTGTTAAGGCGGTACTTTGTTCTTAGCTAACTCTACCATGGCGT‑TAAGAAAA‑GGTACAATGGTTAAAAAAG‑GGATTTTA‑TAACCTTCGAAAAAATAAAAACCAGTCTC‑TTGTTTTTCTTTATAAAT‑ACGTTGCTTAGGACATCCGCCATTGCAGATAGGTTTATATGCACATTGCTGACAtt                                                                                                                                                           <  2:241446/251‑1 (MQ=255)
                                                                                             ggAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTT‑GCAACATCCATAATGTGGAAAAGCGGAACTCTGTTCTTAGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTA‑CTTTAGTAAT‑ACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCg                                                                                                                                              >  1:259027/1‑251 (MQ=255)
                                                                                              gAGTTAAAAGAGAGCGTTATTAATTATTCCCCATTTGTTTT‑GCAACATCCATAGTGTGGTAAAGCGGTACTCTTTTCTTAGTTAACTCTACCATGGCGT‑TCATATAA‑GTTACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTTCGAAAAAATAAAAAACAGTCTCATTGTTTA‑TTTTATTAAT‑ACGATGCTTAGGAAAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCGt                                                                                                                                             <  1:552947/251‑1 (MQ=255)
                                                                                                gtaaaaagaggtgtttataatttattCTCC‑TTTGTTTTTGCACCATCCATATGGTGTGAAAGGGGTCTTCTGTTCTTAGCTAACTCAACCATGGCTT‑TAATAAAA‑GGTACCATGTTTAAAAAAA‑GATTTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCGCATTGTTTA‑TTTTAGTAAT‑AGGATGTTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCGttt                                                                                                                                           <  1:87747/230‑1 (MQ=255)
                                                                                                aatAAAAGAGAGCTTTATAAATTTTTCGCCTNTTGTTTG‑ACATCAACAATAAGGGAGTAAAGGGGTACTCTTTTCTGAGTTAACTTTCCAATGGCGT‑CCATAAAA‑GGTTCAATG‑TTGAAAAAAAGGATTTTA‑AACCCTTCGCAAAAATAAAAAACGGCCTCTTTGTTTA‑CTTTA‑TAATTACGATGGTGGAGAAATCCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCGttt                                                                                                                                           <  2:328850/249‑1 (MQ=255)
                                                                                                     aaGAGAGCGTTATTAATTATTCTCCATTTGTTTT‑GCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTA‑CTTTAGTAAT‑ACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCGtttttttt                                                                                                                                      >  1:51458/1‑251 (MQ=255)
                                                                                                                attaattaTTCTCCATTTGTTTT‑GCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTA‑CTTTAGTAAT‑ACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCGTTTTTTTTGCGCTGTCAGt                                                                                                                           >  1:507199/1‑251 (MQ=255)
                                                                                                                             cctttttttttacacctcccttAATGTGGTAAGGGGGTTCTCGGTTCTTAGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGGTTGAAAAAA‑GGTTTTTA‑AACCCTTCGAAAAAAAAAGAAACTGTCTCATTGTTTA‑TTTTAGTAAT‑ACGATTCTTAGGACACCCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCGTTTTTTTTGCGCTGTCAGTTGTACACTGTTc                                                                                                               <  1:103273/230‑1 (MQ=255)
                                                                                                                              aTTTGTTTT‑GCAACATCAATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTA‑CTTTAGTAAT‑ACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCGTTTTTTTTGCGCTGTCAGTTGTACACTGTTCAt                                                                                                             >  1:583659/1‑251 (MQ=255)
                                                                                                                                  ttttt‑GCAACACCAAAAATGTGGTAAAGCGGTATTTTGTTCTTAGCTACCGATACCATGG‑GTGTAATATAA‑GGCACAATGGTTAAAAAAA‑GGATTTTA‑TATCCTCCACAAAAAAAAGAAGCAGTCCCATTTTTTA‑TTTTATTATT‑CCGATGCTTGGGAAAGCCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCGTTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGt                                                                                                          <  1:551672/249‑1 (MQ=255)
                                                                                                                                        caaaaacccATAATGTGTAAAAGCGGTCCTCTGTTTTTAGTTAATTCTCCCATGGTGT‑TCATTTAA‑GGAACAATG‑TTTAAAAAA‑GGGTTTTAATACCCTACGAAAAAAAAAGAAACAGTCTCATTGTTTC‑GTTTATTAAT‑ACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCGTTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAg                                                                                                    <  1:436612/244‑1 (MQ=255)
                                                                                                                                         cAACATCCATAATGAGGTAAAGCGGTCCTATGTTCTTAGCTACCTCTACCATGGCGT‑TCATATAA‑GGAACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAACCAGTCTCATTGTTTA‑CTTTAGTAAT‑ACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCGTTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGt                                                                                                   >  2:103402/1‑251 (MQ=255)
                                                                                                                                               ccATAATTTGGTAAAGGGGTACTCTGTTCTTAATTACCTCTACCAGGGGGT‑TAATAAAA‑GGTACCGTGGTTAAAAAAA‑GTATTTTA‑TACCCTTCAAAAAAATAAAAAACAGTCTCATTGTTTA‑TTTTAGTAAT‑ACGATGCTTAGAACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCGTTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAg                                                                                               <  1:112154/249‑1 (MQ=255)
                                                                                                                                                           aaaGCTGTACTCTTTTTTTTACTAACTCTACCATGGTGT‑TCA‑ATAAAGGTACCGTGGTGAAAAAAA‑GGTTTTTA‑TAACCTTCGAAAAAAAAAAAAACGGTCTCTTTGTTAT‑TTTTGGTAAT‑ACGTTTCTTGGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCGTTTTTTTTGCGCTGTCAGTTGTACACTGGTCATCGTTTTGAGCTCAGATTTATTAATGttt                                                                                 <  2:597011/251‑1 (MQ=255)
                                                                                                                                                             aGCGGTCCTCTTTTCTTAGTTAACTCTCCCAGGGCGT‑TAATAAAA‑GGTACCATGGTTAAAAAAA‑GGTTTTTA‑TACCCTTCAAAAAAAAAAGAAACAGTCTCTTTTTTTA‑CTTTAGAAAT‑ACGATGTTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTc                                                                                                                                                         <  1:199195/177‑1 (MQ=255)
                                                                                                                                                             aGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTA‑CTTTAGTAAT‑ACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTc                                                                                                                                                         >  2:199195/1‑177 (MQ=255)
                                                                                                                                                                      ctGTTCTTAGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTA‑CTTTAGTAAT‑ACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCGTTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGt                                                                       >  1:112015/1‑250 (MQ=255)
                                                                                                                                                                              aGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTA‑CTTTAGTAAT‑ACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTg                                                                                                                                                     >  2:650344/1‑164 (MQ=255)
                                                                                                                                                                              aGCTAACTCTACCAGGGCGT‑TAATAAAA‑GGTACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTCCGAAAAAATAAGAAACAGTCTCATTGTTTA‑CTTTAGTAAT‑ACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTg                                                                                                                                                     <  1:650344/164‑1 (MQ=255)
                                                                                                                                                                                  acctctCCCATGGCTT‑TAATAAAA‑GGTCCAATGGTAGAAAAAA‑GGTTTTTA‑TACCCTTCGAAAAAAAAAAAAACATCCTATTTGTTTA‑CTTTAGAAAG‑ACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCGTTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAg                                                            <  2:303522/247‑1 (MQ=255)
                                                                                                                                                                                                       ataA‑GTTACCATGGTTGAAAAAG‑GGATTTTA‑TACCCTTCGAAAAAAAAAGAAACAGTTTATTTTTTTA‑TTTTAGTAAT‑ACGATGCTGAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCGTTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTaa                                       <  1:571533/250‑1 (MQ=255)
                                                                                                                                                                                                                                  aTTTTA‑TAACCTACGCACAACTAAGGAACAGTCTAATTGTTTA‑CTTTAGTAAT‑AAGATGCTTAGGACCACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCGTTTTTTTTGCGCTGTAAGTTGTACACTGTTCCTCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTATCCATTACTTTCAACAAcc               >  2:289824/1‑249 (MQ=255)
                                                                                                                                                                                                                                   ttttA‑TAACCTTCGCAACAATAAGAAACAGAATCATTGTTTA‑ATTTAGTAAT‑ACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCGTTTTCTTTGAGCTGTCAGTTGTACACTGTTCAACGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTAGTAAATGTCTCCATTACTTTCAACAAcc               >  2:665489/1‑248 (MQ=255)
                                                                                                                                                                                                                                            ccTTCGAAAAAATAAGAAACATTCTCATTGTTTA‑TTTTATTAAT‑ACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCGTTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCtt    <  1:656234/251‑1 (MQ=255)
                                                                                                                                                                                                                                              ttCGCAAAAATAAGAAACAGTCTCATTGTTTA‑CTTTAGTAAT‑ACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCGTTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGc  >  1:688093/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                        
AGGATTTTCGCAAGGGAAGATGATGATGACGAAACACCCGACAGGAATTTATGTCGGGTGCCTTGTTAAGGTCATAAGAAGGAGGCTAAGAATGGAGTTAAAAGAGAGCGTTATTAATTATTCTCCATTTGTTTT‑GCAACATCCATAATGTGGTAAAGCGGTACTCTGTTCTTAGCTAACTCTACCATGGCGT‑TCATATAA‑GGTACCATGGTTGAAAAAA‑GGATTTTA‑TAACCTTCGCAAAAATAAGAAACAGTCTCATTGTTTA‑CTTTAGTAAT‑ACGATGCTTAGGACAACCGCCATTGCAGATAGGTTTATATGCACATTGCTGACATTTCGCTGGAATCCGTTTTTTTTGCGCTGTCAGTTGTACACTGTTCATCGTTTTGAGTTCAGATTTATTAATGTTTCCAATTTTGTACTGTGGATAGACAAAATGGTCGCATTCGTAAATGTCTCCATTACTTTCAACAACCAGATTATCCTTGC  >  NC_000913/1579520‑1580003

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