Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	2,231,588	G→C	20.9%	pseudogene (595/1526 nt)	mdtQ ←	multidrug resistance outer membrane protein MdtQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	2,231,588	0	G	C	20.9%	91.1 / 2.7	37	pseudogene (595/1526 nt)	mdtQ	multidrug resistance outer membrane protein MdtQ
Reads supporting (aligned to +/- strand):  ref base G (7/21);  new base C (5/3);  total (13/24)
Fisher's exact test for biased strand distribution p-value = 8.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑T‑‑TTG‑CTCGCGTTCCGCCG‑CCCGTG‑CTTTAACCGTACCCAGGCGG‑GCA‑GTAACCTCCGCCCGATTCTTT‑CCCCAGAT‑ATCGAGA‑TGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTA  >  NC_007779/2231358‑2231825                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    gctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑T‑‑TTG‑CTCGCGTTCCGCCG‑CCCGTG‑CTTTAACCGTACCCAgg                                                                                                                                                                                                                                             <  1:616198/239‑1 (MQ=255)  ctgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑G‑‑TTG‑CTCCCGCTCCCCCGCCCCGGC‑CTTTAACCGAAcccggc                                                                                                                                                                                                                                             >  1:1937447/1‑236 (MQ=255)   tgcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑T‑‑TTG‑CTCGCGTTCCGCCG‑CCCGTG‑CTTTAACCGTACCCAGGc                                                                                                                                                                                                                                            <  1:1395734/238‑1 (MQ=255)    gcCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑T‑‑TTG‑CTCGCGTTCCGCCG‑CCCGTG‑CTTTAACCGTACCCAGGCgg                                                                                                                                                                                                                                          <  1:891693/239‑1 (MQ=255)                   cATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑T‑‑TTG‑CTCGCGTTCCGCCG‑CCCGTG‑CTTTAACCGTACCCAGGCGG‑GCA‑GTAACCTCCGc                                                                                                                                                                                                                          <  2:474748/238‑1 (MQ=255)                    aTTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑T‑‑TTG‑CTCGCGTTCCGCCG‑CCCGTG‑CTTTAACCGTACCCAGGCGG‑GCA‑GTAACCTCCGccc                                                                                                                                                                                                                        <  2:612138/239‑1 (MQ=255)                    aTTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑T‑‑TTG‑CTCGCGTTCCGCCG‑CCCGTG‑CTTTAACCGTACCCAGGCGG‑GCA‑GTAACCTCCGccc                                                                                                                                                                                                                        <  1:1471552/239‑1 (MQ=255)                    aTTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCGCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑T‑‑TTG‑CTCGCGTTCCGCCG‑CCCGTG‑CTTTAACCGTACCCAGGCGG‑GCA‑GTAACCTCCGccc                                                                                                                                                                                                                        <  1:1128292/239‑1 (MQ=255)                       aaTATCGTTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑T‑‑TTG‑CTCGCGTTCCGCCG‑CCCGTG‑CTTTAACCGTACCCAGGCGG‑GCA‑GTAACCTCCGCCCGAt                                                                                                                                                                                                                     <  1:772753/239‑1 (MQ=255)                                ttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑T‑‑TTG‑CTCGCGTTCCGCCG‑CCCGTG‑CTTTAACCGTACCCAGGCGG‑GCA‑GTAACCTCCGCCCGATTCTTT‑ccc                                                                                                                                                                                                            <  2:931884/238‑1 (MQ=255)                                ttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑T‑‑TTG‑CTCGCGTTCCGCCG‑CCCGTG‑CTTTAACCGTACCCAGGCGG‑GCA‑GTAACCTCCGCCCGATTCTTT‑ccc                                                                                                                                                                                                            <  1:2235810/238‑1 (MQ=255)                                ttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑T‑‑TTG‑CTCGCGTTCCGCCG‑CCCGTG‑CTTTAACCGTACCCAGGCGG‑GCA‑GTAACCTCCGCCCGATTCTTT‑ccc                                                                                                                                                                                                            <  1:607214/238‑1 (MQ=255)                                ttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑T‑‑TTG‑CTCGCGTTCCGCCG‑CCCGTG‑CTTTAACCGTACCCAGGCGG‑GCA‑GTAACCTCCGCCCGATTCTTT‑ccc                                                                                                                                                                                                            <  1:485591/238‑1 (MQ=255)                                ttCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑T‑‑TTG‑CTCGCGTTCCGCCG‑CCCGTG‑CTTTAACCGTACCCAGGCGG‑GCA‑GTAACCTCCGCCCGATTCTTT‑ccc                                                                                                                                                                                                            <  1:577342/238‑1 (MQ=255)                                  ccACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑G‑‑TTG‑CGCGCGGTCCTCCC‑CCCCGG‑CGTTAACCCGACCCCCGCGG‑GGG‑GGAACCCCTCCCCGCTTCTTT‑tccccc                                                                                                                                                                                                         >  1:1088364/1‑237 (MQ=255)                                    acacCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑T‑‑TTG‑CTCGCGTTCCGCCG‑CCCGTG‑CTTTAACCGTACCCAGGCGG‑GCA‑GTAACCTCCGCCCGATTCTTT‑CCCCAGat                                                                                                                                                                                                       <  1:1143415/239‑1 (MQ=255)                                             aCTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCCACGCTGCCCGCCCGCCATTGGCGGG‑TCGTTG‑CTCGCGTCCCGCCC‑CCCGTT‑TTTTACCTCTCTCCGGGGGG‑GAC‑GAAACCCCCCCCCGTCTTTTT‑CCCCAGAT‑CGCgag                                                                                                                                                                                                >  2:1623711/1‑238 (MQ=255)                                              cTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑T‑‑TTG‑CTCGCGTTCCGCCG‑CCCGTG‑CTTTAACCGTACCCAGGCGG‑GCA‑GTAACCTCCGCCCGATTCTTT‑CCCCAGAT‑ATCGAGA‑TGc                                                                                                                                                                                           <  1:1374901/239‑1 (MQ=255)                                                               ttCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑G‑‑TTG‑CTCGCGTTTCGCCG‑CCCGGG‑CTTTTACCCCACCCCGGGGG‑GGC‑GGAAACTCCGCCCGCATCTTT‑TCCCCGAT‑ATACAGA‑TGCCCGcccccccgtaatc                                                                                                                                                                           >  1:1251236/1‑231 (MQ=255)                                                                tCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑T‑‑TTG‑CTCGCGTTCCGCCG‑CCCGTG‑CTTTAACCGTACCCAGGCgg‑g                                                                                                                                                                                                                                        <  1:2053559/180‑1 (MQ=255)                                                                    atatAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑G‑‑TTG‑GTCCCCGTTCGCCC‑CCCCTGCCTTTACCCCGCCCCAGGGGG‑GGC‑GAAACCCCCCCACGGTGCTTT‑‑CCCAGCT‑CTCGAGAGTGCCGGCCCGCCGtccagacaaca                                                                                                                                                                      >  2:535143/1‑229 (MQ=255)                                                                               ggggggtggggATATTGGTTTTCTTCTTTTTTTCTATTTGCGGCAGGACCGTTTTTAGCGCCCCGGGGTTTTGCCCCCCCCAGTCCAGGCGGGCTCCGCTGCCAGCCAGCAATGGGGGGG‑T‑‑TTG‑CTCGCTTCCCCCCG‑CCCGTT‑TTTTAACCGTCCCCGGGGGG‑GCA‑GTACCCCCCCCCCGATTCTTT‑CCCCAGAT‑ATCGAGA‑TGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGgt                                                                                                                                                           <  2:1122474/228‑1 (MQ=255)                                                                                     tCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑T‑‑TTG‑CTCGCGTTCCGCCG‑CCCGTG‑CTTTAACCGTACCCAgg                                                                                                                                                                                                                                             <  1:1659456/155‑1 (MQ=255)                                                                                                                                tgtgttacccgccgcGTGGGTTTCCCCCCCCCGGCAGGGGGGGGCCCCGCTGCCCCCCAG‑AATGGGGGGGTT‑‑TTC‑CCCCCGTCCCCCCG‑CCCGTG‑CTTAACCCGTCCCCGGGGGG‑GCA‑GTCCCCCCCCCCCGTTTCTTC‑CCCCAGAT‑ATCGAGA‑TGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGc                                                                                                         <  2:754873/229‑1 (MQ=255)                                                                                                                                  tgtTTAACGCCGCCTGGGTTTGCCACTCCCAGTACTGGCGGGCTACGCTGCCAGCCCGCAATTGGCGGG‑G‑‑TTT‑GTCCCGGTCCCCCC‑CCCCCGCCTTTAACCCACCCCCGGGGG‑GCG‑GAACCCCCCCCCCCCTTTTTT‑CCCC‑CAA‑ATCCAGAGGGCCACCCCCCCCTTAACCCAAAAAAACCCTTGGGGAA‑CCCCGGGCCGG‑CGGACCCCCGGCCGGGGCCGTCTGAGGAAACAGccc                                                                                                       >  2:687190/1‑238 (MQ=255)                                                                                                                                     tccaCGCCCCGTGGGTTGCCCACTCCGCGTCCGGGGCGGCACCGCGCCCGGCCAGCGA‑TGGGGGG‑T‑‑TTGCCTCGCGTCCCCCCC‑CCCGTT‑CTTTAACCGCACCCGGGGGG‑ACA‑GCACCCTCCCCACTATTCTTC‑CCCAAGAT‑ATCGAGT‑TGCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCAtt                                                                                                   <  2:446427/236‑1 (MQ=255)                                                                                                                                                                       gCGGGCTACGCTGCCAGCCAGCAATTTGCGGG‑G‑‑TTT‑CTCCCGGTCCCCCC‑CCCCCGCCTTTAACC‑CACCCCCGCGG‑GCA‑GTAAACCCCCCCCCCCTCTTT‑TCCCCCAT‑ATCGCGA‑AGCCCGCCCCCCCCTAAACCCAAAATACCCCTGGTGGACCCCCCGGCCGGGCGGACCCCGGGCCGGTCCGTTCGACCCAAACGGCCCCCTTAAGCCTCCTCCCGCCTTTTTTTCCCGccccctc                                                                  >  2:2193159/1‑236 (MQ=255)                                                                                                                                                                               ccgtGCCCGCCAGGATGTGGGCGG‑T‑‑TTGTCTCTCGTCCCGCC‑‑CCCGTT‑CTTTAACCCCCCCCCGGGGGAGCA‑GCCCCCCCCCCCCGTTTTTTT‑CCCAAGAT‑ATCGAGA‑CGCCCGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGa                                                          <  2:954081/236‑1 (MQ=255)                                                                                                                                                                                                                       ccgccc‑CCCTTG‑CTTAAACCCTCCCCGGCCGG‑CCA‑G‑CCCCCCCCCCCCTTTCTTTCCCCCAGTT‑ATCGTGG‑TGCCTCCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGcc                     <  1:1347913/239‑1 (MQ=255)                                                                                                                                                                                                                         ccccccccGTT‑GTTTTACCCTACCCCGGGGG‑GCAGGTACCCCCCGCCGGA‑TCTTC‑CCCCAAATAAGCGAGG‑TCCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGcc                     <  1:1253058/234‑1 (MQ=255)                                                                                                                                                                                                                          ccg‑ccCGTG‑CTTTAACCGTACCCAGGCGG‑GCA‑GTAACCTCCGCCCGATTCTTT‑CCCCAGAT‑ATCGAGA‑TGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCtt                   >  2:730574/1‑238 (MQ=255)                                                                                                                                                                                                                              cccGTG‑CTTTAACCGTACCCAGGCGG‑GCA‑GTAACCTCCGCCCGATTCTTT‑CCCCAGAT‑ATCGAGA‑TGCCAGCCCGCCGTTCAGCCAAAAGTACCGTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAACAAACGGCCCCATTAAGCCTTCCGCCGACATTTTTTGCCGCCCCATATCCGCCGAAAAGTCGATCTGCCGACCATCCTGAGTGGCAACTGCCTTCg                 >  1:261889/1‑237 (MQ=255)                                                                                                                                                                                                                              cccGTG‑CTTTAACCGTACCCAGGCGG‑GCA‑GTAACCTCCGCCCGATTCTTT‑CCCCAGAT‑ATCGAGA‑TGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGc                >  2:898713/1‑238 (MQ=255)                                                                                                                                                                                                                                         aaCCGTACCCAGGCGG‑GCA‑GTAACCTCCGCCCGATTCTTT‑CCCCAGAT‑ATCGAGA‑TGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAg                                                      >  1:1290243/1‑190 (MQ=255)                                                                                                                                                                                                                                          cccGTCCCCGGGCGG‑TCA‑GCCACCCCCGCCCGTTTCTTC‑CCCCAGAT‑ACGGAGA‑GCCCGGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTgttttccgct    <  1:1340490/238‑10 (MQ=255)                                                                                                                                                                                                                                          aCCGTACCCAGGCGG‑GCA‑GTAACCTCCGCCCGATTCTTT‑CCCCAGAT‑ATCGAGA‑TGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGCCATTTTTTGCCGCTCCATATCCGCCCAAAAGTCGATCTGCGGACCCTCCTGAGGGGCAACTGCCTTCCCCTGGGCTTCag     >  2:2201389/1‑238 (MQ=255)                                                                                                                                                                                                                                          aCCGTACCCAGGCGG‑GCA‑GTAACCTCCGCCCGATTCTTT‑CCCCAGAT‑ATCGAGA‑TGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTCAGTGGCAACTGCCTTCGCCTGGGCTTCag     >  2:1458630/1‑238 (MQ=255)                                                                                                                                                                                                                                             gTACCCAGGCGG‑GCA‑GTAACCTCCGCCCGATTCTTT‑CCCCAGAT‑ATCGAGA‑TGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCagcta  >  2:184756/1‑238 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                    GCTGCCGGGTTTTGCTGGCATTAATATCGGTTTCCACACCTTCAACTGAAGAAGTAATCCCGTTCTGATATAGCTGGCGATCGGTCGCGATAATGGTGTTCTGCTCTTTTTCTATTTGCTGCAAGACCGTGTTTAACGCCGCCTGGGTTTGCCACTCCCAGTACAGGCGGGCTACGCTGCCAGCCAGCAATTGGCGGG‑T‑‑TTG‑CTCGCGTTCCGCCG‑CCCGTG‑CTTTAACCGTACCCAGGCGG‑GCA‑GTAACCTCCGCCCGATTCTTT‑CCCCAGAT‑ATCGAGA‑TGCCAGCCCGCCGTTAAGCCAAAAGTACCGTTGGTGTA‑CCACGGGCCGGTCGTACCTGCGGCCGGATCGTTCAGAGCAAACGGCCCCATTAAGCCTTCTGCCGACATTTTTTGCCGCTCCATATCCGCCGAAAAGTCGATCTGCGGACCATCCTGAGTGGCAACTGCCTTCGCCTGGGCTTCAGCTA  >  NC_007779/2231358‑2231825

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