Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	524,839	G→A	D527N (GAT→AAT)	rhsD →	protein RhsD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	524,839	0	G	A	98.1%	170.3 / NA	53	D527N (GAT→AAT)	rhsD	protein RhsD
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (24/28);  minor base C (1/0);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GGACAGAGTACGGTCTGAATGTGGTGTCCGGCGATATCACGGACATCACCACACCGGACGGGCGGGAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCATCAGCCTTTACCGCCGCTATGACAACCGTGGCCGGTTAACCTCGGTGAAAGACGCACAGGGCCGTGAAACGCGGTATGAATACAACGCCGCAGGCGACCT  >  NC_000913/524594‑525077                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               ggaCAGAGTACGGTCTGAATGTGGTGTCCGGCGATATCACGGACATCACCACACCGGACGGGCGGGAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCg                                                                                                                                                                                                                                           <  2:22892/251‑1 (MQ=255)          aCGGTCTGAATGTGGTGTCCGGCGATATCACGGACATCACCACACCGGACGGGCGGGAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGgcagc                                                                                                                                                                                                                                  >  1:358338/1‑251 (MQ=255)             gTCTGAATGTGGTGTCCGGCGATATCACGGACATCACCACACCGGACGGGCGGGAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCAcc                                                                                                                                                                                                                               >  1:319743/1‑251 (MQ=255)                        gtgtCCGGCGATATCACGGACATCACCACACCGGACGGGCGGGAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGAc                                                                                                                                                                                                                    <  2:324776/251‑1 (MQ=255)                          gtCCGGCGATATCACGGACATCACCCCACCGGACGGGCGGGAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCt                                                                                                                                                                                                                  <  2:316378/251‑1 (MQ=255)                                cGATATCACGGACATCACCACACCGGACGGGCGGGAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGcc                                                                                                                                                                                                            <  2:227959/251‑1 (MQ=255)                                 gATATCACGGACATCACCACACCGGACGGGCGGGAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCg                                                                                                                                                                                                           >  1:124213/1‑251 (MQ=255)                                       aCGGACATCACCACACCGGACGGGCGGGAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACgg                                                                                                                                                                                                     >  1:330886/1‑251 (MQ=255)                                       aCGGACATCACCACACCGGACGGGCGGGAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCCACCGCTACGCTGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCCGCAGACGAACTGGAGGCGCTACgg                                                                                                                                                                                                     >  2:21270/1‑251 (MQ=255)                                         ggACATCACCACACCGGACGGGCGGGAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGGGAGTTACCGGCGGCGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTCCGGGc                                                                                                                                                                                                   >  2:104735/1‑251 (MQ=255)                                          gccATCCCCACACCGGACGGGCGGGAGACGAAATTTTACTATAACGACGGGAACCAGCGGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCa                                                                                                                                                                                                  <  2:73864/249‑1 (MQ=255)                                            cATCACCACACCGGACGGGCGGGAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGt                                                                                                                                                                                                <  2:358338/251‑1 (MQ=255)                                             aTCACCACACCGGACGGTCGTGAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGTCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCGACGATGACGCGCACAGTGAGTTACCGGCGACGACAACTAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGtt                                                                                                                                                                                               <  1:359114/251‑1 (MQ=255)                                             aTCACCACACCGGACGGGCGGGAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGtt                                                                                                                                                                                               <  2:72643/251‑1 (MQ=255)                                                    cacCGGACGGGCGGGAGACGAAATTTTACTATAACGACGGGACCCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCg                                                                                                                                                                                        <  2:235873/251‑1 (MQ=255)                                                     cccGGACGGGCGGGAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGt                                                                                                                                                                                       <  2:324696/250‑1 (MQ=255)                                                       cGGACGGGCGGGAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCAACCGGCAGTTGCTGGAGttt                                                                                                                                                                                     >  2:206370/1‑250 (MQ=255)                                                        ggACGGGCGGGAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCa                                                                                                                                                                                    <  2:355765/251‑1 (MQ=255)                                                          aCGGGCGGGAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCAGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCAcc                                                                                                                                                                                  <  1:8923/251‑1 (MQ=255)                                                            gggcgggAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCg                                                                                                                                                                                 >  1:337823/1‑250 (MQ=255)                                                                  gagaCGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTc                                                                                                                                                                          >  1:77385/1‑251 (MQ=255)                                                                      cGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGc                                                                                                                                                                      <  1:323180/251‑1 (MQ=255)                                                                              tACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGAc                                                                                                                                                              >  1:205594/1‑251 (MQ=255)                                                                                          gggAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGACGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAAt                                                                                                                                                   >  1:47753/1‑250 (MQ=255)                                                                                                  gCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCt                                                                                                                                          >  1:330245/1‑251 (MQ=255)                                                                                                      aCGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACTAATGCGACGGGCAGCACCCTGCAGATGACCGTGAGCCTCTACGGGCAGTTGCTGGCGTTCACCGACTGCGCGGGCTACCAGACCCGTTAGGAAGACGACCGCTTCgg                                                                                                                                      <  2:192835/251‑1 (MQ=255)                                                                                                               gtgtCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGAc                                                                                                                             >  1:42531/1‑251 (MQ=255)                                                                                                                 gtCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGAcgg                                                                                                                           <  1:78586/251‑1 (MQ=255)                                                                                                                 gtCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGAcg                                                                                                                            <  1:180582/250‑1 (MQ=255)                                                                                                                   ccccGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGAGGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGAcggcg                                                                                                                         <  2:280983/251‑1 (MQ=255)                                                                                                                     ccGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCGGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGt                                                                                                                       <  2:49928/251‑1 (MQ=255)                                                                                                                                      ccgccgGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGc                                                                                                       >  2:96165/1‑250 (MQ=255)                                                                                                                                              aaTATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTCTGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGCGGAACGCCTCCGCCtt                                                                                              >  2:332236/1‑251 (MQ=255)                                                                                                                                                   gatgaACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCATCAGCCTTTAcc                                                                                          <  2:14781/250‑1 (MQ=255)                                                                                                                                                      gaACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCATCAGCCTTTAccgcc                                                                                       <  2:337823/250‑1 (MQ=255)                                                                                                                                                          cGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCATCAGCCTTTACCGCCGCTAt                                                                                  >  2:21723/1‑251 (MQ=255)                                                                                                                                                            ggCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCATCAGCCTTTACCGCCGCTATg                                                                                 >  1:86530/1‑250 (MQ=255)                                                                                                                                                                     gTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCATCAGCCTTTACCGCCGCTATGACAACCGTgg                                                                       >  2:284709/1‑251 (MQ=255)                                                                                                                                                                       atgggagccatcccgcgggggggAGACAGTACGCTCCCGCTACGAGGACGCGCCAAGTGAGTTCCCGCCACCGACAACGAACGCGACGGGAAGCCCCCGGCAGATGCCCGGGAGCGGCTACGGGCAGTTGCTGGCGTTCCCCGACTGCTCGGCCTACCAGGCCCTTTATGAATGCGCCCGCTTCGGCCAAATGACGCCGTTCCGCCGGGAGGAAGTCATCAGCCTTTACCGCCGCTATGACAACCGTGGc                                                                      <  1:362975/231‑1 (MQ=255)                                                                                                                                                                                atcgcgaggCGGGGAGACAGTACGCTCCCGCTACGATGACGCGCACAGTGATTTACCGGCGACGACAACGAAGGCGACGGGCAGCACCCGGCAGATGCCCTGGGGCCGCTACGGGCAGTGGCGGGCGTTCACCGTCTGCTCGGGCTCCCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCCCCGCGAGGAAGGCATCAGCCTTTACCGCCGCTATGACAACCGTGGCCGGTTAACCt                                                            <  1:187311/243‑1 (MQ=255)                                                                                                                                                                                  cgcgCAGCGGGGAGACAGTACGCTACCGCTACGCAGCAGCGCCCAGTGTGCTCCCGCAGACGACAACGAATGCGACGCGCAGCACCCGGACGATGACCCGGAGCCGCTCCGGGCAGCTGCTGGCGTTCACCGCCTGCTAGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGCAGGCATCAGCCTTTACCGCCGCTATGACAACCGTGGCCGGTTAACCTcc                                                          >  1:353023/1‑250 (MQ=255)                                                                                                                                                                                  cgcgCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCATCAGCCTTTACCGCCGCTATGACAACCGTGGCCGGTTAACCTCg                                                          >  2:70030/1‑251 (MQ=255)                                                                                                                                                                                        ggggagggaaatgccgccgccgctacAGTGACGCGCACGGTGAGTGCCCGGCGACGACAACGAAGGCGCGGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGCCAGTGGCTGGCGTTCCCCGACTGCGGGGGCCACAAGACCCGTAATAAATCCGACCGCTTGGGCCAGATGAGGGGGGTCCCCGGCGAGGAAGGCACCAGCTTTTACCGCCGCTATGACAACCGTGGCCGGTTAACCTCGGTGaaa                                                    <  2:203998/232‑1 (MQ=255)                                                                                                                                                                                             cagaCAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCGGGCGTTCACCGACTGCTCGGGCTACCAGCCCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGGCCACCGCGAGGAAGGCATCAGGCTTTACCGCCTCTATGACAACCGTGGCCTGTTAACCTCGGTGCAAGACGc                                               >  2:224477/2‑251 (MQ=255)                                                                                                                                                                                             cagaCAGTACGCTACCGCCACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCATCAGCCTTTCCCGCCGCTATGGCAACCGTGGCCGGTTAAACACGGTGAAAGaagc                                               >  2:224492/2‑248 (MQ=255)                                                                                                                                                                                                                       cgcgCACAGTGAGTTCCCGGCGACGACAACGAATGCGACGGGCGGCACCCGGCAGATGACCTGGAGCCGCTGCGGGCATTTGCGGGCGTTCCCCGACTGCTCGGGCTACCAGACCCTTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCCCCGCGAGGAAGGCATCAGCCTTTACCGCCGCTATGACAACCGTGGCCGGTTAACCTCGGTGAAAGACGCACAGGGCCGTGAAACGCGGTATGAAt                     <  1:327937/251‑1 (MQ=255)                                                                                                                                                                                                                        gcgcACAGTGGGTGCCCGGGGACGACAACGAATGCGCCGGGCAGCACCCGGCAGATGCCCGGGGGCCGCGCCGGGCAGTTGCTGGCTTCCCCCGACTGCCCGGGCTACCAGACCCGTTATGAATCCGCCCGCTTCGGCCAGATGCCGGCGGTCCACCGCGAGGAAGGAATCAGCTTTTACCGCCGCTATGACAACCGTGGCCGGTTAACCTCGGTGAAAGACGCACAGGGCCGTGAAACGCGGTATGAATa                    <  2:68242/251‑1 (MQ=255)                                                                                                                                                                                                                        gcgcACAGGGATTTACCGGCGCCGACAACGAAGGCGGCGGGCAGCCCCGGGCAGATGCCCGGGACCCGCTACGGGCATTGGCGGGCGTTCACCGACTGCGCGGGCTCCCAGCCCCTTTATGAATCCGCCCGCTTCGGCCAGATGCCGGCGTTCCACCGCGAGGAAGGCATCAGCCTTTACCGCCGCTATGACAACCGTGGCCGGTTAACCTCGGTGAAAGACGCACAGGGCCGTGAAACGCGGTATGAATa                    <  1:66551/251‑1 (MQ=255)                                                                                                                                                                                                                              aGTGAGTTACCGGCGACGACAACGCATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCATCAGCCTTTACCGCCGCTATGACAACCGTGGCCGGTTAACCTCGGGGAAAGACGCACAGGGCCGTGAAACGCGGGATGAATACAAcgc              >  2:219725/1‑251 (MQ=255)                                                                                                                                                                                                                                   gttcccGGCGACGACAACGAAGGCGAGGGCCAGCCCCCGGCAGATGACCGGGAGCCGCTACGGGCAGTTGCGGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCATCAGCCTTTACCGCCGCTATGACAACCGTGGCCGGTTAACCTCGGTGAAAGACGCACAGGGCCGTGAAACGCGGTATGAATACAACGCCGCAg         <  1:21270/247‑1 (MQ=255)                                                                                                                                                                                                                                   gTTACCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCATCAGCCTTTACCGCCGCTATGACAACCGTTGCCGGTTAACCTCGGTGAAAGACTCCCCGGCCCGTGAAACGCGGTAGGGATACAAcgccgc           >  2:99553/1‑249 (MQ=255)                                                                                                                                                                                                                                      aCCGGCGACGACAACGAATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCATCAGCCTTTACCGCCGCTATGACAACCGTGGCCGGTTAACCTCGGTGAAAGACGCACAGGGCCGTGAAACGCGGTATGAATACAACGCCGCCGGcc      >  2:197527/1‑250 (MQ=255)                                                                                                                                                                                                                                          gCGACGACAACGAATGCGACGGGCAGCCCCCGGCAGAGGCCCTGGAGCCGCTCCGGGCAGTTGCGGGCTTTCCCCGACTGCTCGGGCTACCAGACCCTTTATGAATACGCCCGCTTCGCCCAGATGACGGGGGTCCACCGCGAGGAAGGCATCAGCTTTTACCGCCGCTATGACAACCGTGGCCGGTTAACCTCGGTGAAAGACGCACAGGGCCGTGAAACGCGGTATGAATACAACGCCGCAGGCGACCt  <  1:102582/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                               GGACAGAGTACGGTCTGAATGTGGTGTCCGGCGATATCACGGACATCACCACACCGGACGGGCGGGAGACGAAATTTTACTATAACGACGGGAACCAGCTGACGGCGGTGGTGTCCCCGGACGGGCTGGAGAGCCGCCGGGAATATGATGAACCGGGCAGGCTGGTATCGGAGACATCGCGCAGCGGGGAGACAGTACGCTACCGCTACGATGACGCGCACAGTGAGTTACCGGCGACGACAACGGATGCGACGGGCAGCACCCGGCAGATGACCTGGAGCCGCTACGGGCAGTTGCTGGCGTTCACCGACTGCTCGGGCTACCAGACCCGTTATGAATACGACCGCTTCGGCCAGATGACGGCGGTCCACCGCGAGGAAGGCATCAGCCTTTACCGCCGCTATGACAACCGTGGCCGGTTAACCTCGGTGAAAGACGCACAGGGCCGTGAAACGCGGTATGAATACAACGCCGCAGGCGACCT  >  NC_000913/524594‑525077

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