Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	2,153,826	G→T	33.9%	intergenic (+43/+22)	Y75_RS10870 → / ← pphC	YegJ family protein/protein‑serine/threonine phosphatase PphC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	2,153,826	0	G	T	33.9%	24.6 / 9.5	42	intergenic (+43/+22)	Y75_RS10870/pphC	YegJ family protein/protein‑serine/threonine phosphatase PphC
Reads supporting (aligned to +/- strand):  ref base G (20/6);  new base T (10/4);  total (30/12)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

GCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAA‑TTGTTGTC‑TT‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATC‑ATCATCCGTACG‑CT‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATC  >  NC_007779/2153635‑2154054                                                                                                                                                                                                 |                                                                                                                                                                                                                                            gcaaGGAGTGAAAAAAAAAGAAAGAAGATATTTCAGAGTGGATGTATACGCGGGGGGATAAAATTTACGGGGGTTAACCCATTGATCCTTTATTGGTCACCTACCAAAAAAAAGAGGCCGACGAGTAAAGAGAAAAATTATTTCGCAAATCAGCCGTTATATTTATTTTTAAAACGGGCGGTTGAA‑TTGTTGTC‑TT‑AAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGcc                                                                                                                                                                                             <  1:1072758/235‑1 (MQ=255)  ccAGGAGGGAAAAAACAAAAAAGAAGATATTGCAGACGGGATGTACCCGGGCGGGAATAAAATTTCCGGGGTTTATACTATTGTTCTTTTATGGTTCACCTAACAAAAAGAAGAAGCCGAGGAGCTAAGAGAAAAATTATTTCGTTAATCAGCCTTTTTATTTATTGTTAAAACGGGCGGTTGAA‑TTGTTGTC‑TT‑AAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCa                                                                                                                                                                                            <  1:1864762/238‑1 (MQ=255)             aaaaCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGGGGGGAAAAAATTTAGGGCGGTTATACCTTTGATCCTTTATTGGTCACCTATCAAAAAAAAAAAGCCGACGAGCTAAGAACAAAATTAGTTCGCTAATCAGCTTTTATATTTTTTGTTTAAACGGGGGTTTGAA‑TTTTTGCC‑TT‑CAGCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTatc‑at                                                                                                                                                                               <  1:640671/239‑1 (MQ=255)                                                                 ttACGGCGGTTATACCATTGACCTTTTATTGGCCACCTACCAAAAAGAAGAAGCCGACGAGAAAAGAAAAAAATTATTCCGCTAATCACCCGTTATTTTTTTTGTTAAACCGGGCGTTTAAT‑TTGTTGTC‑TC‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATC‑ATCATCCGTACG‑CT‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTc                                                                                                                        <  1:567033/239‑1 (MQ=255)                                                                                   tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGTA‑TTTTTTTT‑TT‑CAACCGCCCTGAGAAAATCACCCTTCCCCCAACGCCACTGCCCACGTTTTATC‑ATCATCCCTACC‑CC‑CGTTTTCCCCCCGGGTGG‑ACAAAAACTTCCTTTACCAATCCGGTAAAAAAAAGAg‑ttg                                                                                                     >  2:1293976/1‑239 (MQ=255)                                                                                   tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAA‑TTTTTGTC‑TT‑CAACCGCCCTGAGAAACTCACTCTGCCCACAGCGCCCGTGCCAACGTTTTTTT‑ATCATCCGTACG‑CC‑CGTTTACCGCCGGGCTGG‑GCAAAAACTGCTTTAACAATTCAGGTAATAAATCGAg‑ttg                                                                                                     >  1:704017/1‑239 (MQ=255)                                                                                   tGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAA‑TTGTTTCC‑TT‑CAACCGCCATGAGAAAATCACTCTGCCCCCCGCGCCCGGGCCCAAGTTTTTTC‑ATCATCCCTACG‑CT‑CGTTTACCGCCCGGCTGG‑AAAAAAAATGCTTTAACAAATCAGGGAATAAAACGAg‑ttg                                                                                                     >  1:584397/1‑239 (MQ=255)                                                                                       ccTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAA‑TTTTTGTC‑TT‑CAACCGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATC‑ATCATCCGTACG‑CT‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAg‑ttgtt                                                                                                   >  1:605047/1‑237 (MQ=255)                                                                                          ttATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAA‑TTGTTGCC‑TT‑CAACCGCCATGAGAAAATCACTCTGCCCCAAGCGCCAGTGCCAACGTTTTATC‑ATCATCCCTACG‑CT‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGc                                                                                              >  1:1017739/1‑239 (MQ=255)                                                                                          ttATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGATTAA‑TTTTTTTC‑TT‑TAACCGCCCTGAGAAAATAACTCTGCCCCCAACGCCAGTGCCAACGTTTTATC‑ATCATCCGTACG‑CT‑CGGTTACCTCCGGGCTGG‑GCAAAACATTCTTTAACAAAATAGGGAATAAAACGAG‑GTTTTCCTg                                                                                               >  2:1024404/1‑238 (MQ=255)                                                                                                gTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGTA‑TTTTTTTT‑TT‑CAACCGCCCTTAGAAACTCACTCCTCCCACAGCGCCCGTGCCAACATTTTATC‑ATCATCCCTACC‑CT‑CGTTTACCCCCCGGCTGG‑GCAAAAAATGCTTTTACAAATCAGGGAATAAAAAAAG‑GTGTTCCTGCGttttc                                                                                        >  2:706687/1‑236 (MQ=255)                                                                                                     cTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAA‑TTTTTTTT‑TT‑TAACCGCCCTTAGAAACTAACTCTGCCCACAACGGCAGTGGCACCGTTTTATA‑CTCATCCCTACG‑CC‑CGTTTTCCGCCGGGCTGG‑GCAAAAACTCCTTTTACCATTCAGGGAATAAAACGAG‑TTGTTCCTTCGTTGTTgtcg                                                                                    >  2:212369/1‑238 (MQ=255)                                                                                                     cTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAA‑TTTTTTTT‑TT‑TAACCGCCATGAGAAACAAACTCCTCCCACCACCCCCGTGCCAACGGTTTTTC‑ATCATCCGTACG‑CC‑CGTTTTCCGCCGGGCCGG‑GCAAAAACTGCTTTTACAATTCCGGTAATAAATCGAG‑GTGTTTCTTCCGTTCTgccg                                                                                    >  2:624591/1‑238 (MQ=255)                                                                                                        tCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAA‑TTGTTGCC‑TT‑CAACCGCCATGAGAAACTCACCCTCCCCCCAGCGCCCGTGGCAACCGTTTATC‑ATCATTCGTACGCCC‑CGTTTACC‑CCCGGCCGG‑GCAAAAAACGCTTTTACAATTTCGGTAATAAAACGAGTTTTTTCCTCCG‑TGTTGCCGCCAg                                                                                >  2:3581/1‑239 (MQ=255)                                                                                                                  gaagCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAA‑TTGTTGCC‑TT‑CAACCGCCCTGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGGTTTATC‑ATCATCCGTACG‑CT‑CGTTTACCCCCGGGCCGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGCGTTGCTGCCGCCAGTCCATTGaaa                                                                      >  1:837137/1‑239 (MQ=255)                                                                                                                      ccagcgagaaaaGAGCAAAATTATTTCATTAACCAGCCTTTATATTTATTGTTAAACCGGGCGGTTGAT‑TTTTTGTT‑TC‑CACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATC‑ATCATCCGTACG‑CT‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGCGTTGCTGCCGCCAGTCCATTGaaaaa                                                                    <  1:767858/227‑1 (MQ=255)                                                                                                                         acgaGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAA‑TTGTTTCC‑TT‑CCACCGCCCTGAGAAAATCACTCTTCCCCAAGCGCCCGTGCCAAAGTTTTATC‑AACATCCGTACG‑CT‑CGTTTACCGCCCGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGCGTTGCTTCCGCCAGCCCATTGAAAAACgggg                                                               >  1:659079/1‑239 (MQ=255)                                                                                                                                aaGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGGTGAA‑TTTTTGCC‑TT‑CAACCGCCCTGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATC‑ATCATCCGCACG‑CT‑CGTTTACCGCCGGGCTGG‑ACAAAAACAGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGCGTTTCTTCCGCCCGTCCATTGAAAAACGGGGGaaagaa                                                        >  2:1199566/1‑239 (MQ=255)                                                                                                                                       aaaTTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTTAA‑TTGTTTTT‑TT‑CCACCGCCCTGAGAAACTAAATCTTGCCACAAGGCCAGTGGCCAAGGTTTATC‑ATCATTCGGACG‑CC‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTAAGGTAATAAATAGAG‑GTGTTTCTGCGTTGCTGCCGGCCGTCCATTTAAAAAAGGGGAAAAGAAAAGgaca                                                 >  2:794281/1‑239 (MQ=255)                                                                                                                                           tgtttctcaaaTCAGCTTTTATTTTTTTTGTTTAACCTGGCGTTTAAA‑TTGTTGCC‑TT‑CACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATC‑ATCATCCGTACG‑CT‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATgag                                             <  1:279746/230‑1 (MQ=255)                                                                                                                                           tgtttcctaaaTCAGCGTTTATTTTTATT‑TTTAAACCGGGTGTTTAATTTGTTCCC‑TA‑CACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATC‑ATCATCCGTACG‑CT‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATgag                                             <  1:1186026/230‑1 (MQ=255)                                                                                                                                             gTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAA‑TTGTTTCC‑TT‑CAACCGCCCTGAGAAACTCACCCTGCCCCCAGCGCCCGTGCCAACGTTTTTTC‑ATCATCCGTACG‑CC‑CGTTTACCGCCGGGCTGG‑ACAAAAAATGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAAAGGGGGAAAAAAAAGCACATgagga                                           >  2:176897/1‑239 (MQ=255)                                                                                                                                             gTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAA‑TTGTTGTT‑TT‑CAACCCCCACGAGAAAATCACTCTTCCCACAGCCCCAGTGCCCACGTTTTTTC‑ATCATCCGTACG‑CC‑CGTTTTCCGCCCGGCTGG‑GCAAAAACTGCTTTAACAATTCAAGTAATAAATCGAG‑GTGTTCCTGCGTTGCTGCCGCCAGTCCCTTGAAAAAAGGGGGAAAGAAAAGCACATgggga                                           >  2:230276/1‑239 (MQ=255)                                                                                                                                                cGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGGTTCA‑CTATTGTC‑TT‑CTAACCCCCTTAGAAAATCACTCTGTCCCCCTCGCCCGGTCCCACGTTTTGTCTTTTATCC‑ACCT‑CTCCCTTTCCCCCCGGGC‑GG‑GCAAAAAATGCTTTAACCAATCAGGGAATAAATCCAG‑TTGTTCCTGGGTTGCTGCCCCCCGTCCATTGAAAAACGGGGGAAAGAAAGGcacaagtggctcg                                        >  2:1790611/1‑230 (MQ=255)                                                                                                                                                  cTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAA‑TTGTTGCC‑TT‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATC‑ATCATCCGTACG‑CT‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAAtt                                      >  1:1170994/1‑239 (MQ=255)                                                                                                                                                         agaagattttttatt‑TTTTAACCGGGGGTTGAA‑TTTTTTTC‑TTCCACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATC‑ATCATCCGTACG‑CT‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGc                               <  1:498680/230‑1 (MQ=255)                                                                                                                                                          cttttATTTTTTTTGTTAAAACGGGGGGTTAAT‑TTTTTCCC‑TT‑AACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATC‑ATCATCCGTACG‑CT‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCa                              <  1:1199566/236‑1 (MQ=255)                                                                                                                                                          cttttATATTTTTTG‑TAAAACCGGGGGTTTAATTTGTTTTC‑TC‑CACCCGCCATGAGAAACTTACTCTGCCCACAGCGCCAGTGCCAACGTTTTATC‑ATCATCCGTACG‑CT‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCa                              <  1:1463427/236‑1 (MQ=255)                                                                                                                                                             ttATATTTATTGTTTAAACTGGCGGTTGAA‑TTTTTGTC‑TT‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCCCCAGGGCCAACGTTTTATC‑ATCATCCGTACG‑CC‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGGAATAAATCGAG‑TTGTTCCTGCGGTGCCGCCGCCAGTCCATTGAAAAAAGGGGGAAAGAAAGGCACATGAGGAGAATTATCCAGCatat                           >  1:1772377/1‑239 (MQ=255)                                                                                                                                                                     attGTTTAAACTGGCGGTTGAA‑TTGTTGTC‑TT‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGGTTTATC‑ATCATCCGTACG‑CC‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAAACGAG‑TTGTTCCTGCGTTGCTTCCCCCAGTCCATTGAAAAACCGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGcc                   >  1:1056024/1‑239 (MQ=255)                                                                                                                                                                     attGTTTAAACTGGCGGTTGAA‑ATGTTGTC‑TT‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTat                                                                                                                                                                                   >  1:135027/1‑84 (MQ=255)                                                                                                                                                                       tGTTTAAAGCGGCGTTAAAT‑TTGTTTTA‑TA‑CACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATC‑ATCATCCGTACG‑CT‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAg                 <  1:1539779/239‑1 (MQ=255)                                                                                                                                                                       tGTTTAAACTGGCGGTTGAA‑TTGTTGTC‑TT‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATC‑ATCATCCGTACG‑CT‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGGAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAg                 >  1:239928/1‑239 (MQ=255)                                                                                                                                                                       tGTTTAAACTGGCGGTTGAA‑TTGTTGTC‑TT‑CAACAGCCATGAGAAAATCACTCTGCCCACAGCGCCAGTGCCAACGTTTTTTC‑ATCATCCGTACG‑CC‑CGTTTAACCCCCGGCTGG‑ACAAAAAATTCTTTAACAAATCAGGTAATAAATCGAG‑TTGTTGCCGCGTTGCTGCCGCCAGTCCACTTAAAAACGGGGGAAAGAAAGGCACATGAGGAGAATTATACAGCATATTTAAAGCCAg                 >  2:1170363/1‑239 (MQ=255)                                                                                                                                                                       tGTTTAAACTGGCGGTTGAA‑ATGTTGTC‑TT‑CAACCGCCATGAGAAACTCACTCTGCCCCCAGCGCCAGTGCCAACGTTTTATC‑CTCATCCGTACG‑CT‑CGTTTACCGCCCGGGTGG‑ACAAAAACTGCTTTAACAATTCAGGGAATAAATCGAG‑TTGTTCCTGCGTTTCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACACCAg                 >  1:1624396/1‑239 (MQ=255)                                                                                                                                                                          taaaccgggCGGTT‑AA‑TTTTTTTCTTT‑CACCCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATC‑ATCATCCGTACG‑CT‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCg              <  1:230276/232‑1 (MQ=255)                                                                                                                                                                           tAAACTGGCGGTTGAA‑TTGTTTTC‑TT‑CAACCGCCATGAGAAACTCACTCTGCCCCCAGCGCCAGTGCCCACGTTTTTTC‑ATCATCCGTACG‑CT‑CGTTTACCCCCCGGCTGG‑ACAAAAAATTCTTTAACAAATCAGGGAATAAAACGAG‑TTGTTTCTGCGTTGCTGCCGCCAGTCCCTTGAAAAACGGGGGAAAGAAAAGCACATGAGGGGAAATATCCAGCATATTTAACGCCAGTCg              >  2:648304/1‑238 (MQ=255)                                                                                                                                                                           tAAACTGGCGGTTGAA‑TTGTTGCC‑TT‑CAACCGCCCTGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATC‑ATCATCCGTACG‑CT‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTc               >  1:1438901/1‑237 (MQ=255)                                                                                                                                                                            aaaCTGGCGGTTGAA‑TTTTTGCC‑TT‑CAACCGCCCTTAGAAACTCACTCTGCCCACAGCGCCCGTGCCAACGTTTTATC‑ATCATCCGTACG‑CC‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGCGTTGCTGCCGCCCGTCCATTGAAAAACGGGGGAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAAAGCCAGTc               >  1:1432622/1‑236 (MQ=255)                                                                                                                                                                                   cGGTTGAA‑TTGTTGCC‑TT‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATC‑ATCATCCGTACG‑CT‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATAcc     >  1:864418/1‑239 (MQ=255)                                                                                                                                                                                   cGGTTGAA‑TTGTTGCC‑TT‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATC‑ATCATCCGTACG‑CT‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGGAATAAATCGAG‑TTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAAAGGGGGAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATAcc     >  1:1097688/1‑239 (MQ=255)                                                                                                                                                                                      ttGAA‑‑TGTTGTC‑TT‑CAACCCCCCAGAGAAAACCACTCTTCCCACAGCGCCAGTGCCAACGTTTTATT‑ATCATCCTTACT‑CC‑CGTTTACCCCCCGGCTGGCAAAAAAACT‑CTTTAACAAATCAGGTAATAAATCGAG‑TTGTTCCTGCGGTGCTGCCCCCAGTCCCTTGAACAAAGGGGGAAAGAAAGGCACATTAGGGAAATTAACCAACATTATTAACGCCAGATGCACGATACCacc  >  2:1459320/1‑236 (MQ=255)                                                                                                                                                                                                 gtC‑TT‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATC‑ATCATCCGTACG‑CT‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGc                               >  2:159307/1‑200 (MQ=255)                                                                                                                                                                                                 gtC‑TT‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATC‑ATCATCCGTACG‑AT‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAgg                                                     >  2:738550/1‑178 (MQ=255)                                                                                                                                                                                                 gtC‑TC‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATC‑ATCATCCGTACG‑CT‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAgg                                                     <  1:738550/178‑1 (MQ=255)                                                                                                                                                                                                 |                                                                                                                                                                                                                                            GCCAGGAGTGGAAAATCAAGAAAGAAGATATTTCAGACTGGATGTATACGCGCGGGGATAAAATTTACGGCGGTTATACCATTGATCCTTTATTGGTCACCTATCCAAAAGAAGAAGCCGACGAGCTAAGAGCAAAATTAGTTCGCTAATCAGCCGTTATATTTATTGTTTAAACTGGCGGTTGAA‑TTGTTGTC‑TT‑CAACCGCCATGAGAAACTCACTCTGCCCACAGCGCCAGTGCCAACGTTTTATC‑ATCATCCGTACG‑CT‑CGTTTACCGCCGGGCTGG‑ACAAAAACTGCTTTAACAATTCAGGTAATAAATCGAG‑TTGTTCCTGCGTTGCTGCCGCCAGTCCATTGAAAAACGGGGTAAAGAAAGGCACATGAGGAGAATTATCCAGCATATTTAACGCCAGTCGCTGGATACCATC  >  NC_007779/2153635‑2154054

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