Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,628,598	0	G	A	37.5%	80.4 / 23.6	56	S114L (TCA→TTA)	yfgF	cyclic‑di‑GMP phosphodiesterase, anaerobic
Reads supporting (aligned to +/- strand):  ref base G (21/14);  new base A (1/20);  total (22/34)
Fisher's exact test for biased strand distribution p-value = 3.98e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGG‑AAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAAT  >  NC_000913/2628450‑2628747                                                                                                                                                        |                                                                                                                                                        gAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGG‑AAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGa                                                                                                                                                        >  2:983306/1‑150 (MQ=255)       aaaaGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGG‑AAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACgttg                                                                                                                                                 >  1:752656/1‑151 (MQ=255)        aaaGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGG‑AAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACgttgt                                                                                                                                                <  2:120385/151‑1 (MQ=255)         aaGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGG‑AAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACgttgtt                                                                                                                                               >  2:1092782/1‑151 (MQ=255)                   acCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGG‑AAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCgggg                                                                                                                                     >  2:1134811/1‑150 (MQ=255)                    ccGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGG‑AAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGa                                                                                                                                    <  1:594166/150‑1 (MQ=255)                    ccGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGG‑AAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAg                                                                                                                                   <  1:955880/151‑1 (MQ=255)                       aCCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGG‑AAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTa                                                                                                                                 >  2:887727/1‑150 (MQ=255)                        ccAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGG‑AAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTaaa                                                                                                                               <  1:92983/151‑1 (MQ=255)                         cAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGG‑AAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTAAAg                                                                                                                              <  2:731786/151‑1 (MQ=255)                          aGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGG‑AAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTAAAGa                                                                                                                             <  1:1006576/151‑1 (MQ=255)                               ttCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGG‑AAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTAAAGACCCt                                                                                                                         <  2:474341/150‑1 (MQ=255)                                   ctGCTCGCCAGTAATCCTACAAATGCAGCAAACTGG‑AAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTAAAGACCCTG‑Ta                                                                                                                     <  1:884915/149‑1 (MQ=255)                                         gCCAGTAATCCTACAAATGCAGCAAACTGG‑AAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTAAAGACCCTG‑TAACCGcccc                                                                                                             >  1:869924/1‑151 (MQ=255)                                            aGTAATCCTACAAATGCAGCAAACTGG‑AAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAAc                                                                                                          >  2:1064329/1‑151 (MQ=255)                                                   caaaaaaTGCAGAAACCGGA‑AAAATTATCAGAAACAGCGTTGAAGAAAACACAATCTGCAAGAATAGCCGCGGGGAAATTAACCGGG‑TATCACC‑ATGGAAAACGTTGTTA‑CGCCGGGGAGAAAAGCCCCTG‑TAACCGCCCCAACAGAGaa                                                                                                    <  1:176852/146‑1 (MQ=255)                                                   cTACAAATGCAGCAAACTGG‑AAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGaat                                                                                                   >  2:635009/1‑151 (MQ=255)                                                     aCAAATGCAGCAAACTGG‑AAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCCCC‑ATGTAAAACGTTGTTA‑CGCCGGGGAGTAAAGCCCCTG‑TAACCGCCCCAACAGAGaatg                                                                                                  <  2:715729/150‑1 (MQ=255)                                                     aCAAATGCAGCAAACTGG‑AAAAGTATCAGAAACAGCGTTGCAGAAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGT‑TATCACC‑ATGTAAAACGTGGTTA‑CGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGaatga                                                                                                 <  1:33701/151‑1 (MQ=255)                                                       aaaTGCAG‑AAACTGGAAAAATTATCAAAAAAAGCGTGGCAGAAAACACAATCTGCAAGATTATCCGCTGGGAAATAAAACGGT‑TATCCCC‑ATGTGAAACGTTGTTA‑CGCGGGGGAGAAAAGCCCCTG‑AAACCGCCCCAACAGAGaatgata                                                                                               <  1:395275/151‑1 (MQ=255)                                                          tGCAGCAAACTGG‑AAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGaatgataatg                                                                                            >  2:103387/1‑151 (MQ=255)                                                             agcaAACTGG‑AAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGaatgataatgat                                                                                          >  2:1049779/1‑150 (MQ=255)                                                                aaaCGGG‑AAAAGTATCAAAAACAGGGTTGAAGGAAACAAATTCTGCAAGAAAAGCCGTTGGGAAATTAAACGGG‑AATCCCA‑TTGTAAAAGTTTATTC‑CGCCGGGGAATAAAGCCCTTT‑TATCCGCCCCAACAGAGaatgataatg                                                                                            <  1:48529/145‑1 (MQ=255)                                                                 aaCGGG‑AAAGTTATCGAAAACAGCGTAGCAGGAAACACAATCTGCAATAATACAGGCTGGGAAATTAACCGGG‑CATCCGC‑ATGTAAAATTTTATTA‑CGCGGGGGATTAAAGCCTGTG‑TAACCGCCCCAACAGAGAATGATAATGataata                                                                                      <  1:635626/150‑1 (MQ=255)                                                                        aaaaGTATCAGAAACAGCTTTGCAGGAAACACAACCGCCCAAAATATACGCTGGGAAATTAAACGGG‑TATCCCA‑ATGTAAAACGTTGTTA‑CGCGGGGGAGAAAAGACCCTG‑AAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAg                                                                               <  1:288338/151‑1 (MQ=255)                                                                        aaaaGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCGGGGAAATTAACCGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCGGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTa                                                                                <  2:908479/150‑1 (MQ=255)                                                                        aaaaGTATAAGAAACAGCTTTGCAGAAAACACAATCTGCCAGAATATACGCTGGGAAATAAAACGGG‑TATCACC‑ATGTGAAACGTTTTTA‑CGCGGGGGAGAAAAGCCCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAAt                                                                                  <  2:259555/148‑1 (MQ=255)                                                                          aaGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGcc                                                                             >  2:1124143/1‑151 (MQ=255)                                                                          aaGTATCAGAAACAGCGTGGAAGAAAACCCAATCGGCAAGAATATCCGCTGGAAAATAAAACGTG‑TACCCCC‑ATGTAAAGCGTTGTTC‑CGCGGGGGAGAAAAGCCCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGcc                                                                             <  2:222669/151‑1 (MQ=255)                                                                                cAGAAACAGCGATCAAGGAACCATATTCTTCAAAAATATACGCTGGAAAATAAAACGGG‑ATCCACC‑TTGTGAAACGTTATGC‑CGCTGGGGAGAAAAGCCCCTT‑GGCCCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGtt                                                                          <  1:512492/148‑1 (MQ=255)                                                                                 aaaaaCACCGTGGCAGAAACCCCACTCCGCCAGAATATACGCGGGGAAATTAAACGGG‑CAACACC‑ACGTAAAACGTTTTTC‑CCCCGGGGAGAAAAGCCCTTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATg                                                                      <  2:243161/149‑1 (MQ=255)                                                                                   aaccAGCGTGGAAGAAAACACATTATGCCAGATTATGCGGTGGAAAATAAACGGGG‑TACCACC‑TTGTAAAACGTTTTAA‑CGCGGGTGAG‑AAAAACCCTGAGAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGAt                                                                    <  1:1034998/148‑1 (MQ=255)                                                                                   aaaCAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGa                                                                     >  1:376754/1‑150 (MQ=255)                                                                                    aaCAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGAt                                                                    >  1:481367/1‑150 (MQ=255)                                                                                        gggtggCAGGAAACAAATTTTGCAAGAATATACGCGGGGAAATAAAACGGT‑TATCCCC‑ATGGAAAACGTTGTTA‑GGCGGGGGAGAAAAGACCTGG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCg                                                                <  1:775766/145‑1 (MQ=255)                                                                                        gCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGt                                                               >  2:308615/1‑151 (MQ=255)                                                                                           ttGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCCCC‑ATGTAAAACGTTGTTA‑CGCCGGGGATAAAAGACCCTT‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGtt                                                              <  1:769465/149‑1 (MQ=255)                                                                                           ttGCAGGAAACACAACCGGCCAAAATATACGCGGGAAAATAAAACGGG‑TATCCCC‑AGGTAAAACGTTTTTA‑GGCGGGGGATAAAAGCCCTTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGtt                                                              <  1:1045823/149‑1 (MQ=255)                                                                                             gCAGGAAACACAATCTGCCAGAATATACCCTGGGAAATTAAACGGG‑GATCACC‑ATGGAAAACGTTGTTA‑CGCCGGGGAGTAAAGGCCCTG‑TAACCGCCCCAATAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGt                                                          >  1:373992/1‑151 (MQ=255)                                                                                                      acaATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAccc                                                  >  2:548414/1‑150 (MQ=255)                                                                                                      acaAGCTGCCAGAATAGCCGCTGGGAAATAAAACGGG‑TATCACC‑ATGTAAAACGTGGTTA‑CGCTGGGGAGAAAAAACCTTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAAGTAGCCGTTAATGATAGCGTTTCGTAAAAACCCa                                                 <  1:439866/151‑1 (MQ=255)                                                                                                       caATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAAccc                                                  >  2:14043/1‑149 (MQ=255)                                                                                                        aaccgcccAGAATCTCCGCGAGAAAAATAAACGGT‑AATCACC‑ATGTAAAATGGTATTA‑CAACGAGAAATAAAGTCCCTA‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAg                                                <  1:611074/144‑1 (MQ=255)                                                                                                            tgcaataaTATAGGCT‑GGAAAATAAACGGG‑GAACCCCAATGCGAAACGTTGTTA‑CGCGGGGGAGAAAAATCCCTC‑TACCCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAATCCCAGctct                                            <  2:349408/144‑1 (MQ=255)                                                                                                                aaaaTATACGCGGGAAAATTAAACGGT‑TCCCCCC‑TTGTAAAAC‑TTGTTACCGCCGGGGAGAAAAACCCCTT‑AAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTg                                           <  1:982819/145‑1 (MQ=255)                                                                                                                     ataCGCGGGAAAATTAACCGGG‑TATCCCC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTAAAGCCCCTT‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTg                                   <  2:990852/150‑1 (MQ=255)                                                                                                                      taCGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTAAAGACCTTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAATCCAGCTCTGCATATTTGCg                                 >  1:785606/1‑151 (MQ=255)                                                                                                                      taCGCGGGGAAATAAAGCGGT‑TATCACC‑ATGTAAAACGTTGTTA‑CGCCGGAAAATAAAAACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGc                                  <  1:501045/150‑1 (MQ=255)                                                                                                                             ggAAAATTAACGGTTTATCCCC‑TTGTAAAATTTTTTTG‑CGCTGGGGTGAAAAGACCCTG‑AAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACaa                            <  2:982390/150‑1 (MQ=255)                                                                                                                                aaTTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATc                       >  1:461409/1‑151 (MQ=255)                                                                                                                                    aaaCGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCt                    >  2:1007386/1‑150 (MQ=255)                                                                                                                                        ggg‑TATCACC‑ATGTAAAACTTTGTTA‑CGCGGGGGAGAAAAGCCCCTG‑TAACCGCCCCAACAGAGAATGATAACGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGc               <  2:102723/151‑1 (MQ=255)                                                                                                                                         gg‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGCGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg              >  1:222320/1‑151 (MQ=255)                                                                                                                                          gttaccccc‑ATGTAAAACTTTGTTA‑CGCCGGGGAGAAAAGCCCCTG‑CAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCg              <  2:913865/144‑1 (MQ=255)                                                                                                                                              tcccc‑ATTTAAAACTTTGTTA‑CCCCGGGGAGTAAAGCCCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGAt            <  2:1066119/146‑1 (MQ=255)                                                                                                                                                 cc‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTg       >  2:903163/1‑151 (MQ=255)                                                                                                                                                       tgAAACGTTGTTA‑CGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAAt  >  2:762973/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                        GAGGTTAAAAGGCATGACACCGACCAGATTTTCTCTGCTCGCCAGTAATCCTACAAATGCAGCAAACTGG‑AAAAGTATCAGAAACAGCGTTGCAGGAAACACAATCTGCCAGAATATACGCTGGGAAATTAAACGGG‑TATCACC‑ATGTGAAACGTTGTTA‑CGCCGGGGAGTAAAGACCCTG‑TAACCGCCCCAACAGAGAATGATAATGATAATAAAATTAGCCGTTAATGATAGCGTTTCGTAAAAACCCAGCTCTGCATATTTGCGGACAAAAATCCCTAAGGCGATCCCTGGCAAT  >  NC_000913/2628450‑2628747

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