Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,388,044	T→A	10.0%	D149V (GAC→GTC)	PP_1209 ←	cold‑shock protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,388,044	0	T	A	10.0%	212.0 / 13.1	80	D149V (GAC→GTC)	PP_1209	cold‑shock protein
Reads supporting (aligned to +/- strand):  ref base T (35/37);  new base A (4/4);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.02e-01

TCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGC  >  NC_002947/1387896‑1388186                                                                                                                                                         |                                                                                                                                                tCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGt                                                                                                                                                 >  7:307619/1‑149 (MQ=255)   aaCGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGtct                                                                                                                                               >  5:214935/1‑149 (MQ=255)     cggcggTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCg                                                                                                                                             >  5:369731/1‑149 (MQ=255)      ggcggTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCgg                                                                                                                                            <  3:252556/149‑1 (MQ=255)          gTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGgaga                                                                                                                                        <  2:59055/149‑1 (MQ=255)           tttACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCTCCCGAGACTCGGGCGAt                                                                                                                                       >  1:349417/1‑149 (MQ=255)           tttACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGAt                                                                                                                                       <  5:328152/149‑1 (MQ=255)           gttACTTCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCGCCTAACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCggg                                                                                                                                           <  1:117294/148‑1 (MQ=255)             tACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCTCCCGAGACTCGGGCGATg                                                                                                                                      <  2:349417/148‑1 (MQ=255)             tACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATg                                                                                                                                      <  1:151845/148‑1 (MQ=255)               cTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCTCCCGAGACTCGGGCGATGaca                                                                                                                                   <  8:183074/149‑3 (MQ=255)                 gCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAAc                                                                                                                                  >  3:133124/1‑148 (MQ=255)                 gCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAAc                                                                                                                                  >  4:36745/1‑148 (MQ=255)                  cAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTACGGAGATGAAAcc                                                                                                                                 >  6:189763/1‑148 (MQ=255)                         gTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCtt                                                                                                                         >  8:244832/1‑149 (MQ=255)                          tCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCT‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCttt                                                                                                                        <  7:354106/149‑1 (MQ=255)                            ttCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCg                                                                                                                                             >  8:161087/1‑126 (MQ=255)                            ttCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCg                                                                                                                                             <  7:161087/126‑1 (MQ=255)                                gCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCttt                                                                                                                        >  3:333442/1‑143 (MQ=255)                                gCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCttt                                                                                                                        <  4:333442/143‑1 (MQ=255)                                 ccTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACgtg                                                                                                                 <  8:307619/149‑1 (MQ=255)                                  cTGCAGGCCCTTGTCTCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGTCCCTCACCGCGTATGGCCCTATAGTTCACGAAGATGTCATCTCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACgtgt                                                                                                                <  5:189763/149‑1 (MQ=255)                                  cTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACgtgt                                                                                                                >  4:96970/1‑149 (MQ=255)                                   tGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGtt                                                                                                               >  2:215315/1‑149 (MQ=255)                                    gCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGtt                                                                                                               >  8:69707/1‑148 (MQ=255)                                        gCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAAc                                                                                                           >  6:40997/1‑148 (MQ=255)                                        gCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAAc                                                                                                           >  7:267449/1‑148 (MQ=255)                                               tCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGa                                                                                                    <  2:246585/148‑1 (MQ=255)                                               tCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGAc                                                                                                   <  7:244832/149‑1 (MQ=255)                                                cgcgGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACgtgt                                                                                                                <  4:143357/135‑1 (MQ=255)                                                cgcgGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACgtgt                                                                                                                >  3:143357/1‑135 (MQ=255)                                                cgcgGTGCATCACCGAAAACTCCACGCGCTGCCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTAACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCACGAGTCTCGGGAGATGAAACCAAAGCCATTCGACGTGTTGAACCATTTGACg                                                                                                  >  5:377792/1‑149 (MQ=255)                                                 gcgGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCttt                                                                                                                        >  3:40106/1‑126 (MQ=255)                                                 gcgGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCttt                                                                                                                        <  4:40106/126‑1 (MQ=255)                                                   ggTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑g                                                                                              >  7:227186/1‑149 (MQ=255)                                                       cATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCg                                                                                           <  8:227186/148‑1 (MQ=255)                                                           aCCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTAt                                                                                       <  4:175033/148‑1 (MQ=255)                                                           aCCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTAt                                                                                       <  6:369731/148‑1 (MQ=255)                                                            ccGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTc                                                                                                                       <  8:123488/116‑1 (MQ=255)                                                            ccGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTc                                                                                                                       >  7:123488/1‑116 (MQ=255)                                                                   cTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGc                                                                              >  4:109128/1‑149 (MQ=255)                                                                    tCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGcc                                                                             >  3:156472/1‑149 (MQ=255)                                                                      cACGCTCGGGGCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTTACCGCGGATGGACCGAAAGTGCACGAAGATGTGATCGCCCGAGTCTCGGGAGATGAAACAAAAGCCTTTCGACGTGATGAACCATTTGACGGT‑GCCGGTATCACGGTTGcc                                                                             <  6:377792/147‑1 (MQ=255)                                                                       aCGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGc                                                                          <  6:214935/149‑1 (MQ=255)                                                                         gcgcTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCGTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCAt                                                                        <  3:109128/149‑1 (MQ=255)                                                                            cTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCtt                                                                     <  4:156472/149‑1 (MQ=255)                                                                            cTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCtt                                                                     <  4:133124/149‑1 (MQ=255)                                                                               gacTTAG‑‑‑‑ACCTGAATGCGGTGGCTCTCTGCGCGGATGGCCCGATTGTGCACGAAGATGTGATCGCCGGAGTCTCGGGAGTTGTTACGAAATCCTTTCGTCGTGTTTAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTtgt                                                                   <  4:238667/146‑1 (MQ=255)                                                                               gCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGTGGCCGGTATCACGGTTGCCGGCATCTtgt                                                                   <  1:215315/149‑1 (MQ=255)                                                                                ccTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGc                                                                 <  8:267449/149‑1 (MQ=255)                                                                                          aCCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACgtg                                                           >  5:126963/1‑149 (MQ=255)                                                                                            cAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGtt                                                         >  2:182716/1‑149 (MQ=255)                                                                                               aaTGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAg                                                      >  3:65409/1‑149 (MQ=255)                                                                                                   cGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCAGTTGACGGT‑GCCGGTATCACGGTTGCCCGCATCATGTGCCACCTGTTGAGTGcc                                                  >  6:382614/1‑149 (MQ=255)                                                                                                           cTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCt                                          >  6:82073/1‑149 (MQ=255)                                                                                                           cTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACC‑AAGCCTTTCGACGTGtt                                                                                                               <  7:363059/80‑1 (MQ=255)                                                                                                           cTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACC‑AAGCCTTTCGACGTGtt                                                                                                               >  8:363059/1‑80 (MQ=255)                                                                                                                  gtgagggCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGt                                   <  5:132013/144‑1 (MQ=255)                                                                                                                       ggCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGt                                                     >  4:8256/1‑126 (MQ=255)                                                                                                                       ggCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGt                                                     <  3:8256/126‑1 (MQ=255)                                                                                                                        gCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGtgctgc                              <  3:96970/148‑1 (MQ=255)                                                                                                                          ccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTgc                                       <  8:108082/137‑1 (MQ=255)                                                                                                                          ccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCggg                           >  2:103302/1‑149 (MQ=255)                                                                                                                          ccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCggg                           >  7:15729/1‑149 (MQ=255)                                                                                                                          ccGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCggg                           <  7:69707/149‑1 (MQ=255)                                                                                                                                gTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGc                     >  2:49208/1‑149 (MQ=255)                                                                                                                                     cGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACg                <  4:65409/149‑1 (MQ=255)                                                                                                                                      tttggtttCATCTACCGAGACTCGGGCGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACgg               >  6:102789/8‑149 (MQ=255)                                                                                                                                       aaGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGc              <  6:126963/149‑1 (MQ=255)                                                                                                                                               aTCTCCCGAGACTCGGGCGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCAGCATCTTGTGc                                                                 <  5:99732/90‑1 (MQ=255)                                                                                                                                               aTCTCCCGAGACTCGGGCGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCAGCATCTTGTGc                                                                 >  6:99732/1‑90 (MQ=255)                                                                                                                                               aTCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGca      >  4:369621/1‑149 (MQ=255)                                                                                                                                                tCGCCCGAGTCTCGGGTGATGAAACCAAAGCCTTTCGTCGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCAGGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGcac     <  5:382614/149‑1 (MQ=255)                                                                                                                                                tCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCGTTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGcac     >  6:19069/1‑149 (MQ=255)                                                                                                                                                 cGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGccc                   <  4:325649/134‑1 (MQ=255)                                                                                                                                                 cGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGccc                   >  3:325649/1‑134 (MQ=255)                                                                                                                                                  gCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGcacg    <  1:103302/148‑1 (MQ=255)                                                                                                                                                   cccGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGcacggc  >  4:207378/1‑149 (MQ=255)                                                                                                                                                    ccGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGcacggc  >  3:188052/1‑148 (MQ=255)                                                                                                                                                    ccGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCCTGCTGCGGGCCAAGCCCACGCCGAGATGcacgga  >  6:369679/1‑147 (MQ=255)                                                                                                                                                    ccGAGACTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTgc                                       <  5:298835/111‑1 (MQ=255)                                                                                                                                                    ccGAGACTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTgc                                       >  6:298835/1‑111 (MQ=255)                                                                                                                                                         |                                                                                                                                                TCAACGGCGGTTTACTGCAACCACGTCTTCGGCCTGCAGGCCCTTGTCGCGGTGCATCACCGAAAACTCCACGCGCTGGCCTTCG‑‑‑‑ACCAGAATGCGGTGGCCCTCACCGCGGATGGCCCGAAAGTGCACGAAGATGTCATCGCCCGAGTCTCGGGAGATGAAACCAAAGCCTTTCGACGTGTTGAACCATTTGACGGT‑GCCGGTATCACGGTTGCCGGCATCTTGTGCCACGTGTTGAGTGCCGCGTGCCTTGCGCGTGCTGCGGGCAAAGCCCACGGCCAGATGCACGGC  >  NC_002947/1387896‑1388186

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