Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,554,539	T→G	Y485* (TAT→TAG)	malT →	DNA‑binding transcriptional activator MalT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,554,539	0	T	G	98.1%	167.4 / NA	52	Y485* (TAT→TAG)	malT	DNA‑binding transcriptional activator MalT
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (24/27);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTATAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGAAAAAGCATTCCAGCTGATCAACGAGCAGCATCTGGAACAGCTGCC  >  NC_000913/3554292‑3554784                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      acatcTTGGTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGGGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCGAGCCCGTGCTGAACATGAAATCAAGGACAGCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGATCGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCTGGCTGGTTCTAGAGc                                                                                                                                                                                                                                                    <  2:397872/247‑1 (MQ=255) aCAGCTTGCTGGAAAAGCCGCAGTGGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGGTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGAGTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCTCCGGGCTGTTTCTAGAGc                                                                                                                                                                                                                                                    <  2:397889/251‑1 (MQ=255)           ggAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGCTGCCAAGCCAACATCGCTACAGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGCAGACACCCTGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCCTTGTgg                                                                                                                                                                                                                                          >  1:400226/1‑251 (MQ=255)                           gtgtTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTggg                                                                                                                                                                                                                          <  2:365118/251‑1 (MQ=255)                                  tGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGGAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCCACCTCGGTGCTGGGGGAAggg                                                                                                                                                                                                                   >  2:83396/1‑249 (MQ=255)                                      ggCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGtgctgc                                                                                                                                                                                                               <  1:357592/251‑1 (MQ=255)                                            gCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACAGCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGGCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGTGCTGGTTCTATAGCCGCATTGTGGCAACCTCGGTGCTGGGTGGAGTGctgcactgca                                                                                                                                                                                                         <  1:129841/251‑1 (MQ=255)                                                 tGCAAAGCCAACATCGCTACGGCGAAGTTAACACACTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGc                                                                                                                                                                                                    >  1:61824/1‑251 (MQ=255)                                                  gCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAAGCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCTCTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGc                                                                                                                                                                                                    <  2:304027/250‑1 (MQ=255)                                                       gCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTg                                                                                                                                                                                              >  1:171807/1‑251 (MQ=255)                                                        ccAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGa                                                                                                                                                                                             <  1:230228/251‑1 (MQ=255)                                                            cATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGACATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCg                                                                                                                                                                                         >  2:101630/1‑251 (MQ=255)                                                                    cGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTgg                                                                                                                                                                                 >  1:13900/1‑251 (MQ=255)                                                                     ggCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGgc                                                                                                                                                                                <  2:100690/251‑1 (MQ=255)                                                                          aGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTa                                                                                                                                                                            >  1:367217/1‑250 (MQ=255)                                                                             tAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATgc                                                                                                                                                                        >  1:223773/1‑251 (MQ=255)                                                                                      gCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGCTGAAGCGGACCGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGGGGCAACCTCGGTGCTGGGAGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCCGCAAACCg                                                                                                                                                               >  2:406059/1‑251 (MQ=255)                                                                                      gCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAAGTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGTATCGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCTAAGTCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAAc                                                                                                                                                                 <  1:381886/249‑1 (MQ=255)                                                                                                gCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGc                                                                                                                                                     <  1:175278/251‑1 (MQ=255)                                                                                                      cATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGcc                                                                                                                                                >  1:202743/1‑250 (MQ=255)                                                                                                         gAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGATGGCAAAAATGGCACTGGAAGAGCCGCCGCCGGGATGGATCTAGAGCCGCATTGTGGCAACATCGGTGCTGGGGGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCAAGCCAGc                                                                                                                                             >  1:35265/1‑250 (MQ=255)                                                                                                                     aTCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGc                                                                                                                                 >  1:372757/1‑250 (MQ=255)                                                                                                                         gagaAGACACGATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGACGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGAGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTAc                                                                                                                            >  2:124717/1‑251 (MQ=255)                                                                                                                                   cATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGATCAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGa                                                                                                                  <  1:240941/251‑1 (MQ=255)                                                                                                                                           cAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACGGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATc                                                                                                          <  1:358056/251‑1 (MQ=255)                                                                                                                                             gAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCa                                                                                                        <  1:101630/251‑1 (MQ=255)                                                                                                                                                   aaCGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGGGCTGCACTGCAAAGGCGAATTTACCCGCTCACTGGCTCTAATGCAGCAAACCGAACATATGGCACGCCAGCGCGATGTCTGGCACTACGCTTTGTGGTGTTTAATCCAGCAAAg                                                                                                  <  2:419301/251‑1 (MQ=255)                                                                                                                                                       ctctGTGCGCCCAGGTGGCGATTAACGCTGGTAATCCGGATGAAGCGGAACGGCTGGCAACACTGGCACTGGAAGAGATGCCGCCGGGCTGGTCCTCGAGCCGCCTTGTGGCGACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAAATGGCCCGCCAGCACGCTGTCTGGCCCTAGGCTTTGTGGAGTTTAATCACGCAAAGTGa                                                                                               >  2:365872/1‑250 (MQ=255)                                                                                                                                                             gcgcCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTg                                                                                        >  1:52794/1‑251 (MQ=255)                                                                                                                                                              cgcCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGAGGCAACTTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGt                                                                                       >  1:226031/1‑251 (MQ=255)                                                                                                                                                               gcCCAGTTGGCGATTAACGATGGTAATCCGGGTGAAGGGGACCTGCTGGCAAAACGGGGACTGGAAGAGTTGCCGCCGGGCGGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCGCACTGGCGCGAAGGCAGCAAAGCGAACAGATGGCAGGCCAGCACGATTTCTGGCACTACGCTTTTTGGAGTTTAATGCAGCAAAGTGAAATTCTGtt                                                                                      <  2:426001/251‑1 (MQ=255)                                                                                                                                                                cccAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGttt                                                                                     <  1:17242/251‑1 (MQ=255)                                                                                                                                                                    ggTGGCGATTAACGATGGTAATTCGGATGAAGGGGGACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGGGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGTTCACTGGCGCTAATGCAGCAAGCCGAACAGATGGCTCGCCAGCACGATGTCTGGCACTACGCTGTGATGAGTTTAATCCAGCAAAGTGAAAGTCTGTTTGccc                                                                                 <  2:439769/251‑1 (MQ=255)                                                                                                                                                                    ggTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCGAATGCAGTAAACCGAACAGATGGCACGCCAGCAGGATGTCTGGCACTACGCGTTGGGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGccc                                                                                 <  2:36690/251‑1 (MQ=255)                                                                                                                                                                          gATTAACGATGGTAATCCGGATGAAGGGGAACGGCTGGCGAAAGTGGCACTGGAAGGGCTGCCGCCGGGCTGGTTGTAGAGGCGCATTGTGGCAACCGCGGTGCTGGGTGAAGTGCTGCACTGCAAGGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTGTAATCCAGCAAAGTGAAATTCTGTTGGCCCAAGGGt                                                                           <  2:171318/251‑1 (MQ=255)                                                                                                                                                                              agtgatgggaatcaggagtaatcGGAACGGCTGTCAAAACGGGCACTGGAAGAGCGGCCGCCGGGCTGGTTCTAGATCCGCACTGGGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTTCAAAGGCGAATTGACCCGCTCACTGTCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTGTGTGGAGTTTAATCCGGCAAAGTGAAATTCAGTTTTTCCCAGGGTTcc                                                                        <  2:150521/228‑1 (MQ=255)                                                                                                                                                                                  aTGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGGGGCAACCTCGGGGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTATGGCACTACGCCTTGTGGAGCTTAATCCAGCAAGGTGAAATTCTGTTTGCCCAAGGGTTCCTGCac                                                                   >  2:7161/1‑250 (MQ=255)                                                                                                                                                                                          ccGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGt                                                            <  1:36626/250‑1 (MQ=255)                                                                                                                                                                                             gATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTTGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGACACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGa                                                        <  1:99376/251‑1 (MQ=255)                                                                                                                                                                                                      gAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGGTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGa                                               <  1:119744/251‑1 (MQ=255)                                                                                                                                                                                                        aCGGCTGGCAAAACTGGCACTGGAAGAGCTCCCGCCGGGATGGTTCTAGAGCCGCCTTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCAATGCAAAGGCGAATTGACACGCTCCCTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCCTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGaa                                              >  1:264018/1‑250 (MQ=255)                                                                                                                                                                                                           gCTGGCAAAACTGGCACTGGAAGAGCTGCAGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGAGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCCAAGTGACATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGaaaa                                            >  2:105293/1‑249 (MQ=255)                                                                                                                                                                                                                         cACTGGAAGAGCTGCCGCCGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGAAAAAGCATTCCAGCTGATc                            >  1:369395/1‑251 (MQ=255)                                                                                                                                                                                                                            tGGAAGCGCCGCCGCCGGGCTGGTCCTAGAGCCGCATTGTGGCAACCTCGGTGCCGCCTGAAGTGCTGCCCTGCAAAGCCGCATTGACCCGCTCACTGGCGCTAACGCAGCAACCCGAACAGATGCCACGCCCGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCCACGTGAAATTCTGTTTGCCCAAGGGTTCGTGCAAACAGCGTGGGAAACGCAGGAAAAAGAATTCAAGCTGATAAAc                         >  2:355376/1‑251 (MQ=255)                                                                                                                                                                                                                              gAAGAGCTAACGCCGGGCGGGATCTAGAACCGCATTGCGGCAACATCGCAGCTGGGCGCAGAGCCGCACTGCAAAGGCGAATTGAACCGCTCACTGGCGCTAATCCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGCGGAGTTTACTCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGATCCTGCAAACCGCGTGGGCAACGCAGGAAAAAGCCTTCCAGCTGATCAACgc                       >  2:409819/1‑250 (MQ=255)                                                                                                                                                                                                                                      gccgccGGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCGTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACTCAGGAAAAAGCATTCCAGCTGATCAACGAGCAGCAt                <  1:365872/250‑1 (MQ=255)                                                                                                                                                                                                                                         gccgGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTTGGTGAAGTGCTGCGATGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCCTGCCAGCTCGATGTCGGGCACTACGTTTTGGGGAGGTTAATCCAGCAAAGTGAAATGCTGGTTGCCCAAGGGTTCGTGCAAACCGCGTGGGAAACTCAGGAAAAAGCATTCCAGCTGATCAAAGAGCAGCATCTgg            <  2:252622/251‑1 (MQ=255)                                                                                                                                                                                                                                         gccgGGCTGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGAAAAAGCATTCCAGCTGATCAACGAGCAGCATCTgg            <  1:83396/251‑1 (MQ=255)                                                                                                                                                                                                                                           cgGGCAGGTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAACAGCAGGAAAAAGCATTCCAGCTGATCAACGAGAAGCATCTGGa           >  2:334514/1‑250 (MQ=255)                                                                                                                                                                                                                                                 ggTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGAAAAAGCATTCCAGCTGATCAACGAGCAGCATCTGGAACAGCTg    >  1:337514/1‑251 (MQ=255)                                                                                                                                                                                                                                                  gTTCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGAAAAAGCAGTCCATCTGATCAACGAGCAGCATCTGGATCAGCTGc   <  2:41638/251‑1 (MQ=255)                                                                                                                                                                                                                                                    tCTAGAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGAAAAAGCATTCCAGCTGATCAACGAGCAGCATCTGGAACGGCTGcc  <  1:124717/250‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                      ACAGCTTGCTGGAAAATCCGCAGTTGGTGTTATTGCAGGCGTGGCTGATGCAAAGCCAACATCGCTACGGCGAAGTTAACACCCTGCTAGCCCGTGCTGAACATGAAATCAAGGACATCAGAGAAGACACCATGCACGCAGAATTTAACGCTCTGCGCGCCCAGGTGGCGATTAACGATGGTAATCCGGATGAAGCGGAACGGCTGGCAAAACTGGCACTGGAAGAGCTGCCGCCGGGCTGGTTCTATAGCCGCATTGTGGCAACCTCGGTGCTGGGTGAAGTGCTGCACTGCAAAGGCGAATTGACCCGCTCACTGGCGCTAATGCAGCAAACCGAACAGATGGCACGCCAGCACGATGTCTGGCACTACGCTTTGTGGAGTTTAATCCAGCAAAGTGAAATTCTGTTTGCCCAAGGGTTCCTGCAAACCGCGTGGGAAACGCAGGAAAAAGCATTCCAGCTGATCAACGAGCAGCATCTGGAACAGCTGCC  >  NC_000913/3554292‑3554784

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