Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,482,648	0	A	C	28.6%	132.1 / 11.1	70	V397G (GTT→GGT)	valS	valine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (36/14);  new base C (0/20);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 8.60e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAA‑CGTTACCCGCTT‑CGTCATACCATGCCGGGATACGG‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTAT  >  NC_000913/4482402‑4482894                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       gaagaGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAA‑CGTTACCCGCTT‑CGTCATACCATGCCGGGATACGG‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGACTGTC‑GCCGTTCTCAACCGc                                                                                                                                                                                                                                                    >  2:186767/1‑251 (MQ=255)  acgagaACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCAAACATAAA‑CGTTACCCGCTT‑CGTCATACCATGCCGGGATACGG‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCTACCGCt                                                                                                                                                                                                                                                   >  2:16522/3‑251 (MQ=255)           aGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAA‑CGTTACCCGCTT‑CGTCATACCATGCCGGGATACGG‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAAccgcc                                                                                                                                                                                                                                          <  2:446872/251‑1 (MQ=255)                 cGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAA‑CGTTACCCGCTT‑CGTCATACCATGCCGGGATACGG‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGACGT‑ACATGTTTTCGTACTGCTTCGGTACGACC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTc                                                                                                                                                                                                                                    >  2:35659/1‑251 (MQ=255)                  gagaACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAA‑CGTTACCCGCTT‑CGTCATACCATGCCGGGATACGG‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCg                                                                                                                                                                                                                                   >  2:219174/1‑251 (MQ=255)                    aaccgtcttctacctgccacAGGACAACATCAGCACCGAGGTTATTTTCTTTCCGCACTTCGTCTTCGTTGCGGCCACCATAAA‑CGTTACCCGCTT‑CGTCATACCATGCCGGGATACGG‑TGCCCCCCCC‑ACAACTGACGA‑GAGATACCCCAGTCCTGCATACCGCGCATCCAGGAGAATT‑ACATGTTTTCGTACTGCTTCGGTACGACC‑GGAATGTC‑TCCGTTCTCACCCGCTTCAACCGCCGGTTTCGcc                                                                                                                                                                                                                                 <  2:425914/232‑1 (MQ=255)                     caCGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGCTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAA‑CGTTACCCGCTT‑CGTCATACCATGCCGGGATACGG‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑CCATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGcc                                                                                                                                                                                                                                 >  2:89283/2‑250 (MQ=255)                      aCGTCTTCGTCCTGACGCAGGACAACACCAGCCCCGAGTTTATTTTCTTTACGCACTTCTTCTTCTTTGCGGCCAACATAAA‑CGTTCCCCGCTT‑CGTCATACCATGCCGGGATACGG‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATTTTTTCGTCCTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAg                                                                                                                                                                                                                               <  1:342803/251‑1 (MQ=255)                       cGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAA‑CGTTACCCGCTT‑CGTCATACCATGCCGGGATACGG‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAgg                                                                                                                                                                                                                              >  2:458308/1‑251 (MQ=255)                                   gACGAAGGACAACATCACCACCGAGGTTATTTTTTTTACGCCCTTCTTCTTCGTTGGGGCCAACATAAC‑CGTTCCCCGCTT‑CGTCATACCAGGCCGGGATACGG‑TGACCCCCCC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCACCCTCTTCAACCGCCGGTTTCGCCAGGACATCGgcac                                                                                                                                                                                                                    <  1:252718/249‑1 (MQ=255)                                     cccggaacaaaaccgccACCGAGGTTTTTTTTTTTACGCACTTCGTTTTCGTGGCGGCAATCAAAAT‑GGTTACCCGCTG‑CGTCATACCATGCCGGGATACGG‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTCCGGTACGAAC‑TGAATGTC‑GCCGTTCTCACCCGCTTCAACCGCCGGTTTCGCCAGGACATCGgcacgcac                                                                                                                                                                                                                <  2:388155/235‑1 (MQ=255)                                          gACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAA‑CGTTACCCGCTT‑CGTCATACCATGCCGGGATACGG‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTAcc                                                                                                                                                                                                           >  1:264630/1‑251 (MQ=255)                                              acaTCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAA‑CGTTACCCGCTT‑CGTCATACCATGCCGGGATACGG‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTg                                                                                                                                                                                                       >  1:260653/1‑251 (MQ=255)                                                 ccagcccagtgttttttttttttCCGCATTTCGTCTCCGTTGCGCCCAAAATAAA‑CGTTACCCGCGT‑CGTCATACCAGGCCGGGATACGG‑TGACCCCCCC‑AAAACGGACGA‑GAGATACACCAGCCTTGAATACCGCGCATCCAGGAGAAGT‑ACATGTTTTCTCCCTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCACACCCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTggtc                                                                                                                                                                                                    <  2:201498/236‑1 (MQ=255)                                                          gggTTATTTTCTTTCCCCCTTTCGTCTTCGTTGCGGCCAACATAAC‑CGTCACCCGCTT‑CGTCACCCCAGGCCGGGAACCGG‑TG‑CCCCACCCAAAACTGACGA‑GAGATACACCAGTCTTGAATATCGCGCACCCGGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAGC‑TGAATGTC‑GCCGTTCTCACCCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCAt                                                                                                                                                                                           <  2:39284/250‑1 (MQ=255)                                                            gTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAA‑CGTTACCCGCTT‑CGTCATACCATGCCGGGATACGG‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTg                                                                                                                                                                                         >  1:285792/1‑251 (MQ=255)                                                             tttttttCTTTACGCACTCCGTTTCCGTGGCGCCCAACATAAA‑CTTTCCCCGTTT‑CCTCATCCCTTCCCGTGATCCGG‑GCCCCCCCCC‑ACAACTGACGA‑GGGATACCCCAGTCCTGAAAATCGCGCATCCAGGAAAATT‑ACATTTTTCCTTCCTGCTTCGGTCCAAAC‑TGAATGTC‑GCCGTTCACACCCCCTTCACCCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTg                                                                                                                                                                                         <  2:390504/247‑1 (MQ=255)                                                                      ttccgccctccgtcttcgtGGGGGCCAACATAAC‑CGTTCCCCCCTT‑CGTCATACCAGGCGGGGTTACGT‑GGCCCCCCCC‑ACAACTGTCGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAATT‑ACATGTTTTCTTTCTGCTTCGTTACGACC‑TGAATGTC‑GCCGTTCTCACCCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATa                                                                                                                                                                                 <  1:299174/239‑1 (MQ=255)                                                                      gtacccccttcgtcttcgtGGCG‑CCAACAAAAAACGTT‑CCCGCTTTCGCCAACCCAGGCCGGGATACGG‑GGCCCCCCCC‑ACACCGGAGAA‑GAGATACACCAGTCCTGAAAACCGCGCACCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGACC‑TAAATGTC‑GCCGTTCCCAACCGCTACAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGAt                                                                                                                                                                                  <  2:50863/241‑1 (MQ=255)                                                                        ccgctcttcgtcttcttGGCGGCCAACATAAA‑CGTTCCCCTATT‑CGTCATCCCAGGCCGTGATACGG‑TGACCCCCCC‑ACAACTGACGA‑GATATACACCACTCCTGAATATCTCGCCTCCGGGATAAGT‑CCATTTTTTCGTACTGCTTCGGTACGACC‑TAAATGTC‑GCCGTTCTCACCCGCTACAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTa                                                                                                                                                                             <  2:243519/246‑1 (MQ=255)                                                                                       gtTGCGGCCAACATAAA‑CGTTACCCGCTT‑CGTCATACCATGCCGGGATACGG‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGc                                                                                                                                                              >  1:85822/1‑251 (MQ=255)                                                                                            ggCCAACATAAA‑CGTTACCCGCTT‑CGTCATACCATGCCGGGATACGG‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTa                                                                                                                                                          >  1:144513/1‑250 (MQ=255)                                                                                                aaCATAAA‑CGTTACCCGCTT‑CGTCATACCATGCCGGGATACGG‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACACCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGa                                                                                                                                                      >  2:234062/1‑250 (MQ=255)                                                                                                    aaaaacttaccccGCTT‑CTTCATACCAGGCCGGGATACGG‑TGCCCCCACC‑ACAACGAACGA‑GAGATACCCCAGTCCTGAATATCGCGAAGCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACAACC‑TGAATGTC‑GCCGTTCTCACCCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACgg                                                                                                                                                  <  2:52459/241‑1 (MQ=255)                                                                                                     aaa‑CGTTACCCGCTT‑CGTCATACCATGCCGGGATACGG‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTc                                                                                                                                                >  1:124455/1‑251 (MQ=255)                                                                                                                   tt‑CGTCATACCATGCCGGGATACGG‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGtt                                                                                                                                   >  1:203014/1‑251 (MQ=255)                                                                                                                     tcttCATACCATGCCGGGATACGT‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGt                                                                                                                                    <  1:458308/246‑1 (MQ=255)                                                                                                                     cgcgcATACCATGCGGGGATACGG‑TGCCCCCACC‑AAAACTGAGGA‑GAGATACACAAGCCCTAAATACCGCGCATCCAGGAGAAGA‑ACTTTTTTTCGTACTGCTTCGGTAGGAAC‑TGATTGTC‑GCCGTTCTCAACCGCTTCACCCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGAGAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGttt                                                                                                                                  <  2:330280/247‑1 (MQ=255)                                                                                                                       gtaaTACCTTGCGGGGTTACGG‑TGACCCCACC‑ACACCTGAATA‑GAAATACACCATTCCTGAATATCGCGCAGCTAGGAGAAGT‑ACTTGTTTTCTTACTGCTTGGTTACAAAC‑TGAAGGTG‑GCTTTTCTCAACCGCTTGAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGGTCGATAACTACGACGCCACGGTCGACGTAAGGAACGGTCAGGTCGTGCGGTTTaa                                                                                                                                <  1:63392/248‑1 (MQ=255)                                                                                                                            cccgtGCCGGGA‑ACGGGTGACCCCCCC‑CCAACGGACGG‑GGGATACCCCAGCCCTGAAAACCGCGCACCCGGGAGAAGT‑ACATGTTTTGGTACTGCTGCGGTCCGAAC‑TGAAGGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATttct                                                                                                                           <  2:72775/247‑1 (MQ=255)                                                                                                                              cATGCCGGGATACGG‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCACCCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATttcttc                                                                                                                         >  1:165555/1‑251 (MQ=255)                                                                                                                                    gggATACGG‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCcagc                                                                                                                     >  1:416775/1‑249 (MQ=255)                                                                                                                                          ccgtggcccccccC‑ACAACTACCGG‑GAGATACCCCAGCCCTAAATTCCGCGCAACAGGGAAAATT‑ACATGTTTTCTTTCTGCTTCGGTACGACC‑TGAATGTC‑GCCGTTCCCACCCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCaa                                                                                                              <  2:287663/244‑1 (MQ=255)                                                                                                                                           gg‑TGACCCCCCC‑ACAACTG‑CGAAGAGATACACCGGGCATTAATATCGCGCACCCAGGAGAAGT‑CCCTTTTTTCGTACGGCTTCGGTACGACC‑TGAATGCC‑GCCTTTCTCACCCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAgc                                                                                                            <  2:408735/251‑1 (MQ=255)                                                                                                                                           gg‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAgc                                                                                                            >  2:285768/1‑251 (MQ=255)                                                                                                                                             ggtgccccccccaaaccTGACGG‑GAGATCCCCCAGCCCTAAAAATCGCGCATCAAGGAGAAGT‑ACTTGTTTCCGTCCTGCTCGGGTACAAAC‑TAAATGTC‑GCCTTTCCCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAgc                                                                                                            <  2:333854/235‑1 (MQ=255)                                                                                                                                                cccccACC‑ACAACTGACGA‑GAGATCCACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGt                                                                                                        <  1:17674/250‑1 (MQ=255)                                                                                                                                                     ccc‑acAACTGACGA‑GAGATACCCCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACTTGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTCCTCACCCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACt                                                                                                   <  1:186767/250‑1 (MQ=255)                                                                                                                                                       cccccaccggacagaagaaacccctgTCCTGAATATCGGGCACCCGGGA‑AAGTTAAATGTTTCCGTACTGCTTCGGT‑CGAACCGGAATTCC‑CCCGTCCTCACCCGCTACAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACt                                                                                                   <  2:260653/226‑1 (MQ=255)                                                                                                                                                            aCTGACGA‑GAGATGCACCGCCCCTGAATACCGCGCACCCAGGAGAAGT‑ACATGTTTCCTTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCACCCCCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCa                                                                                             <  2:357600/251‑1 (MQ=255)                                                                                                                                                                gacgaggtatacccaaGTCCTGAATATCGCGCCCCCAGGGGACGT‑ACATGTTTTCTTACTGCTTCGGTCCGAAC‑GGAATGCC‑GCCGTTCTCACCCGCTTCAGCCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAAGTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAAc                                                                                           <  2:132979/235‑1 (MQ=255)                                                                                                                                                                  ga‑gagaTACCCCGGTCTGGAATACCGCGCACCCAGGAGAAGT‑CCTTTTTTTCGTCCTGCTCCGTTACGAAC‑TGAATGTC‑CCCGTTCTCAACCCCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACt                                                                                       <  2:264630/251‑1 (MQ=255)                                                                                                                                                                       gaTACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCtt                                                                                   >  2:37585/1‑251 (MQ=255)                                                                                                                                                                               aGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTgc                                                                            >  1:338686/1‑250 (MQ=255)                                                                                                                                                                                 tCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTCAGGAACGGTCAGGGCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCCACTGCGGCAACGACTGCTTTACGTgcagc                                                                         >  2:440648/1‑251 (MQ=255)                                                                                                                                                                                       atatCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCTCCACGGTCGACGCAAGGAACGGTCGGGTCGTGCGGTTTAGTTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCACAGACTGCTTCACGTgtcgcaacacg                                                                   >  2:225997/1‑241 (MQ=255)                                                                                                                                                                                       ataccgcgCACCCAGGAGAAGT‑ACTTGTTTCCGTCCGGCTGCGAAACGAAC‑TGAA‑GTCCGCCGTTCTCACCCGCTTCACCCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACg                                                                   <  2:254454/247‑1 (MQ=255)                                                                                                                                                                                       aaaaccgcgcaacaaggaaaaggtcATGTTTGCGTACTGTTTCGTTGCGCAC‑GGAATGCC‑GCCTTTCCCACCCGCCTAATCCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAAc                                                                    <  2:203014/227‑1 (MQ=255)                                                                                                                                                                                         atCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCt                                                                 >  1:183949/1‑251 (MQ=255)                                                                                                                                                                                             gcgcgcccaggagagttacagtttttCGTCCTGCTTCTGTCGGAAC‑TGACTTCC‑GCCTTTCTAAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCa                                                              <  2:416775/229‑1 (MQ=255)                                                                                                                                                                                                aTCCAGGAGCAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGttt                                                          >  2:135725/1‑251 (MQ=255)                                                                                                                                                                                                  ccAGGAGAAGC‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCCCCAGTTTCt                                                        >  1:80588/1‑251 (MQ=255)                                                                                                                                                                                                  ccAGGAAAAGT‑ACATGTTTTCGTCCTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCACCCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCt                                                        <  1:173141/251‑1 (MQ=255)                                                                                                                                                                                                     ggAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTgg                                                      >  2:263509/1‑250 (MQ=255)                                                                                                                                                                                                     ggAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGGAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTgg                                                      >  2:50461/1‑250 (MQ=255)                                                                                                                                                                                                        gaAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAAct                                                  >  2:216675/1‑251 (MQ=255)                                                                                                                                                                                                        gaAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAAcc                                                  >  2:183010/1‑250 (MQ=255)                                                                                                                                                                                                          aGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAActct                                                >  2:290625/1‑251 (MQ=255)                                                                                                                                                                                                          aGT‑ACATGTTTTCGTACTGCGTCGGACCGAAC‑GGAATGTC‑GCCGTTCACAACCCCTACAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCGTTACGTGCAGCAAAACGCTCCAGTTTCTGGAActct                                                <  2:109056/251‑1 (MQ=255)                                                                                                                                                                                                           gT‑ACATGTTTTCGGACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAActatg                                               >  2:236413/1‑248 (MQ=255)                                                                                                                                                                                                                attttttCGTCCGGCTGCGGCACGAGC‑TAAATGTC‑GCCGTTCCCAACCGCTCAAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAgg                                           <  2:102800/248‑1 (MQ=255)                                                                                                                                                                                                                   tttttctgcctgtgtcggccctatttGAATGCC‑GCGGTTCCCACCCGGTCCAACCGCCGGTTTCGCCAGGACATCGGCCCGCACGGACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCGGGAACTCTGCAggg                                          <  2:413799/225‑1 (MQ=255)                                                                                                                                                                                                                          ccTGCTTAGGTACGAAG‑TGCATGTC‑GCCGGTCTCAACCGCTTCCACCGCCGGTTTCGCCAGGACATCGGCGCGCACGTACCACTGGTCGACCAGCATTGGTTCGATAACTACGCCCCCCCGGTCGCCGTAAGGAACGGTTCGTTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCCTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATATCGCTg                                 >  2:190371/2‑251 (MQ=255)                                                                                                                                                                                                                           cTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTgg                                >  1:52169/1‑251 (MQ=255)                                                                                                                                                                                                                                 cGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATaa                           >  2:83900/1‑250 (MQ=255)                                                                                                                                                                                                                                         aC‑TGAAGGTC‑GCCGTTCTCACCCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGAt                  <  1:35659/251‑1 (MQ=255)                                                                                                                                                                                                                                            tGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTc                >  1:15796/1‑251 (MQ=255)                                                                                                                                                                                                                                                 gTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTAc           >  1:255166/1‑251 (MQ=255)                                                                                                                                                                                                                                                        gTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTgg     >  1:42060/1‑251 (MQ=255)                                                                                                                                                                                                                                                             cAACCGCTTCACCCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTAt  <  1:78939/249‑1 (MQ=255)                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                       GAAGAGAACCAGGTATCGAGAACGTCTTCGTCCTGACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAA‑CGTTACCCGCTT‑CGTCATACCATGCCGGGATACGG‑TGACCCCACC‑ACAACTGACGA‑GAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGT‑ACATGTTTTCGTACTGCTTCGGTACGAAC‑TGAATGTC‑GCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGCCGCCACGGTCGCCGTAAGGAACGGTCAGGTCGTGCGGTTTAATTTCTTCCAGCAGGCCAAGCGCGTCAACTGCGGCAACGACTGCTTTACGTGCAGCAAAACGCTCCAGTTTCTGGAACTCTGCAGGGATTTCGCTGGAATAAACGTCAGATTCGTTACCTTTGGTAT  >  NC_000913/4482402‑4482894

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