Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,639,764	G→A	S181S (TCC→TCT)	arcA ←	DNA‑binding transcriptional dual regulator ArcA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,639,764	0	G	A	92.0%	127.5 / ‑1.5	53	S181S (TCC→TCT)	arcA	DNA‑binding transcriptional dual regulator ArcA
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base A (21/25);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

GAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGTAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGT  >  NC_000913/4639524‑4640003                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   gAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAAttt                                                                                                                                                                                                                                          <  1:309043/251‑1 (MQ=255)        ggTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑gg                                                                                                                                                                                                                                  <  1:644765/251‑1 (MQ=255)            gAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑GCATTTTCTTCAGCAGTTCAGCCCGCGGCTGAATTTTTACT‑GGGttt                                                                                                                                                                                                                              >  2:184132/1‑251 (MQ=255)            aaaaaaTAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGttt                                                                                                                                                                                                                              <  2:144326/250‑1 (MQ=255)                 aTAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGGGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTcac                                                                                                                                                                                                                          >  1:554439/1‑250 (MQ=255)                    aaaaCGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGGCTACAGGACGGGCGTGGGGTTTCATCCCACGGCCGGATC‑TTTTCTTCCGCCGTTCAGCCCGGGACTGAATTTTGGCT‑TGGTTTTcaccga                                                                                                                                                                                                                       >  2:134791/1‑247 (MQ=255)                               aaaaaGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATgg                                                                                                                                                                                                           <  1:507914/251‑1 (MQ=255)                                  aaGCGCCGTTTTTTTTGCCGGTGGTAAAGCCGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCTAGGGTTTT‑ACAGAAGTGAAGCATGgcgc                                                                                                                                                                                                        >  1:537919/1‑251 (MQ=255)                                     cgccggtgttttttACGGTGTTACAGCCGTTTAATTTGCCAGGTCACCACAGCATCGTTATCCCTCACCTTGGATGGTGGCGATGATTTTCGGCGTATCCGGCGTAGATTCAAAATTTTTCCTAATCCGGCGGATCGTCCCGTCTTCAGTACTGTCCTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GCGTTTTCACTGAAGTGAAGCATGGCGCGGa                                                                                                                                                                                                     <  2:176875/243‑1 (MQ=255)                                                     ggtggtAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGGTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCT‑GGGTTTTCACAGAAGGGAAGCATGGCGCGGAACTAGCTGggcggc                                                                                                                                                                                       >  2:524625/1‑249 (MQ=255)                                                                cGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTATATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTCTCACAGAAGTGAAGCATg                                                                                                                                                                                                            <  2:97576/217‑1 (MQ=255)                                                                cGATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATg                                                                                                                                                                                                            >  1:97576/1‑217 (MQ=255)                                                                 gATTAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTc                                                                                                                                                                         >  1:331189/1‑251 (MQ=255)                                                                    tAATCTTCCAGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGGCCCGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTTAGCAGTTCAGGACGAGACTGAATTTTTCCT‑GGGTTTTTACCGGAGGTAAGCATGTCGCGGGACACACTGCGCGGcggtttttacgcctcgcc                                                                                                                                                                      >  2:418005/1‑234 (MQ=255)                                                                             aGATCACCACAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGcc                                                                                                                                                            >  1:641416/1‑251 (MQ=255)                                                                                     acAGAAGCGATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCaacg                                                                                                                                                    <  1:484887/250‑1 (MQ=255)                                                                                             gATAACCTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGACACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCCAGCGACTGAATTTTGCCT‑GGGTTTTCCCAGAAATGAAGCATGGCGCGGGACTCGCTGCGCGGCCGCTTGTACTGCTCGCCCT‑CAGGGCCGGTCAACGCACGGCTg                                                                                                                                            >  1:375605/1‑251 (MQ=255)                                                                                                   cTTCACCGTGGATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATg                                                                                                                                      >  1:315858/1‑251 (MQ=255)                                                                                                           gggaTGGTGGCGATGCTTTCCGGCGTATCCGGCGTAGATTCGAAACGTTTACGAATACGGCGGATCGTCACGTCTTCAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCATCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATTGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTc                                                                                                                              <  1:560877/248‑1 (MQ=255)                                                                                                            ggatgtgggCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCTCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTcc                                                                                                                             <  2:107832/244‑1 (MQ=255)                                                                                                                 gtggCGATGATTGCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCTG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GTGTTTTCACAGAAGTGAAGCATGGCGCGGAACGCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAAcc                                                                                                                        <  1:609585/251‑1 (MQ=255)                                                                                                                   ggCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAATTCAGCACGAGACTGAATTTTGCCT‑GGGGTTTCACAGACATGAAGCATGGCCAGGAACTCGCTGCGCGGCCGCGTGTACGGATCGCCATTCAGT‑CCGATCAACGTACGGCCGTTGATGTACAGTTCCCCACCAt                                                                                                                      >  2:315418/1‑251 (MQ=255)                                                                                                                    gCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCAtt                                                                                                                     >  1:558460/1‑251 (MQ=255)                                                                                                                      gatgatTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTGCGGTCGTGTGGCTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGGTCAGCACGGGAATGAATTTTGCCT‑GGGTTTTCACATAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTCCTGATCGCCAT‑CCGGGCCGATCAACGAACGGCTGATTAAGTCCAATTCCCCACCATTg                                                                                                                    >  2:622466/1‑250 (MQ=255)                                                                                                                             tccggcgtatccggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGt                                                                                                             <  1:524625/250‑1 (MQ=255)                                                                                                                                ggcgtatccggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCt                                                                                                        >  1:603927/1‑251 (MQ=255)                                                                                                                                 gcgtatccggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTTAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCGCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCtt                                                                                                       <  1:120998/251‑1 (MQ=255)                                                                                                                                  cggatccggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGTGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCttt                                                                                                      <  2:641416/248‑1 (MQ=255)                                                                                                                                        cggcgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGc                                                                                                >  1:416708/1‑251 (MQ=255)                                                                                                                                           cgtaGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTT‑TATGCATTCAACGCTac                                                                                             >  1:540172/1‑251 (MQ=255)                                                                                                                                                   cGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGGA‑GCTTTAAACGCTACGACGttct                                                                                     >  1:399669/1‑251 (MQ=255)                                                                                                                                                    gAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGttctt                                                                                    <  1:184132/251‑1 (MQ=255)                                                                                                                                                     aaaTGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGttcttc                                                                                   <  1:241149/251‑1 (MQ=255)                                                                                                                                                             aCGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACa                                                                            <  2:280951/250‑1 (MQ=255)                                                                                                                                                                    cgtcggttcttatcGTCTACAGGTCGGGCGTGTGGTTTTAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCTGCTTTTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGTCGTTCTTCGCTGACAGTACCCAg                                                                    <  2:358760/238‑1 (MQ=255)                                                                                                                                                                    cgggggATCGTCACGTCTACAGTACGGGCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAg                                                                    <  2:244225/247‑1 (MQ=255)                                                                                                                                                                    cggcggATCGTAACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGATCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGGTGTCCAGATGCCAACCATTGAGCTTGTC‑GCTTTCAACGCTACGACTTTTTTCGCTGCCAGTAccc                                                                      >  2:195702/1‑249 (MQ=255)                                                                                                                                                                        ggATCGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTc                                                                >  1:634926/1‑251 (MQ=255)                                                                                                                                                                            cGTCACGTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATg                                                             <  1:612807/250‑1 (MQ=255)                                                                                                                                                                                  gTCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGACCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGATGACAGTACCCCGATTCATGGTAAggg                                                      >  2:569550/1‑251 (MQ=255)                                                                                                                                                                                   tCTACAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGCAAGTACCCAGATTCATGGTACgggg                                                     >  2:640817/1‑250 (MQ=255)                                                                                                                                                                                      aCAGTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAg                                                  <  1:158408/251‑1 (MQ=255)                                                                                                                                                                                         gTACGGTCGTGTGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAgc                                                 >  1:510294/1‑248 (MQ=255)                                                                                                                                                                                                ggtgtggTTTCAGCTCACGGCCGG‑TCAGTTGCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGGTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGt                                        <  2:634926/246‑1 (MQ=255)                                                                                                                                                                                                           aGCTCTCGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CATGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCTTGCACGAATCg                              <  2:537919/250‑1 (MQ=255)                                                                                                                                                                                                                 cGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGCGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACGGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCAAGAATAGTCAGTTc                       >  2:340215/1‑251 (MQ=255)                                                                                                                                                                                                                 cGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTc                       >  1:647191/1‑251 (MQ=255)                                                                                                                                                                                                                   gCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCAc                     <  1:499732/251‑1 (MQ=255)                                                                                                                                                                                                                   gCCGG‑TCATTTTCTTCAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGACAGCTTGTACTGCTCGACAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGATGCGTGCACGACTCGTCAGTTCAc                     >  1:276839/1‑251 (MQ=255)                                                                                                                                                                                                                                  tcAGCAGTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACg         <  1:239023/249‑1 (MQ=255)                                                                                                                                                                                                                                    agcagTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCATCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGGCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTTAACGGttt     <  2:127403/251‑1 (MQ=255)                                                                                                                                                                                                                                     gcagTTCAGCACGAGACTTAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGCAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTg    <  2:15708/251‑1 (MQ=255)                                                                                                                                                                                                                                       agTTCAGCACGAGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCCGATTCATGGTACGGGACAGCAGGTTGCGGGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTggg  >  2:532891/1‑250 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   GAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTCTACAGTACGGTCGTGCGGTTTCAGCTCACGGCCGG‑TCATTTTCTTCAGCAGTTCAGCACGGGACTGAATTTTGCCT‑GGGTTTTCACAGAAGTGAAGCATGGCGCGGAACTCGCTGCGCGGCAGCTTGTACTGCTCGCCAT‑CAGGGCCGATCAACGAACGGCTGTTGATGTCCAGTTCCCAACCATTGAACTTGTA‑GCTTTCAACGCTACGACGTTCTTCGCTGACAGTACCCAGATTCATGGTACGGGACAGTAGGTTGCGTGCACGAATCGTCAGTTCACGCGGGTTGAACGGTTTGGT  >  NC_000913/4639524‑4640003

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