Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,862,661	T→G	26.5%	pseudogene (943/1617 nt)	glvC ←	pseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,862,661	0	T	G	26.5%	38.0 / 5.4	40	pseudogene (943/1617 nt)	glvC	pseudogene, arbutin specific enzyme IIC component of PTS;enzyme; Transport of small molecules: Carbohydrates, organic acids, alcohols; PTS system, arbutin‑like IIB component; PTS system, arbutin‑like IIC component
Reads supporting (aligned to +/- strand):  ref base T (9/19);  new base G (9/2);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 1.07e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

CCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑AATCA‑GC‑AAG‑CCCG‑CCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAA  >  NC_000913/3862425‑3862889                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              cccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTT‑g                                                                                                                                                                                                                                            <  1:262285/239‑1 (MQ=255) cccAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCA‑GTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAC‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTAAGGTT‑Gc                                                                                                                                                                                                                                           >  1:885447/1‑239 (MQ=255)           aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTTCCCCCAGCATGGC‑GGTTAAGGGT‑GCGGG‑GAtgat                                                                                                                                                                                                                                 >  1:17323/1‑236 (MQ=255)           aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGT‑GCGGG‑GAtac                                                                                                                                                                                                                                  >  1:1025217/1‑236 (MQ=255)           aTTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGA‑ACATTACGGTCGACATTGAGGCCGCCCGCAC‑GGCGTGTACCGCAAACAGCAACCGTGAAATGAACAGGAAGGTAAA‑TTCCACCGGTTCGGTCATTCCCCCCAGCATTGC‑GGTTAAGGTT‑Gctgtgattcg                                                                                                                                                                                                                                  >  2:1112886/1‑229 (MQ=255)                           tatCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGAT‑GCGGG‑GATTA‑GCAAAG‑CCCC‑CCAt‑ttt                                                                                                                                                                                                             >  1:202776/1‑239 (MQ=255)                                gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCCAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTTAGGGT‑GGGGG‑AATCACGC‑AA‑‑CCCG‑CGCC‑TTTTAccc                                                                                                                                                                                                        >  1:954372/1‑238 (MQ=255)                                gcgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGTT‑GGGGG‑GAACA‑GC‑CAG‑CCCG‑CCCC‑TTTTAccc                                                                                                                                                                                                        >  1:1138245/1‑238 (MQ=255)                                 cgTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGT‑GCGGG‑GATTA‑GA‑GAA‑CCCC‑CCCC‑TTTTTcccca                                                                                                                                                                                                      >  1:436527/1‑237 (MQ=255)                                                    gATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑GAACA‑GC‑AAC‑ACCGCCCCC‑TTTTACCC‑ATTTTTCGGTGCgggaggga                                                                                                                                                                                   >  1:797896/1‑239 (MQ=255)                                                      tCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑GATCA‑GC‑AAA‑CCCG‑CCCC‑TTTTTCCCCCTTTTCCCGGGCGGCAGGGaaa                                                                                                                                                                                 >  1:750574/1‑238 (MQ=255)                                                        cAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGATAAGGTT‑GCGGG‑GAACC‑GC‑AAA‑CCCC‑CCCC‑CTTTAACCGCGTTTTCTGTGCggtcgtggagga                                                                                                                                                                               >  2:140429/1‑229 (MQ=255)                                                         attttcGCGGAAAAACATAGCAAAACACCCCCCCGCCCGTGTCCCCCACCCCACAAAAGAGAACCATTACGGCCGACATTGAGGCCGCCAGCAC‑GGCGTGTCCCGCAAACAGCAGAGGTAAAATGAACAGGAAGGTAAA‑TTCAAGCGGTTCGGTAATCCCCCCCAGCATGGC‑GGTTAGGGTT‑GCGGA‑AATAA‑GC‑AAG‑CCCG‑CCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTAc                                                                                                                                                                              <  2:291425/237‑1 (MQ=255)                                                           ttttGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑GATCA‑GC‑AAA‑CCCG‑CCCC‑TTTTTCCCCATTTTTCGGTGCGGCAGGGAAGTACa                                                                                                                                                                             >  1:998706/1‑238 (MQ=255)                                                                                          cGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTCAAGGG‑GCGGG‑GAGCA‑GC‑ACA‑CCCC‑CCAC‑TTTTACCCCTATTTTCGCGGCGGCACGGACGGACATCGCTAAAGAACTACCCA‑CGGGCCCAAAGagt                                                                                                                                            >  2:498663/1‑237 (MQ=255)                                                                                           cccagttttGCCCACCCCCCCAAAGAAAAACTTTAGGGGCGAAGTTGA‑GCCCCCCGCGCGGGCGTGTACCGCAAACAACAAGGGTAAAAGGAACAGGAGTGAAAA‑TTCCGGTGGTTCGGAATTCCCCACGAGTATGGC‑GGTTAAGGTT‑CGGGG‑AATCA‑GC‑AAG‑CCCG‑CCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCA‑CGGCACCAAAGAt‑tt                                                                                                                                          <  2:599458/232‑1 (MQ=255)                                                                                              caggTTGCCCACCCCACAAAAAAAATACATTACGGCCAACATGAAGCCCCCCAGCAC‑GGCGTGTCCCGCAAACAGAACCGGGGAAAGGAACAGAAAGTTAAA‑TTCCGGCGGTTCGGTAATCCCCCCCAGCATGGC‑GGTTAAGTTT‑GCGGG‑AATCA‑GC‑AAG‑CCCG‑CCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTg                                                                                                                                       <  1:449534/236‑1 (MQ=255)                                                                                                        accacacaaAAGAGATACATAACGGTCGACATTGAGGCCCCCAGCAC‑GGGGTTTCCCGCAAACAGCAACGGTAAAAGGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATCCCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑AATCA‑GC‑AAG‑CCCG‑CCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTg                                                                                                                             <  1:335246/239‑1 (MQ=255)                                                                                                          cacaccaAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGGTAAAGGT‑GCGGG‑GAACA‑GC‑AAGACCCC‑CCCT‑TTTTACCCGATTTT‑CGGTGGGGCAGTGAAGGACAACGCTAAAAAAAAGCCCA‑CCGGCCCCAAGAA‑TTTTTAGtttaccgggc                                                                                                                           >  1:1182898/1‑231 (MQ=255)                                                                                                                      aTACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGT‑GCGGG‑AATAA‑GC‑AAA‑CCCC‑CCACTTTTTACCCGTTTTTCGGGGCGGGCGGTG‑AGAACAAACCTAAAAAAAAGCCCC‑CAGGCCCAAAGAT‑TTTTTAGTGTGCGGGCCAGGCAAAACCg                                                                                                                >  2:1110160/1‑238 (MQ=255)                                                                                                                          aTTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑AATCA‑Gc                                                                                                                                                                                                                               <  1:926555/129‑1 (MQ=255)                                                                                                                                gTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTTAAGGTG‑CGGGG‑AACAA‑CC‑AAG‑CCC‑‑CCCC‑CTTTACCCCCTTTTCCGGTGGGGCAAGGAAGAACACCCCTAAAAAAAAACCCC‑CGGCACCAAAGATGTTTTGAGTTACC‑TGGCGGGCAAAACCCCCCTCCcggca                                                                                                     >  2:644905/1‑237 (MQ=255)                                                                                                                                 tCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGGT‑GCGGG‑AATCG‑AC‑AAG‑CCCC‑CCCC‑TTTTTCCCGTTTTTCTGGGGCGGCCGAGAAAGACAAAACTTAAAAAATTGCCA‑CAGCACCCAAGAT‑TTTTAAGGTACCGTTCTCTGCATAAACGCCCTCCggggaa                                                                                                    >  2:954132/1‑236 (MQ=255)                                                                                                                                        tGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGCGGGTAAAGGTT‑CCGGG‑AA‑CA‑AC‑AAGCCCCC‑CCCC‑TTTTCCCCGTTTTTCGGGGGC‑GCAGGGGAGGACATACCTCAACAAATGCCCACCGGACCCAAAAAT‑TTTTTGTTTCCCGCGCGGGGAACAACCCCCTCCCGGGaacacgcat                                                                                              >  2:840190/1‑234 (MQ=255)                                                                                                                                                     ccgcgggtgtACCGAAAACAACAACTGTAGAATAAACAGGAATGTAAA‑TCCCAGC‑TTTCCGTAACTCCCCCCAAGGAGGG‑GGCTTAGGTT‑GGGGGCAACAA‑GC‑AAC‑CCCG‑CCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCgg                                                                                 <  2:1107366/233‑1 (MQ=255)                                                                                                                                                     cc‑GGCGTGTACCGAAAACAGCAAGGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTGGGTAATCCCCACCAGCATGGC‑GGTTAAGGTT‑GGGGG‑AATCA‑GC‑AAG‑CCCG‑CCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTc                                                                                   <  1:140429/235‑1 (MQ=255)                                                                                                                                                        ggCGTGTACCGAAAACAGCAACGGTGAAATAAACAGGAAGGTAAA‑TCCCAGCGGTTCGTTAATCCCCACCAGCAGGGC‑GGTAAGGTTT‑GGGGA‑AATCA‑GC‑AAG‑CCCG‑CCAC‑TTTTACCCGATTTTCCGGTg                                                                                                                                                                                            <  1:761077/131‑1 (MQ=255)                                                                                                                                                                   cAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTT‑GCggg                                                                                                                                                                                                                                        <  1:650812/82‑1 (MQ=255)                                                                                                                                                                   cAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAAGGT‑GCggg                                                                                                                                                                                                                                        >  2:650812/1‑82 (MQ=255)                                                                                                                                                                                 gAAATGAAGAGGAAGGTAAA‑TCCCGGCGTTTCGTAAATTCCCACCAGTATGCG‑GGTTAAGGTT‑GGGGA‑AATCA‑GA‑AAA‑CCCG‑CCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                      <  1:851352/238‑1 (MQ=255)                                                                                                                                                                                 aaaagaaaaaggaaggaAAA‑TTCCGGGGGGTGTGTAATCCCCCCCCGCGGGGG‑GGTTAGGGGG‑GCGGA‑AACCA‑GC‑AAG‑CCCG‑CCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                      <  2:797896/229‑1 (MQ=255)                                                                                                                                                                                  aaaGGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATCCCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑AATCA‑GC‑AAC‑CCCG‑CCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAg                                                           <  1:367797/232‑1 (MQ=255)                                                                                                                                                                                       aaCAGGAAGGAAAA‑TTCCGGGGTTTCGGTAATTCCCCCAAGCATGGC‑GTTAAAGGTG‑GCGGA‑ACTCA‑GA‑AAG‑CCCG‑CCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGcc                                                 <  2:1182898/237‑1 (MQ=255)                                                                                                                                                                                         cgaagaagTAAATTTCCGGCGGTTCGAAAACCCCACCCGGCAGGGG‑GGTTAAGTGT‑GGGGT‑AATCA‑AC‑GAG‑CCCG‑CCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAg                                              <  2:113379/232‑1 (MQ=255)                                                                                                                                                                                                       tcgatgggTACGGTAATCCCCCCCAGCAGGGC‑GGTTAAGGTG‑GCGGG‑AATCA‑GC‑AAG‑CCCG‑CCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATgcc                                <  1:498663/233‑1 (MQ=255)                                                                                                                                                                                                                      acccccACCGGCAGGGC‑GTTAAAGGTG‑GCGGA‑AACCA‑AC‑AGG‑CCCG‑CCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGc                 <  1:209973/236‑1 (MQ=255)                                                                                                                                                                                                                           caccaGCA‑GGG‑GGTTTAGGTTGGCGGG‑AACAA‑GC‑AAG‑CCCG‑CCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGAc            <  1:1115232/239‑1 (MQ=255)                                                                                                                                                                                                                                 cATGGC‑GGTTAAGGTG‑GCGGG‑AATCA‑GC‑ACG‑CCCG‑CCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGAtgct                                                      <  1:981100/191‑1 (MQ=255)                                                                                                                                                                                                                                       ggttaaaGGTG‑GCGGG‑AAGCA‑GC‑AGC‑CCCG‑CCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa  <  1:398086/233‑1 (MQ=255)                                                                                                                                                                                                                                       cgtttAGGGTG‑GGGGA‑AATCA‑GA‑ACG‑CCCG‑CCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATaaa  <  1:1124370/235‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                              CCCAATGGCGATTTGGGTCAGCATCATATCCGCGTGGTTGCTGAACATCGGGATCCAGTTTTGCGGTAAAACCTAGTCAATCAGACCTCCGCCCATGTTGCCCACCACACCAAAGAGATACATTACGGTCGACATTGAGGCCGCCAGCAC‑GGCGTGTACCGCAAACAGCAACGGTGAAATGAACAGGAAGGTAAA‑TTCCAGCGGTTCGGTAATTCCCACCAGCATGGC‑GGTTAAGGTT‑GCGGG‑AATCA‑GC‑AAG‑CCCG‑CCAC‑TTTTACCCGATTTTCCGGTGCGGCAGTGAAGTACATCGCTAAAGAAATGCCCA‑CGGCACCAAAGAT‑TTTTGAGTTACCGTGCAGGGCAAAACCGCCTTCCGGGAACAACGATTTCAGCGGCTCGGCGCTCAAACTGAACTCTTGCAGATGCTGCGCCCAGTACATCTGAATGCCGCCTTCAACAGCTGCCGGACCAAAGATAAA  >  NC_000913/3862425‑3862889

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