Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,331,890	C→T	L463L (CTC→CTT)	proP →	osmolyte:H(+) symporter ProP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,331,890	0	C	T	97.6%	130.2 / NA	43	L463L (CTC→CTT)	proP	osmolyte:H(+) symporter ProP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (12/29);  minor base G (0/1);  total (12/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTT  >  NC_000913/4331644‑4332130                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   gACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGTCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCttgt‑cg                                                                                                                                                                                                                                               <  1:341055/251‑1 (MQ=255) gACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCttgt‑cg                                                                                                                                                                                                                                               <  1:227154/251‑1 (MQ=255)               aCAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGataa                                                                                                                                                                                                                                  <  2:183307/250‑1 (MQ=255)                cAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGataata                                                                                                                                                                                                                                <  2:249206/251‑1 (MQ=255)                cAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGataata                                                                                                                                                                                                                                <  2:110159/251‑1 (MQ=255)                    gcgcTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGa                                                                                                                                                                                                                            <  1:167479/251‑1 (MQ=255)                      gcTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGc                                                                                                                                                                                                                          <  2:229812/251‑1 (MQ=255)                         ggcggcggcATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCa                                                                                                                                                                                                                         <  1:434087/249‑1 (MQ=255)                           cggcggcATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGaaa                                                                                                                                                                                                                     <  1:501223/251‑1 (MQ=255)                              cggcATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATc                                                                                                                                                                                                                  <  1:332781/251‑1 (MQ=255)                                gcATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCga                                                                                                                                                                                                                <  2:249535/251‑1 (MQ=255)                                    ttAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGc                                                                                                                                                                                                                          <  2:344874/237‑1 (MQ=255)                                     tAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCgatg‑at                                                                                                                                                                                                           <  1:31513/250‑1 (MQ=255)                                      aaTATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑atat                                                                                                                                                                                                         <  2:472123/251‑1 (MQ=255)                                             cGGTGCTTGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCa                                                                                                                                                                                                   <  1:94176/250‑1 (MQ=255)                                                ggCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCACgag                                                                                                                                                                                               <  1:264050/250‑1 (MQ=255)                                                 gCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCACgag                                                                                                                                                                                               <  2:576300/250‑1 (MQ=255)                                                                            ggcggcCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGAGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGCCCACGAGATTGCCGATTTGCAGGCGAAACGTAccc                                                                                                                                                                   <  2:401821/251‑1 (MQ=255)                                                                                         gTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTgcagca                                                                                                                                                      <  1:370671/251‑1 (MQ=255)                                                                                                      aGAATCTGATGATGCCGGCATATTCCCTGATGGTAGTGGCGGTGGTTGGTTTAAGCACCGGCGTAACCATGAAAGAGACGGCAAAGCGTCCGT‑TG‑AAAGGGGCGACACCGGCGGCGTAAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAAtt                                                                                                                                         <  2:438136/251‑1 (MQ=255)                                                                                                             gatgatGCTTGCCTATTACCTGATGGTAGGGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑GG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAAGCGATG‑ATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAAt                                                                                                                                  <  1:260804/251‑1 (MQ=255)                                                                                                                                              gttgtttggttgaagccccGGCGAAACCATAAAAGAGACGGCAGATCGTCCGT‑TG‑AAAGGCGCGACGCCGGCGTCGTCAGATAAAAAGGAAGCGAGGGAAATACTGTT‑CGAACATGACGATAATATCG‑GTAGTAAAAAGTTGCAGATTGCCCGCGAAAGTGCCAATTTGACGGCGAAACGTCCCCTCCTGGTGCAGCAACATCCGCGAAGTGAGGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGtt                                                                                                 <  2:269483/234‑2 (MQ=255)                                                                                                                                                           aTCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAAt                                                                                    >  1:455964/1‑251 (MQ=255)                                                                                                                                                             cACCGGCGTAACCATAAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCATATATACAGGAAGCGAAGGAAATTTTTGG‑CGAGCATTACGATAATATCGAGCAGAAAATCGTTG‑ATTTTGACCACGAGATTGCCGATTTGCAGGCGAAACTTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCTTGATGTGACGCTGTCTTATCAGGCCAAtt                                                                                   <  1:201055/250‑1 (MQ=255)                                                                                                                                                                cggcgaacccATGAAAGAGAAGGCAAATCGTCCGT‑TG‑CAAGGTGCGACACCGGCGGCGTCAGATCTACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGCTAATATCGAGCAGAAAATCGATG‑ATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAAc                                                                               >  2:38107/9‑251 (MQ=255)                                                                                                                                                                cGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAAc                                                                               >  1:512292/1‑251 (MQ=255)                                                                                                                                                                      aaCCATGAAAGAGCCGGCAAAGCGTCC‑T‑TGTAAAGGTGGGACACCGGCGGCGTCAGATAGACAGGAAGCGAAGGAAATTCT‑GTGCGGGCATGAGGATAATATCGAGCAGAAAATCGATG‑ATATTGAGCACTGGATTGCCTATTGGCAGGCGAAGCGTACCCGCCTGGTGCATCAACATGCGCTGATGTATGAATATGTTGAAGCGGATGGCCTGATGTGGCGCTGTCTTATCAGGCCAATTGAACTCTTa                                                                          <  1:56699/250‑1 (MQ=255)                                                                                                                                                                          aTGAAAGGGACGGAAAAACGTC‑GTGTG‑AAAGGGGCGATCCCGGGGGCGTCAGATAAACAGAAAGCGAAGGAAATGTTTGT‑AGAGCATTACGATAATATCGAGCAGAAAAGCGAGA‑ATATTAACAACGAGATTGCCGTTTTGAGGGCGAAAAGAGCCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGtt                                                                     <  2:423992/251‑1 (MQ=255)                                                                                                                                                                              aaGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACt                                                                 <  1:532817/251‑1 (MQ=255)                                                                                                                                                                                 agaCGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGCATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGTCCAATTGAACTCTTAAGGTTCACTTa                                                               >  2:575735/1‑250 (MQ=255)                                                                                                                                                                                 agaCGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCACGAGATCGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTa                                                               >  2:397699/1‑250 (MQ=255)                                                                                                                                                                                            cGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGAcgcg                                                   <  2:386679/251‑1 (MQ=255)                                                                                                                                                                                                  gtgg‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATAc                                              <  2:537224/248‑1 (MQ=255)                                                                                                                                                                                                        aaGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTctcc                                         >  1:299092/1‑251 (MQ=255)                                                                                                                                                                                                                       cggcgTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACgg                            >  2:455464/1‑249 (MQ=255)                                                                                                                                                                                                                       cggcgTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACgg                            >  2:72685/1‑249 (MQ=255)                                                                                                                                                                                                                        ggcgTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACggagg                         >  2:322549/1‑251 (MQ=255)                                                                                                                                                                                                                          cgTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTCCCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGgaggag                       <  2:234599/251‑1 (MQ=255)                                                                                                                                                                                                                           gTCAGATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCTGGTTAACggagg                         >  1:231141/1‑248 (MQ=255)                                                                                                                                                                                                                               gATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAg                                                                                                                                                                                                                           >  1:139902/1‑51 (MQ=255)                                                                                                                                                                                                                               gATATACAGGAAGCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAg                                                                                                                                                                                                                           <  2:139902/51‑1 (MQ=255)                                                                                                                                                                                                                                           gCGAAGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCg      <  1:416955/251‑1 (MQ=255)                                                                                                                                                                                                                                             aaaGGAAATTCTTGT‑CGAGCATTACGATAATAGCGAGCAGAAAATCGATG‑ATATTGACCACGAGATTGCCGATTTGCAGGCGAAACTTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCgtg    <  1:524546/250‑1 (MQ=255)                                                                                                                                                                                                                                              aaGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCgtgt   >  2:314818/1‑251 (MQ=255)                                                                                                                                                                                                                                               aGGAAATTCTTGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGtt  <  1:173815/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   GACGCATATCCGTTACAGCGCGCTGGCGGCGGCATTTAATATTTCGGTGCTGGTTGCCGGTCTGACGCCAACGCTGGCGGCCTGGCTGGTCGAAAGCTCGCAGAATCTGATGATGCCTGCCTATTACCTGATGGTAGTGGCGGTGGTTGGTTTAATCACCGGCGTAACCATGAAAGAGACGGCAAATCGTCCGT‑TG‑AAAGGTGCGACACCGGCGGCGTCAGATATACAGGAAGCGAAGGAAATTCTCGT‑CGAGCATTACGATAATATCGAGCAGAAAATCGATG‑ATATTGACCACGAGATTGCCGATTTGCAGGCGAAACGTACCCGCCTGGTGCAGCAACATCCGCGAATTGATGAATAAGCTGAAACGGATGGCCTGATGTGACGCTGTCTTATCAGGCCAATTGAACTCTTAAGGTTCACTTAATCTCTGACGCGCATACTCTCCTCCAGGTTAACGGAGGAGAGTGCAATGAAAAACCGTGTT  >  NC_000913/4331644‑4332130

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