Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,379,638	0	G	A	38.1%	53.0 / 14.8	42	F45F (TTC→TTT)	frdC	fumarate reductase membrane protein FrdC
Reads supporting (aligned to +/- strand):  ref base G (22/4);  new base A (0/16);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 2.91e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGATTACGATGGTGGCAACCACAGTTACCGCCCAGAGACTTTT‑GATAATTGGCTCTGGTCCCATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGG‑TGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAG  >  NC_000913/4379400‑4379876                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  ggATTACGATGGGGGCAACCACAGTTACCGCCCAGAGACTTTT‑GATAATTGGCTCTGGTCCCATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGG‑TGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACAGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTc                                                                                                                                                                                                                                        >  2:269582/1‑250 (MQ=255)                cAACCACAGTTACCGCCCAGAGACTTTT‑GATAATTGGCTCTGGTCCCATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGG‑TGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACc‑acaca                                                                                                                                                                                                                      <  2:45701/251‑1 (MQ=255)                             cGCCCAGA‑ACTTTTTGATAATTGCCTCTGGTCCCTTTTTTTCGTGTTTTAAAATGAGATTGGCCGTTTGCGGGGCCGGTTAAAACCAGGTTTTGG‑TGGGCAGCAGGGCTGCCGCAAGAGTGATCAGGTAAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAACCCCGCCCAGGC‑TTCGGGGCCATTTTTC‑AGGG‑AAAA‑CAGCCCGAAAATCAGTTAAAT‑GCTGAACC‑ACACAGCCGGAA‑ccgccg                                                                                                                                                                                                        <  2:30124/251‑1 (MQ=255)                              gCCCAGAGACTTTT‑GATAATTGGCTCTGGTCCCATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGG‑TGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACC‑ACACAGCCGGAA‑ccgccgg                                                                                                                                                                                                       >  2:183412/1‑250 (MQ=255)                                                aaTTGGCTCTGGTCCCATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGG‑TGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGt                                                                                                                                                                                      >  2:461351/1‑251 (MQ=255)                                                 aTTGGCTCTGGTCCCATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGG‑TGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGt                                                                                                                                                                                      >  1:375758/1‑250 (MQ=255)                                                             cccATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGG‑TGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATaaa                                                                                                                                                                          >  1:152167/1‑250 (MQ=255)                                                              ccATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGG‑TGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAAc                                                                                                                                                                        >  1:508325/1‑251 (MQ=255)                                                                    tttCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGG‑TGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGTC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAACGCGATACAACGGCAAt                                                                                                                                                                  >  2:18530/1‑251 (MQ=255)                                                                              taaaaTGATATGGGCCGCTTTCGGTGCCAGTTAAACCCAGTTTTTGG‑TGTGCAGGAGAGCGGCCGCGAGGGTGATCAGTTAAATGACCACGATAACCGGGTTTT‑AAAAAAAGTCGCCGAACCCCGCCCGGGC‑TTCCGGGCCTTTTTAC‑GGGG‑CAAA‑CAGCCCAAAAATCAGTTAAAT‑GC‑GAACCCACCCAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTccac                                                                                                                                                         <  1:18530/247‑1 (MQ=255)                                                                               aaaaTGATTTTGGCCGCTTTGGGTGCCAGTACAAACCAGGTTTGGG‑TGACCAGCAGAGCGGCCCCCAGAGTGATCATTTTAAGAATCGCAATACCCGGTTTTT‑TTAAAAAGTCGCCAAACCCCGCCCAGCC‑TTCCGGCCTTTTTTTC‑AGGG‑CAAC‑CCCCCCAAAAATCATTTCAAT‑GCTGACCC‑CCACACCCGAAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTccacca                                                                                                                                                       <  1:140577/249‑1 (MQ=255)                                                                                  aTGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGG‑TGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAggg                                                                                                                                                    >  2:63572/1‑250 (MQ=255)                                                                                     atatTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGG‑TGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTgg                                                                                                                                                  >  1:438690/1‑250 (MQ=255)                                                                                       tagggcccgtttcgggcccaattCAAACCAGGTTTTGGGTGTGCAGCAGGGCGGCCGCAAGAGTGATCAGGTAAATGATCACGATA‑CCGGGTTTTTGAAAAAAGTCGACGAACCCCCCCCAGGC‑TCCCGGGCATTTTTTA‑AGGG‑AAAC‑CAGCCCAAAAATAAGTAAAAT‑GCTGAACC‑ACACAGCGGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGa                                                                                                                                                 <  2:152167/229‑1 (MQ=255)                                                                                                 ttgggtccctgttcaaacagtgTTTTGG‑TGTGCATAAGACCTGCCGCCAGATTCATCAGGTTAATAATCACGAGATCCGGTTTTG‑GTAAAAAGTCAACGAACCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGC‑CAAC‑CACCCAAAAAATCAGTTCAAT‑GCGAAACC‑ACAAAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGcc                                                                                                                                      <  2:385314/229‑1 (MQ=255)                                                                                                             caaaaCCGGGTTTGGG‑TGTCCAGCAGAGCGGCCGCCAGAGTGATCAGGTTAATGACCACAATACCCGGGTTTT‑GAAAAAAGTCGACAAACCCCGCCCAGGC‑TTCCGGCCCTTTTTTC‑AGGG‑AAAA‑CAGCCCAAAAATCATTTCAAT‑GC‑GAACCCACACAGCCGAAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGttt                                                                                                                         <  1:286051/249‑1 (MQ=255)                                                                                                                                      agagCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTc                                                                                                 >  2:63071/1‑251 (MQ=255)                                                                                                                                       gagCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAACAAGTCGACGAATCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTcc                                                                                                >  1:406875/1‑251 (MQ=255)                                                                                                                                          cTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGATTACGTTTAGTCGTCATGTTGCACTCCtt                                                                                              >  2:408625/1‑250 (MQ=255)                                                                                                                                           tGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCtt                                                                                              >  1:503174/1‑249 (MQ=255)                                                                                                                                                caagggtgtacagttaaaTGATCACAATA‑CCGGGTTTTTGAAAAAAGCCGACGACCCCCCCCCAGGT‑TGCCGGGCCATTTTT‑‑AGGG‑CAAAACAGCCCAAAAATCAGTACAAT‑GCTAAACC‑CCACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCATCATGTTGCACTCCTTAGCGTgg                                                                                       <  2:120314/235‑1 (MQ=255)                                                                                                                                                      ttagtaGGTTAATAATCCCGATACCCGGTTTTT‑GAAAAAAGCCGACGAACCCCCCCCGGCC‑TTCGGGGCTATTTTCC‑AGGG‑CAAA‑CAGCCCAAAAATCATTTCAAT‑GTGGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCa                                                                                  <  1:376026/245‑1 (MQ=255)                                                                                                                                                      aaagcgggTTAACGTTCACGATA‑CCGGGTTTTTGAAAAAAGTGGACAAACCCCGCCCGGGC‑TTCGGGACTATTTTAC‑GGGT‑CAAA‑CCCCCCGAAAAACA‑TGAAATTCCTGACCC‑ACACAGCAGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAg                                                                                 <  1:183412/245‑1 (MQ=255)                                                                                                                                                       aaTCAGTTTATTGATCCCGAAAACGGGGTTTT‑GAAAAAAGTCGACGAACCCCGCCCAGTC‑TTCGGGGCCTTTTTTC‑AGGG‑CAAC‑CGGCCCAAAAAGCATTTCAAT‑GC‑GAACCCACACAGCGGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAg                                                                                 <  1:195656/249‑1 (MQ=255)                                                                                                                                                                 aatgaccaagataccgggtttttaaaaaaaGGCGACGACCCCCGCCCAGGC‑TCCCGGGCCTTTTTTC‑GGGG‑CAAA‑CACCCAAAAAACCATTTAAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATa                                                                       <  1:520005/226‑1 (MQ=255)                                                                                                                                                                       acaatacccGGGTTTT‑GTAAAAAGTCGACGA‑TCCCCCCCCGGCCTCCCGGGCC‑TTTTTC‑AGGGGCAAA‑CAGCCCAAAAATCAGTCCAAT‑GC‑GAACCCACACAGCCGAAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGt                                                                <  1:258672/244‑1 (MQ=255)                                                                                                                                                                        agataccggggTTTT‑GAAAAAAGTCGACGACCCCCGCCCAGGC‑TCCCGGGCCATTTTTC‑AGGG‑CAAA‑CACCCCAAAAATCAGTTCAAT‑GCTAACCC‑ACACGGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTc                                                               <  1:26239/243‑1 (MQ=255)                                                                                                                                                                         gATAACCGGGTTTT‑TTAAAAAGTCGACGAATCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCAGAAAATCAGTTCAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGCGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTGTCCACCAGGTGGACGTCATTGGCCGTACATACGGGTTACGTTTAGACGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGACGCGATAAGAAAgact                                                              >  2:270514/1‑248 (MQ=255)                                                                                                                                                                             aCCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCg                                                          >  2:42700/1‑251 (MQ=255)                                                                                                                                                                                   tttt‑GTAAAAAGTCGACGAATCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACttt                                                    >  2:515614/1‑251 (MQ=255)                                                                                                                                                                                          caaaaGTCGACGAATCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACC‑CCACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGAAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGCCGGCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGCCTTTCGAACTTTCTCCtt                                              >  2:399997/2‑251 (MQ=255)                                                                                                                                                                                             aagccgacaaaccccGCCCAGGC‑T‑CCGGGGCATTTTTCCGGGG‑CAAA‑CAGCCCAAAAAGAAGTTAAAT‑GCGAAACC‑CCACA‑CCGGAAACCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGc                                           <  1:27094/239‑1 (MQ=255)                                                                                                                                                                                             aaGTCGACGAATCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGc                                           >  1:98705/1‑251 (MQ=255)                                                                                                                                                                                                ccgacgacccccGCCCAGGC‑TCCCGGCCCTTTTTAC‑AGGG‑CAAA‑CACCCCAAAAATCATTTCAAT‑GCTAACCC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCct                                        <  1:63071/242‑1 (MQ=255)                                                                                                                                                                                                   ccgactCCCCCCCAGGC‑TTCCGGCCCTTTTTTC‑AGGG‑CAAC‑CAGCCAAAAAATAAGTTAAAG‑GCTAAACC‑ACACAGCCGGAC‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCctgct                                     <  2:508325/246‑1 (MQ=255)                                                                                                                                                                                                     gaATCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCctgctg                                    >  1:117797/1‑250 (MQ=255)                                                                                                                                                                                                     gaATCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCctgctg                                    >  1:262451/1‑250 (MQ=255)                                                                                                                                                                                                                       ccGGGCAATTTTCC‑AGGG‑CAAA‑CAGCCCAAAAATCATTTCAAT‑GCTGACCC‑ACCCAGCGGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGAt                   <  1:452269/250‑1 (MQ=255)                                                                                                                                                                                                                          ggCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAAtttttt                                                                                                                                                             >  1:376818/1‑109 (MQ=255)                                                                                                                                                                                                                           gcctttttcC‑AGGG‑CAAA‑CAGCCCAAAAAGCAGTTAAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACg              <  1:46852/247‑1 (MQ=255)                                                                                                                                                                                                                                        gg‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGCCATTGGCCGTACATACGGTTTACGTTTAGTAGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGACAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGGGTTTCggga    >  2:382676/1‑248 (MQ=255)                                                                                                                                                                                                                                        gg‑CAAA‑CAGCCAGAAAATCATTTCAAT‑GCTGAACC‑ACACAGCGGAAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAg  <  1:500915/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                  GGATTACGATGGTGGCAACCACAGTTACCGCCCAGAGACTTTT‑GATAATTGGCTCTGGTCCCATTTTTTCGTCTTTTACAATGATATTGGCCGCTTTCGGTGCCAGTTCAAACCAGGTTTTGG‑TGTGCAGCAGAGCTGCCGCCAGAGTGATCAGGTTAATGATCACGATAACCGGGTTTT‑GTAAAAAGTCGACGAATCCCGCCCAGGC‑TTCCGGGCCATTTTTC‑AGGG‑CAAA‑CAGCCCGAAAATCAGTTCAAT‑GCTGAACC‑ACACAGCCGGAA‑CCGCCGTGCCTTCGCGCAGCATGTAAAAGCGATAAAACGGCAATTTTTTCCACCAGGTGGACGTCATTGGCCGTACATACGGTTTACGTTTAGTCGTCATGTTGCACTCCTTAGCGTGGTTTCAGGGTCGCGATAAGAAAGTCTTTCGAACTTTCTACTTTGCCCTGCTGAATGGCCGCAGCCGGATCGACGTGTTTCGGGCAG  >  NC_000913/4379400‑4379876

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