Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	4,021,098	T→C	17.2%	G42G (GGA→GGG)	Y75_RS19860 ←	rhs element protein RhsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	4,021,098	0	T	C	17.2%	175.2 / 2.1	78	G42G (GGA→GGG)	Y75_RS19860	rhs element protein RhsB
Reads supporting (aligned to +/- strand):  ref base T (47/14);  new base C (2/13);  total (49/29)
Fisher's exact test for biased strand distribution p-value = 8.66e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCTGTAAGCGGATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACT‑GCTGTAGGTGCGGGAGAGGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAATATATTTAGAA  >  NC_007779/4020906‑4021328                                                                                                                                                                                                     |                                                                                                                                                                                                                                          cAGCTGTAAGCGGATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACT‑GCTGTAGGTGCGGGAGAGGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATAc                                                                                                                                                                                            >  1:976828/1‑239 (MQ=11)     ggTAAGCGGAGATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACT‑GCTGTAGGTGCGGGAGAGGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TAc                                                                                                                                                                                       <  2:988096/238‑1 (MQ=17)        aggaggATTTCCGCAGGCATTGTCTAGGCGGGGCCGAGGCTGCCGACCGGCGCGGGCTTTTTTGTCCGGTAACT‑GCTGTAGGTGCGGGAGAGGATGAAGGCGT‑GTGGGCCGTGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGGGCGGGGGATTGACCGGAGGGCTGGA‑CGTCACACCGCCGGGGCCCACTGAGC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TAccc                                                                                                                                                                                     <  2:1507690/234‑1 (MQ=11)             gtgttcggctggCCTTTTCCAGCCGGGGCCGAGCCTCCCCCCCGGCGCGGGCGTTTTTGTCCGGTAACT‑GCTGTAGGTGCGGGAGAGGATGAACGGCC‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACCCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCt                                                                                                                                                                                <  2:797525/228‑1 (MQ=11)             gATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACT‑GCTGTAGGTGCGGGAGAGGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCt                                                                                                                                                                                <  1:1443165/238‑1 (MQ=32)              atatCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACT‑GCTGTAGGTGCGGGAGAGGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGa                                                                                                                                                                              >  1:386595/1‑239 (MQ=32)              atatCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACT‑GCTGTAGGTGCGGGAGAGGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGa                                                                                                                                                                              >  1:1942137/1‑239 (MQ=32)              atatCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACT‑GCTGTAGGTGCGGGAGAGGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGa                                                                                                                                                                              >  1:1881215/1‑239 (MQ=32)              atatCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACT‑GCTGTAGGTGCGGGAGAGGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGa                                                                                                                                                                              >  1:760807/1‑239 (MQ=32)              atatCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACT‑GCTGTAGGTGCGGGAGAGGACGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGC‑CGGCACCCCGCCGGGGCACACCGGAC‑CGGCCAAGCCCGTCGGGGCACCAAAACC‑TACCCCgggcga                                                                                                                                                                              >  2:41556/1‑235 (MQ=11)                     cAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACT‑GCTGTAGGTGCGGGAGAGGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑a                                                                                                                                                                      >  1:51074/1‑239 (MQ=32)                                  gCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACT‑GCTGTAGGTGCGGGAGAGGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCg                                                                                                                                                          >  1:851520/1‑238 (MQ=32)                                     ggggCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACT‑GCTGTAGGTGCGGGAGAGGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGGGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATCCTACCGCCAt                                                                                                                                                      >  1:1833760/1‑239 (MQ=17)                                       ggCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACT‑GCTGTAGGTGCGGGAGAGGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGG‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATa                                                                                                                                                     >  2:1405012/1‑238 (MQ=17)                                            ggtgcccccACCGGCGCGGGCGTTTTTGTCCGGTAACT‑GCTGTAGGTGCGGGAGAGGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCCTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACCCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAACACG‑TACCCCGGGTGAACCATG‑CACAAGTCTACAGACCTACTCCg                                                                                                                                                >  2:2135377/6‑238 (MQ=11)                                            gggttCCCCGCGGGGGGGGGCTTTTGTGTCGGGTAACT‑GCTGGGGGTGGGGGGGAGGAT‑ACAGGGA‑GGGGGGGGGGGGGGGGGGGGGCGGTTTCCCCGGGGAGGGCCTTTGGCCGCGGGCGGGGGATGGACCGGATGGCGGGC‑GCTCCCTCCCCGGGGCCACCCCGAAC‑AGGCCCCG‑CGGGGGGGGCCCCAGTTCGTTCCCCCGCCTGACCCCTG‑AACAATGCTACCGCCATACTGCg                                                                                                                                                <  2:584676/234‑1 (MQ=2)                                                 ccccACCGGCGCGGGCGTTTTTGTCCGGTAACT‑GCTGTAGGTGCGGGAGAGGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATg                                                                                                                                           >  1:2012004/1‑238 (MQ=32)                                                            cgGGCGTTTTTGTCCGGTAACT‑GCTGTAGGTGCGGGAGAGGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGAcgc                                                                                                                               >  1:1262623/1‑239 (MQ=33)                                                            cgGGCGTTTTTGTCCGGTAACT‑GCTGTAGGTGCGGGAGAGGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAAGGCTACCGCCATACTGCGTCATGTCACCCTGAcg                                                                                                                                >  2:1599385/1‑238 (MQ=18)                                                              ggCGTTTTTGTCCGGTAACT‑GCTGTAGGTGCGGGAGAGGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑CACAATGCTACCGCCATACTGCGTCATGTCACCCTGAcgcg                                                                                                                              >  1:1785530/1‑238 (MQ=18)                                                                   ttttGTCCCGTAACTGGGTGTGGGGGGGAGAGGGGATG‑ACGGCAGGGGGGCCGGGAGGGGGGTGGTGGGTTGCCCGGGGAGGGCCCTTT‑GCACGGGGCGGGGGTTGGCCCGGTGGGCCGGC‑CGCCCCTCCGCCGGGGCACCCCAACC‑GGGCCCCGCCGGGGGGGGCACCAATCCG‑TCCCCCGGCGGACCCCTG‑AAAAAGGCTACCGCCATACTGCGTCATGTCACCCTGACGcgccgc                                                                                                                          <  2:851520/238‑1 (MQ=255)                                                                        tCCGGTAACT‑GCTGTAGGTGCGGGAGAGGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGttt                                                                                                                    >  1:1079784/1‑238 (MQ=35)                                                                                        gggTGCGGGAGAGGATGAACGGCA‑GGGGGCCGGGCGGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGAGGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGACCCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTaaaa                                                                                                    <  1:1887384/238‑1 (MQ=255)                                                                                            gCGGGAGAGGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCt                                                                                                >  1:2049250/1‑239 (MQ=40)                                                                                            gCGGGAGAGGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCCCATCGGAAAATCCt                                                                                                >  2:1772270/1‑239 (MQ=255)                                                                                              gggAGAGGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑AGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCAAACGCCATAATGCGTCATGTCACCCTGCCGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCttt                                                                                              >  2:1336790/1‑239 (MQ=255)                                                                                                   gggagaaaCGGGA‑GCGGGCGGGGCGGGGTGGTGTCGTTCGCCCCGGGCAGGCCCTTT‑GCCCGGCGGAGGGGATTGCCCGGGGGGCCGGG‑CCACACCCCGCCGGGCCCCCCCGAAC‑GGGCCCCGCGGGGGGGGGCACCAATACG‑TCCCCCGGCTGCCCCCGA‑AAAAATTCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCATAAACTCCTTTTCta                                                                                          <  2:2012004/232‑1 (MQ=9)                                                                                                   agGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCtat                                                                                         >  2:1156910/1‑239 (MQ=255)                                                                                                     gATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCt                                                                                                                                                                                >  2:6624/1‑151 (MQ=33)                                                                                                     gATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCt                                                                                                                                                                                <  1:6624/151‑1 (MQ=33)                                                                                                          agggcg‑gcGGGCCGGGAAGGGCGATGTCGGTTTCACCGGGAAGGCCCTTT‑GCACCGAGCAGGGGATTGACCGGAGGGCCGGA‑CGTCACTCCGCCGGGGCACCCCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TCCCCCGGCTGACCCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTAtttt                                                                                   <  1:1324642/236‑1 (MQ=255)                                                                                                          aCGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCATCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTa                                                                                  >  2:411537/1‑239 (MQ=255)                                                                                                          aCGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCt                                                                                                                                                                                >  2:582357/1‑146 (MQ=34)                                                                                                          aCGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCt                                                                                                                                                                                <  1:582357/146‑1 (MQ=34)                                                                                                          aCGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTa                                                                                  >  1:1115791/1‑239 (MQ=255)                                                                                                            ggCA‑GGGGGCGGGGCGGGGAGTTGGCGGTTGCCCCGGGAAGGCCTTTT‑GCCCCGAGCGGGGGATTGACCGGAGGGCGGGA‑CGTCCCCCCGCCGGGGCACCCCGAAC‑AGGCCACGCCGGTGGGGCCACCAATCCG‑TCCCCCGGCTGACCCCTG‑AACAATTCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTaa                                                                                 <  2:31620/238‑1 (MQ=14)                                                                                                              cggcgggggggggcGGGGGGATGCGGCTTTCACGGGGAGGGCCCTTG‑GCCCCGGGCGGGGGATTGACCGGTGGGCCGGA‑CGCCACCCCCCGGGGCCACCCAGAAA‑GGGCCACGCCGGGGGGGGCACAAATACG‑TCCCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAAt                                                                                <  2:1317383/228‑1 (MQ=17)                                                                                                                 gCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACa                                                                            >  1:644991/1‑239 (MQ=255)                                                                                                                         ggcAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATa                                                                                                                                                     >  1:767167/1‑158 (MQ=33)                                                                                                                              ggCGATGTCGGTTTCACTGGGAAGGACCTTT‑GCACCGAGCCGGGGATTGACCGGATGGCCGGA‑CGACACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATCCG‑TACCCCGGCTGAACACTG‑GACAATGCTACCGCCGTACTGCGGCATGTCACCCTGACGCGCCGCCGGTTTCACGCTCACCGACAAGTCCTTTTCTATTATTTTAATGACAGTATACTGATACg                                                               >  2:1080993/1‑239 (MQ=14)                                                                                                                               gCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACACTGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCta                                                                                          >  1:576480/1‑211 (MQ=255)                                                                                                                                   tGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATg                                                                                                                                           >  1:1342794/1‑158 (MQ=33)                                                                                                                                   tGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGaaaaa                                                          >  2:154272/1‑239 (MQ=255)                                                                                                                                    gTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAAt                                                                                >  1:148606/1‑216 (MQ=255)                                                                                                                                        gTTTCCCGGGGAAGGCCCTTT‑GCACCGAGCAGGGGATTGACGGGAGGGCCGGA‑CGTCACTCCGCCGGGGCACCCCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TCCCCCGGCTGACCCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCga                                                      <  1:97074/238‑1 (MQ=18)                                                                                                                                                           ttggcgcggggCAGGGACTGGACGGGAGGGCCGGC‑CGCCACTCCGCCGGGGCACACCGACC‑AGGCCCCGCCGGTGGGGCCCCCAATCCT‑TCCCCCGGCTGAACCCGA‑AACAATCTTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTcaca                                   <  2:897442/229‑1 (MQ=21)                                                                                                                                                           tt‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCt                                                                                                                                                                                >  2:828288/1‑98 (MQ=37)                                                                                                                                                           tt‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCt                                                                                                                                                                                <  1:828288/98‑1 (MQ=37)                                                                                                                                                           tt‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACaa                                  >  1:1364090/1‑239 (MQ=255)                                                                                                                                                            t‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACa                                                                                                                                                                   <  1:2084824/109‑1 (MQ=35)                                                                                                                                                            t‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACa                                                                                                                                                                   >  2:2084824/1‑109 (MQ=35)                                                                                                                                                            t‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTAt                                                                        >  1:1287346/1‑200 (MQ=255)                                                                                                                                                            t‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATAc                                                                >  1:761734/1‑208 (MQ=255)                                                                                                                                                            t‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAAt                                 >  1:899050/1‑239 (MQ=255)                                                                                                                                                                   gAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTg                           >  2:1979143/1‑239 (MQ=255)                                                                                                                                                                    aGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGa                          >  2:2030436/1‑239 (MQ=255)                                                                                                                                                                    aGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGa                          >  1:1041224/1‑239 (MQ=255)                                                                                                                                                                    aGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGGTTTCCGCTCATCGTCAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGCGGCCACAATGCATTGa                          >  2:174199/1‑239 (MQ=255)                                                                                                                                                                    aGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGATTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGa                          >  2:830649/1‑239 (MQ=255)                                                                                                                                                                     gcggggaATGGACC‑GATGGGCGGAGCGCCACCCCCCCGGGGCA‑ACCGGACAAGGCCACGCCGGTGGGGCCCCCAATCCG‑TCCCCCGGCTGAACCCGG‑AACAATGCTCCCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGa                          <  2:837604/235‑1 (MQ=21)                                                                                                                                                                     gCAGGGGATTGACCGGAGGGGCGGA‑GGCCACCCGGCGGGGGCACCCCGACC‑GGGCCACGCGGGTGGGGGACCAAATCCG‑TCCCCCGGC‑GAACCCCGAAACAATTCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGa                          <  2:1360484/238‑1 (MQ=21)                                                                                                                                                                      cAGGGGATTGACCGGATGGCCGGG‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCt                                                                                                                                                                                <  1:181114/88‑1 (MQ=255)                                                                                                                                                                      cAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCt                                                                                                                                                                                >  2:181114/1‑88 (MQ=38)                                                                                                                                                                         gggATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGaa                         >  1:1337848/1‑235 (MQ=255)                                                                                                                                                                            aTTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGAttttt                  >  2:1058150/1‑239 (MQ=255)                                                                                                                                                                             tgctccGGGTGCCCGCA‑CGTCCCCCCGCCGGGCCACACAAAAC‑AGGCCACGCCGGTGGGGGCACCAATACT‑TCCCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCTTCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTa                 <  1:1772270/235‑1 (MQ=25)                                                                                                                                                                              tGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAg                >  1:1835315/1‑239 (MQ=255)                                                                                                                                                                                 cgggaggcccGGA‑CTTCCCCCCGCCGGGGCACCCAGAAC‑AGCCCCCGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGaa             <  1:348779/231‑1 (MQ=25)                                                                                                                                                                                 ccGGAGGGCCGGG‑CGCCACCCCGCCGGGGCCGCCCGAAC‑GGGCCACGCGGGGGGGGACACCATTACG‑TCCCCCGGCTGACCCCTG‑AACAATGTTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGaa             <  2:1906525/239‑1 (MQ=25)                                                                                                                                                                                  gggAGGGCCGCA‑CGTCGCCCCGCGGGGGCACCCCGAAC‑AGGCCACGCCGGTGGGGGCCCCAATACG‑TCCCCCGGCTAAACCATG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAat            <  1:2135377/238‑1 (MQ=25)                                                                                                                                                                                  cggtgggCGGGA‑CGTCCCCCCGCGGGGCCCCCCCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TCCCCCGGCTAAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAat            <  1:275892/234‑1 (MQ=25)                                                                                                                                                                                  cGGATGGCCGGA‑CGTCCCCCCGCCGGGGCACCCCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TCCCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAat            <  1:649132/239‑1 (MQ=255)                                                                                                                                                                                   ggATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAata           >  2:778700/1‑239 (MQ=255)                                                                                                                                                                                     aTGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACtt                                            >  1:745969/1‑204 (MQ=255)                                                                                                                                                                                     aTGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAatata         >  1:216489/1‑239 (MQ=255)                                                                                                                                                                                        cccGGG‑CCTCACCCCGCCGGGGCACCCCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TCCCCCGGCTGACCCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAATATAtt       <  1:814005/237‑1 (MQ=255)                                                                                                                                                                                             a‑CGTCACGCCGCCGGGCCACCCCGAAC‑AGCCCACGCCGGTGGGGGCACCAGTCCG‑TCCCCCGGCTGAACCCGG‑ACCAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTCTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAATATATTTAGaa  <  1:1614678/238‑1 (MQ=37)                                                                                                                                                                                              ccGTCACTCCCCCGGGGCCCACCGACC‑AGGCCACGCCGGTGGGGGCCCAAATCCG‑TCCCCCGGCTAAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAATATATTTAGaa  <  1:461304/237‑1 (MQ=37)                                                                                                                                                                                                     |                                                                                                                                                                                                                                          CAGCTGTAAGCGGATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACT‑GCTGTAGGTGCGGGAGAGGATGAACGGCA‑GCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACCGAGCAGGGGATTGACCGGATGGCCGGA‑CGTCACTCCGCCGGGGCACACCGAAC‑AGGCCACGCCGGTGGGGGCACCAATACG‑TACCCCGGCTGAACCCTG‑AACAATGCTACCGCCATACTGCGTCATGTCACCCTGACGCGCCGCCGGTTTTCCGCTCATCGTAAAATCCTTTTCTATTATTTTAATGACAGTATCTTCATACGAAAAAGCGAAGAAAAACTTGTGGTCACAATGAATTGAAGATTTTTAGGAATATATTTAGAA  >  NC_007779/4020906‑4021328

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