Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	954,388	C→T	23.8%	A27T (GCA→ACA)	focA ←	formate channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	954,388	0	C	T	23.8%	41.4 / 4.3	37	A27T (GCA→ACA)	focA	formate channel
Reads supporting (aligned to +/- strand):  ref base C (18/10);  new base T (2/7);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 5.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

TGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACG  >  NC_000913/954154‑954591                                                                                                                                                                                                                                                  |                                                                                                                                                                                                             tGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATTGCC‑CGATAGAAAAGCTTAACGGGAATGTTT‑GTTGCtttt                                                                                                                                                                                                          >  1:1135772/1‑238 (MQ=255) tGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATAGTC‑AGCTAGAAAAGCTTAA‑ACGGGAGTTTCGTGGCTTta                                                                                                                                                                                                          >  1:490041/1‑239 (MQ=255)    cccccCTCCCCTTTACAACCAACATCAACCCGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACAC‑AAAGACACAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAg                                                                                                                                                                                                       <  2:1094599/238‑1 (MQ=255)                  gCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGG‑TATTCGCC‑AGGTAGAAAAGCTTTA‑GCGGATGTTTCTTTGCTTTATAGA‑CACCCCCCTCTTc                                                                                                                                                                                        >  1:937656/1‑239 (MQ=255)                                     tGGAAGTAAAGAGATCGGCTCCGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGG‑TAAACGCC‑ACATAGAAAGTCTTAA‑ACGGATGTTTCGTTGCTTTTTAGA‑CACCCCCCCCTTCGGCACCTTTGTCCCTTcca                                                                                                                                                                     >  1:818736/1‑239 (MQ=255)                                         ataaaaGAGCCCGGCCCCCCACACCACAC‑AAAAAAACACCCCAAGAGAGAAGCAAAGGCCCCCAACCATCTTTG‑CCAGGACGAAGGTCTTTGT‑GCCGGTGCCAGTGGTTGCTTTGATAGAGAAGACGAATGGTATTGAG‑ATAAAAACCCGGGGGG‑TAATCGCC‑AGAAAGAAAGTCTTAA‑GGGGATGTTTCGTTGCTTTATAGA‑CCCCCGCCTCTTCGGCCACTTTGGCCATTgcagca                                                                                                                                                                  <  2:1134153/235‑1 (MQ=255)                                          agacaagagttggCTCCTCAGAAAAAATAAAAGAATCAGCCCGTGAGAAAACCAAATGCCGCCAACTATTTTGC‑CCAGGCCGAAGGGCGTGGG‑CCATGTCGATGTGGTGGCTGGTATAAAGAAGAGGAGTGGAATTTAA‑ATGAAATCCCCGGGGG‑TACTCCCA‑AGATAAAAATTCGTAT‑GCGGAGTTTTCGTTGTTTTAGAAA‑CCCCCGCCTCTTCGGCCACTTTGGCCATTgcagca                                                                                                                                                                  <  2:185993/227‑1 (MQ=255)                                             aaGAGATCGGCTCCGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAG‑ATGAAAAAACCGGCGG‑TAATCGCC‑AGATAGAAAGTCTTAC‑GCGGATGGTTCGTTGTCTTCTAGA‑CACACGCCTCTTCGGCCACCTTTGCCAtttctcccggtag                                                                                                                                                             >  1:548219/1‑228 (MQ=255)                                                  aTCGGCTCCGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGCG‑ATGAAAACACCGGCGG‑CAAATGCC‑CCGTAGAAAATCTTAA‑ACGGAAGTTTAGTTGTTTTATAAA‑AACCCGCCCCTTCGTCCCCTTTTGACACTtcagcaagggaaaaaa                                                                                                                                                        >  2:1028746/1‑229 (MQ=255)                                                    cGGCTCCGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAAGCTTCA‑GCGGCTGGTTCGTTGCTTTATAGA‑CACCCGCCCCTTCGGCCACTTTTGCCCTTTCAGCAGGAAGTaaaaaa                                                                                                                                                      >  1:657322/1‑237 (MQ=255)                                                       cTCCGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATAGCC‑AGCTTGAAAGGCTTAA‑GACGGTGTTTCGTTGCTTTTTAGA‑CACCCGCCTCTTCGGCCCCTTTTGGCATTGCAGCAGGCAGTAAAAGAtaa                                                                                                                                                   >  1:35702/1‑237 (MQ=255)                                                          cGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGG‑GAATCGCC‑AGATAGAAAGGCTTAA‑GCGGATGTTTCGTTGCTTTATAGA‑CACCCGCCCCTTCGGCCCCTTTGGCCATTGCAGCAGGAAGTCAAAGATAaaaa                                                                                                                                                >  1:269483/1‑239 (MQ=255)                                                                 aacaC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAGTCTTAA‑GCGGATGGTTCGTTGCTTTATAGA‑AACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGGAAAAAATCAAAAGGGGTg                                                                                                                                          >  1:234073/1‑238 (MQ=255)                                                                 aacaC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAGGCTTAA‑ATGGAGGGTTCGTTGCTTTATAGG‑CACCCGCCCCTTTGGTCACCTTGGCCATTGCAGCAGGAAGTAAAAGACGAAAAGGGGTgc                                                                                                                                         >  1:692570/1‑238 (MQ=255)                                                                             tcgccccaAGAGGAAAGCAAGTGCCGCCAACTGTTTTTG‑CCAGGCCGAGGGGTTTTGG‑CCCTGTGCCAGGGGTTGCTGTAAAAAAGAAGACGAATGCGA‑TGAGAATAAAACCCCCGGCGG‑TAACCGCC‑AAATAGAACGTCTTAA‑GCGGATGTTTCGTTGTTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTcaca                                                                                                                              <  2:269483/235‑1 (MQ=255)                                                                              cAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATGGCC‑AGATAGAAAGTCTTAA‑GCGGAAGTTTCGTTGCTTTATAGA‑AACCCCCCCCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAAATCAAAAGGGTTGTAACCTTTCacac                                                                                                                             >  1:796566/1‑239 (MQ=255)                                                                                 ccccAGAGAGAAGCAAATGCCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGGTAGAAAGTCTTAA‑ACGGATGTTTCGTTGCTTTATAGA‑AACCCGCCCCTTCGGCCACCTTTGCCCTTGCAGCCGGAAGGAAAAAATTCAAAGGGGTTTCAGCCTTTCCACTaa                                                                                                                          >  1:1268126/1‑239 (MQ=255)                                                                                           aaGCAAATGCCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAAGCTTAA‑GCGGATGTTTCGTTGCTTTTTAAA‑TCCCCCCCCCTTCGGCCCCTTTTGCCATTGCCGCCGGAAGAAAAAAATCAAAAAGGTGGTCACCCTTTCCCCTAAATCCCTCttt                                                                                                                >  2:776548/1‑239 (MQ=255)                                                                                           aaGAAAATGCCGCCAACCAGTTTTG‑CCATGCCGCAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGCTATAGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGG‑TCACTACC‑CAATAGACAAATTTTA‑GCGAGAGTTTCCTTGCATTATTCG‑CACCACGCTCTTTGTCCACTTTTGCCCTTTCAACCAGAGATAAAAAATTCAAAGGGTTATCAGCCTTCAACAAAACTCTCCCttt                                                                                                                >  2:290818/1‑239 (MQ=255)                                                                                                 aTGCCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCCCTTTGGCCCTTTCAGCCGGAAGTAAAAAATTAAAAGGGGTGTTAGCTTTTACCCTAACTCTCCCTTTATTAAg                                                                                                          >  1:612697/1‑239 (MQ=255)                                                                                                         aacagtttttG‑CCATG‑CGGAGGGGATTGTGGCCTGTGCCGTTGGTGGCTGTGATATAGAAGCCGAATGCGATTGAA‑ATGAAAACACGGGCGG‑TAACCGCC‑AGATAGAAAGTTTTAA‑GCGGATGTTTCGTTGCTTTATAGC‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGcgac                                                                                                  <  1:477529/233‑1 (MQ=255)                                                                                                          acatgTTTTGCCCATGCGAAAGGGTTGTGT‑CCTTGGGACAGTGTTTGCTGAGATAAAGAAGAAAATTGGGATGAAA‑ATCAAAACCCCGGGGT‑TAATCCCA‑AGATAGAAACTTTTAA‑GGGGTTGTTTCGTTGCTTTATAGA‑CCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGcgac                                                                                                  <  2:1268126/235‑1 (MQ=255)                                                                                                                         gCCAAAGGGCATTGT‑GCCTGTGCGAGTGGTTGCTGTGAAATAGAAGACAAATGCGATTGAG‑ATGAAACCCCGGGCGG‑TACTCGCC‑AGATAGAATGTCTAAA‑GCGGATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAg                                                                                   <  1:351088/239‑1 (MQ=255)                                                                                                                                   aTTGG‑GCCGGCGCAAGGGGTTTCTGTGATAAAGAAGGCGAAGGGGATTGAA‑ATGAAAACACGGGGGT‑TAATCGCA‑AGATAGAAAGTCTTAG‑GGGGATGTTTCGTTGTTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACa                                                                                      <  1:334579/226‑1 (MQ=255)                                                                                                                                      tgggcctggaccAGTGGTTGTTGTAATATAGAAGCAGAAGGCGATTGAA‑ATGAAACCCCCGGCGG‑TAACCGCA‑AGAAAGAAAGTCTTAA‑GGGGTTGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATga                                                                       <  1:776548/236‑1 (MQ=255)                                                                                                                                           tgtggctggtggttgtggttATAAAGAAGACGAATGAGATAGAA‑ATAAAAGCCGCGGAGG‑TAAAGGCG‑AAGAAAAAAGTCTAAA‑GCGGTTGTTTTGTGGTTTTAACAA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGGCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGaaa                                                                   <  2:548219/230‑1 (MQ=255)                                                                                                                                             gtgccAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑CGATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGa                                                                                                                                                      >  1:121102/1‑153 (MQ=255)                                                                                                                                                 aaggggTTGCTTGAATATAAAAGACGAATGCGATTGAG‑ATAAAAACCCCGGGGG‑TAATCGCA‑AGATAGAAAGTCTTAA‑GCGGTTTTTTTGTTCCTTTA‑AGACCACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGatat                                                            <  1:701739/235‑1 (MQ=255)                                                                                                                                                     ggTTGCTGTAAAATAGAAGAAGAATGGGTTTGAG‑AT‑AAAACCCCGGGGGTTAACCCCC‑AGATAGAAAGTTTTAA‑GCGGATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGAt                                                        <  1:542082/239‑1 (MQ=255)                                                                                                                                                        tGCTGTGATATAGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGCAGTAAAAGATCAAAAAGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCat                                                     >  1:844605/1‑239 (MQ=255)                                                                                                                                                                  agaaagcagAATGCGAGTGAG‑ATAAACACCCGGGGGG‑TAATCGCA‑AGATAGAAAGTCTTAA‑GGGGTTGTTGCGTGGTTTTATAGA‑CCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTg                                           <  1:195635/231‑1 (MQ=255)                                                                                                                                                                       gaCGAATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAATCTTAA‑GCGGATGTTTCGTTGCTTTATAGA‑AACCCGCCTCTTCGGCCACTTTGGCCATTTCAGCAGGAAGTAAAAGATCAAAAGGGGTGTCAGCTTTTACACTAACTCTCTCTTTTTTAATTCGGCGGCGGGATACTAACAAAACATTTAAAATGAGAAAATGATATAGAACATAAATCCCCTGGCTta                                      >  2:1372602/1‑239 (MQ=255)                                                                                                                                                                        aCGAAGGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGata                                                             <  1:388841/215‑1 (MQ=255)                                                                                                                                                                                 aTTGAG‑ATAAAACCCCGGGGGG‑AAACCGCC‑AGAAAGAAA‑TGTTAAGGGGGTTGTTTCGTGGTTTTAAAGC‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGt                             <  1:1014338/238‑1 (MQ=255)                                                                                                                                                                                           aaaaCCCCGGGGG‑TAAACCCCAAAATAAAAATTCTTAA‑GCGGTTGTTTGGTTGTTTTAAAGA‑CCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCAt                    <  1:1372602/239‑1 (MQ=255)                                                                                                                                                                                             aaCACCGGCGG‑TAATCGCC‑AGATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACCTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGGTGTCAGCTTTCACACTAAACCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCCTTATAGATTAGAAACTGATATAGATCatat                                                   >  1:1138186/1‑205 (MQ=255)                                                                                                                                                                                                             cGCA‑AGATAGAAAGTCTTAA‑GGGGTGTTTTCTTTTCTTTATAGC‑CCCCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACg  <  1:770150/239‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                             TGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACAC‑AAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTG‑CCATGCCGAAGGGCATTGT‑GCCTGTGCCAGTGGTTGCTGTGATATAGAAGACGAATGCGATTGAG‑ATGAAAACACCGGCGG‑TAATCGCC‑AGATAGAAAGTCTTAA‑GCGGATGTTTCGTTGCTTTATAGA‑CACCCGCCTCTTCGGCCACTTTGGCCATTGCAGCAGGAAGTAAAAGATCAAAAGGGTTGTCAGCTTTCACACTAACTCTCTCTTTATTAAGTCGGCGACGAGATACTAACAAAGCATTATAGATGAGAAATTGATATAGATCATATCTCGCCTGGCTTATAGGCCCGTAACTCGCATGGTTTTTATGCAAATACG  >  NC_000913/954154‑954591

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