Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,825,0560GA29.7% 68.0 / 10.5 37intergenic (‑119/+84)ves/spyHutD family protein Ves/ATP‑independent periplasmic chaperone
Reads supporting (aligned to +/- strand):  ref base G (23/3);  new base A (0/11);  total (23/14)
Fisher's exact test for biased strand distribution p-value = 4.26e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCC  >  NC_000913/1824817‑1825302
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
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)
                                                                                                                                                                                                                               ttttttacactggtgtGAGAA‑AAGAAAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATGCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCTGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCTGCATGt                      <  1:450158/235‑1 (MQ=255)
                                                                                                                                                                                                                                   ttacattttttagaaaa‑gaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGc                  <  1:659061/235‑1 (MQ=255)
                                                                                                                                                                                                                                      aaCT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGATTCTTGGCAGGAAGAATTTTAAGCAGATACAGGCATTTTAACTTTTGCCGCCCGACGTTCTCTCAGACCCTTCTCAACAATACCATTTCATTGATTTTTCTCTTTCGGGGCCAGGATGTTGTAAATTTTCTTTTGGGTTTCCATGTCCGCCAgacagttagcgtttc             >  1:671777/1‑237 (MQ=255)
                                                                                                                                                                                                                                        actttttGAGAA‑GAAAAAAAAAAACCGCCGTTCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTgcg            <  1:491178/247‑1 (MQ=255)
                                                                                                                                                                                                                                                agaa‑gaaaaAAAAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGttctt    <  2:253905/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  aa‑gagaaAAGAAAACCGCCGATCCTGTCCACCGAATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAAACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTcc  >  1:197267/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
CAATGCTGGCACGCCAGTAAAAATCACGTTTTGCCGGTGGGAACGTGCAAATTTCGCGCGTTTCACCGGCAGCATTTCGCCACAGATTCACCGACATTTTACGCA‑TATCAAAGTATTCCATACTCGCTCCCGGCCCGTTCTTTTC‑ATTGTTATCGTTCAGGTAACCG‑TGG‑ATGGTTATCGGCAACT‑CGGGTGATTTACGTT‑AGTGGTGATCAGGCTTTTTTTACAACT‑GTTTGAGAA‑GAGAAAAGAAAACCGCCGATCCTGTCCACCGCATTACTGCAAGGTAGTGGACAAGACCGGCGGT‑CTTAAGTTTTTTGGCTGAAAGATTATTCAGCAGTTGCAGGCATTTTACCTTTTGCCGCTGGACGTTCTGTCAGACGCTTCTCAAAATTAGCATTAAATTGCTTTTTCTGTTCCGGCGTCAGGATGTTGTAAATTTTGTTCTGGGTTTCCATGTGCGCCAGCATGTTAGCTTTGCGCTGTTCTTCC  >  NC_000913/1824817‑1825302

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: ‑