Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,825,0310CT27.0% 68.6 / 11.7 37intergenic (‑94/+109)ves/spyHutD family protein Ves/ATP‑independent periplasmic chaperone
Reads supporting (aligned to +/- strand):  ref base C (23/4);  new base T (0/10);  total (23/14)
Fisher's exact test for biased strand distribution p-value = 2.87e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AACAAAGAAAACTCCCCATTCGCCGCGATGGATGCAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCC  >  NC_000913/1824783‑1825266
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              
aaCAAAGAAAACTACCCCTTCGCCGCGATGGATGCAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGTTTTACGTT‑AGTGGTGATCAGGCtt                                                                                                                                                                                                                                              >  2:478225/1‑251 (MQ=255)
                     gccgcGATGGATGCAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTgagaa‑g                                                                                                                                                                                                                       <  1:510646/251‑1 (MQ=255)
                      ccgcGATGGATGCAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTgagaa‑ga                                                                                                                                                                                                                      <  1:52767/251‑1 (MQ=255)
                                  cAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTAGCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAAcc                                                                                                                                                                                                          <  2:202278/251‑1 (MQ=255)
                                      gCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑CTTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGGGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACACCT‑GTTTGAGAA‑GAGAAAAGAAAAccgccg                                                                                                                                                                                                      >  2:534361/1‑251 (MQ=255)
                                                                 ggcgggTGGAAACGTCAAAATTTCGCGCTTTTCCCCGGCATCATTCCCCCAAAGATTCCCCG‑TATTTTTCCCAATATCAAAGTTTTCCTTCCCCCCCCCCGGCCTTTTTTTTTC‑TTTGTTATCGTTCAGGTACCCG‑TGG‑TTGTTTATCGCCAACG‑CGGGTGATTTACGTT‑AGTGGTGATCAGTTTTTTTTTACACTT‑TTTTGAGAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGt                                                                                                                                                                           <  1:167584/247‑1 (MQ=255)
                                                                        gggAACGCGCAAATTTCGCGTTTTCCCCCGGCAGTATTCCCCCACAGTTTCCCCGAAATTTTACGAA‑TAACAAAGTTTCCCATCCCACCTCCCGCCCTTTTCTTTTC‑TTTTTTATCGTTCAGGTACCCG‑TGG‑AGGTTTATCGGCAACT‑CGGGTTTTTTACGTT‑AGTGGTGATCAGTTTTTTTTTACAACT‑TTTTGAGAA‑GAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAc                                                                                                                                                                    <  1:562478/251‑1 (MQ=255)
                                                                                 aaaaTTTCGCGCGTTTCACCGGCAGTTTTTCCCAAAATTTTCCCCGACTTTTTACGCA‑TATAAAAGTATCCCATCCTCGCCCCCGCCCCGTTCTTTTA‑TTTGTTATCGTTCAGGTACCCG‑GGG‑TTGTTTAGCGGAACCT‑GGGGTTTTTTACGTT‑AGTGGTGATCGGTTTTTTTTTACACCT‑GTTTGAGAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggc                                                                                                                                                           <  1:395290/250‑1 (MQ=255)
                                                                                  aaaTTTCGCGGTTTTCACCGGAAGCTTTTCCCCACAGATTCACCGACTTTTTACGCA‑TATAAAAGTTTTCCATACCCCCCCCCGCCCCTTTTTTTTC‑TTTTTTATCTTTC‑GGTAACCG‑GGGGAGGTTTATGGGAAACT‑GGGGTTATTTACGTT‑AGGGGTGATCAGGTTTTTTTTACAACT‑GTTTGAAAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcg                                                                                                                                                          <  1:396121/251‑1 (MQ=255)
                                                                                   aaTTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcgg                                                                                                                                                         >  1:253905/1‑251 (MQ=255)
                                                                                          cgcgTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGc                                                                                                                                                                                >  1:163749/1‑221 (MQ=255)
                                                                                               ttCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGTTCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACcggcgggtttcaggtttt                                                                                                                                             >  2:491178/1‑239 (MQ=255)
                                                                                                      gcagctTTCCGCAACAGATCCCCCGACTTTTTACGCA‑TATAAAAGATTTCAATCCCCCCCCCGGGCCTGTTTTTTTA‑ATTTTTAGCGTTCAGGTATCCG‑TGG‑ATGTTTACCGG‑AACTTGGGGTGATTTCCTTT‑AGTGGTGATCAGGTTTTTTTTAAAACT‑GTGTGAGAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGa                                                                                                                                     <  1:612000/251‑1 (MQ=255)
                                                                                                       cagcaTTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGa                                                                                                                                     >  1:495536/1‑250 (MQ=255)
                                                                                                        agcaTTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGGTTTTTGGCTGaa                                                                                                                                    >  2:4772/1‑250 (MQ=255)
                                                                                                          caTTGCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGATAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGAAAAGACCGGCGGA‑Cttcaatagttttgctgaa                                                                                                                                    >  2:163293/1‑232 (MQ=255)
                                                                                                                    caGATTCACCGACATTTTCCGCA‑TATCACAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAACGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑TTTTAATTTTTTTGCtttaaaatttatctgcc                                                                                                                       >  2:462701/1‑235 (MQ=255)
                                                                                                                    caGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTcagcc                                                                                                                       >  1:306897/1‑250 (MQ=255)
                                                                                                                             cGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTa                                                                                                                                                                                    >  1:388664/1‑182 (MQ=255)
                                                                                                                                 aTTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGG‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCAttt                                                                                                          >  1:576540/1‑251 (MQ=255)
                                                                                                                                       cctaataaaaaagtttcccATACCCTCTCCCGGCCCGTTTTTTTC‑TTTGTTATCGTTCAGG‑AACCGGTGT‑AGTGTTATCGGCAACG‑GGGGTGTTTTACGTT‑AGTGGTAATCGGGTTTTTTTAACACCT‑TTTTGAGAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTAcc                                                                                                      <  1:390729/233‑1 (MQ=255)
                                                                                                                                        cacataacaagatATTCCACACTCGCTCCCGGCCGGTT‑TTTTCTATTGTTATTGTTAAGGTAACCG‑TGG‑AGGGTTATCGGCAACT‑TGGGTGATTTCCGTT‑AGTGGTAAACGGTCTTTTTTTACACTT‑GTTGGAGAA‑AAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCtt                                                                                                    <  1:176367/239‑1 (MQ=255)
                                                                                                                                                  aGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTgg                                                                                          >  1:701850/1‑251 (MQ=255)
                                                                                                                                                              tCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTTAGCAGTTGCAGGCATTTTACCTTTTTCCGCTGGACGTTCTGtt                                                                                >  2:358900/1‑248 (MQ=255)
                                                                                                                                                                ccccccGGCCTGTTTTTTTT‑ATTGTTATCGGTCAGGTAGCCG‑GGG‑TGGGTTATCGGAAACG‑CGGGTGTTTTA‑GTTGAGTGGTGATCGGGTTTTTTTAACAACT‑GTTGGAAAA‑GAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGAc                                                                            <  2:306897/248‑1 (MQ=255)
                                                                                                                                                                      ggCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAAAT‑CGGGTGTTTTACGTT‑AGTGGTGCTCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAACGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGG‑CTTCAGTTTTTTGGCTGAAAGATTATTCAGCAGTTTCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCCAACGCTtct                                                                      >  2:318542/1‑251 (MQ=255)
                                                                                                                                                                         ccGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCa                                                                    >  1:111507/1‑250 (MQ=255)
                                                                                                                                                                                tttC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGc                                                            >  1:94797/1‑251 (MQ=255)
                                                                                                                                                                                 tttcatgtttttcgtcaggttatccG‑TGG‑GTGTTTATCGCCACCT‑GGGGTGATTTACTTT‑GGGGGTGATCATGTTTTTTTTAAACCG‑GTTTGAGAA‑AAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAg                                                             <  1:478225/227‑1 (MQ=255)
                                                                                                                                                                                   c‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTTGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCAtt                                                         >  1:45009/1‑251 (MQ=255)
                                                                                                                                                                                                taaGGTACCCG‑TGG‑AGGTTTATCGGAATCT‑GGGGTGTTTGACGTA‑AGTGGTGATCGGGTTTTTTTTACA‑GT‑GTTTGAAAAAGAAAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCtttt                                              <  1:458847/248‑1 (MQ=255)
                                                                                                                                                                                                     gtaacggggg‑ATGGTTATCGGCATGT‑GGGGTGTTTTCGGTT‑GGTGGTGATCAG‑TTTTTTTTGAAACTTGTTGGAGAA‑AAAAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGt                                         <  1:4772/243‑1 (MQ=255)
                                                                                                                                                                                                                   gTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGCA‑GAGAAAAGAAAAACGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTTGCTGAAAGATTATTGAGCAGGTGCAGGGCTTTTTCCTTTTGCCGCTGGGCGGTCTGTCTAACGCTTATCAAAATTAGCATTAAATTGCTTTTTCTTTTTCGGctgcaggg                            >  1:301170/1‑244 (MQ=255)
                                                                                                                                                                                                                         ggCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCCAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGAtgttgt                      >  1:554632/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aTTTACGTT‑AGAGGTGATCCGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAAACGCCGATCCTGTCCACCGCATTACTGCCAGGTAGTGGACAAGACCGGCGGT‑CTTACGTTTTTTTGCTGAAAGATTCTTCAGCAGTTGCAGGGATTTTATCTTTTGTCGCTGGATGTTTTGTCAGACGCTTCTCAAAATTAGACTTAACTTGCTTTTTCTGTTCCGGCGTCCGGGTGTTGTAGATTTTTTTCTg         >  2:629076/1‑251 (MQ=255)
                                                                                                                                                                                                                                          tACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGAGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGAAGGTGCAGGCATTTTACCTTTTGCCGCTGGACGGTCTGTCAGACGATTCTCAAAATTAGCAATAAATTGCTTTTTTTGTTCCGGCGCCAGGGTGTTGTAAATTTTGTTCTGGGt      >  2:172322/1‑251 (MQ=255)
                                                                                                                                                                                                                                              tt‑AGTGGTGCTGAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAACCCGACGATCCTGTCCACCGCATTACTGCAAGGTAGAGGACAAGAACGGCGGAACTAAAGTTTTT‑GGGTGAAAGCTTATTCAACACTTGCAGGCCTTTTACCTTTTTCCCCTGGACGATCTGTCCAACGCTTCCCAAAATTAGCATTTAATTTTTTTTTCTGTTTCGGCGTCAGGGTGTTGTAAATTTTGTTCTggggtttc  >  2:287800/1‑246 (MQ=255)
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                              
AACAAAGAAAACTCCCCATTCGCCGCGATGGATGCAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCC  >  NC_000913/1824783‑1825266

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