Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	472,849	A→G	G84G (GGA→GGG)	glnK →	nitrogen regulatory protein PII‑2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	472,849	0	A	G	98.0%	159.9 / NA	49	G84G (GGA→GGG)	glnK	nitrogen regulatory protein PII‑2
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (21/27);  minor base T (0/1);  total (21/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 = 9.99e-01

GTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGG  >  NC_000913/472607‑473089                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 gtgtccgtgttaatagaaCCATTAAAGCGAGAAGATGTTTGTGAAGAGTTATCTTCCATTTGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCAATCCTAGTTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCAATTAAGTGAACGATAAAGTCAGTAAGGCGGATTACACCGGTAAAATTga                                                                                                                                                                                                                                          <  2:96604/235‑2 (MQ=255)      cGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACg                                                                                                                                                                                                                                     >  1:844/1‑251 (MQ=255)            taTCAAACAATTCAAGCTGGAAGACGTTAGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCACAAGAAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGTTAACGTCCGTAAGGCGGCTTCAACCTGGAAAATTGACGCCGGCaca                                                                                                                                                                                                                                >  2:37856/2‑248 (MQ=255)                   ccATTCAAGCCGGAAGACGTTCGTGAAGCGTTATCTTCCAGTGGTATGCAGGGCCGGGCCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCCTGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACGCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTtcgt                                                                                                                                                                                                                        <  2:636162/251‑1 (MQ=255)                           ccTGGAAGACGTTCGTGCAGCGTGTTCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGCAAGGTTTCGGGCGTCAGAAAGAGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCTATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAAt                                                                                                                                                                                                                <  2:170318/250‑1 (MQ=255)                            cTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAAtt                                                                                                                                                                                                               >  1:589726/1‑251 (MQ=255)                              ggAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTg                                                                                                                                                                                                              <  1:144427/250‑1 (MQ=255)                              ggAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTg                                                                                                                                                                                                              <  1:579196/250‑1 (MQ=255)                                        cGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCat                                                                                                                                                                                                   <  1:259385/251‑1 (MQ=255)                                         gTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCatt                                                                                                                                                                                                  >  1:491326/1‑251 (MQ=255)                                          tGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGAAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCattc                                                                                                                                                                                                 >  1:319082/1‑251 (MQ=255)                                                gTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTGCGTCGCTGAATTGCAACGCGTCattcgtattc                                                                                                                                                                                           <  2:648662/251‑1 (MQ=255)                                                     cctcctatGGTCTTCGGGGCCTTACCGTCGCCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCACAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAAGTGCAACGCGTCATTCGTATTCGTAc                                                                                                                                                                                       <  1:1407/243‑1 (MQ=255)                                                        ccATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGc                                                                                                                                                                                   <  1:437719/251‑1 (MQ=255)                                                                 ttCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCcgac                                                                                                                                                                          >  1:602358/1‑251 (MQ=255)                                                                     gggCCTGACCGTCACCGAAGGGAAAGGTTTCGGGCGGCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAg                                                                                                                                                                      <  2:406559/251‑1 (MQ=255)                                                                        ccTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATTACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTCCACCGGGAAAATTGGCGGCGGCACAATCTTCGTCGCTGAATTGCAACGCGACATTCGTATTCGTACCGGCGAAGCCGACGAAacgg                                                                                                                                                                   >  2:489847/1‑251 (MQ=255)                                                                        ccTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAgcgg                                                                                                                                                                   <  1:235357/251‑1 (MQ=255)                                                                        ccTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAgcgg                                                                                                                                                                   >  1:145572/1‑251 (MQ=255)                                                                             ccGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTg                                                                                                                                                              <  1:190001/251‑1 (MQ=255)                                                                                 cACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGAAGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTAAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAACCTTCGTCGCTGACTTGCAACGCGTCATTTGGATTCGTCACGGCGAACCCGACGAAGCGgcgcgtaaat                                                                                                                                                          >  2:617469/1‑245 (MQ=255)                                                                                     gAAGTGAAAGGTGTCGGGCGTCAGAAAGGGCATGCCGAGCTGTATCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACGCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTTTAATctct                                                                                                                                                      <  2:138412/251‑1 (MQ=255)                                                                                           aaaGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGTCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGTATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGcaca                                                                                                                                                <  1:615949/251‑1 (MQ=255)                                                                                                    gggCGTCAGAAAGGGCATGCCGAGCTGTGCCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAg                                                                                                                                       <  1:91439/251‑1 (MQ=255)                                                                                                     ggCGTCAGAAAGGGCATGCCGAGCGGGACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATGGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAgg                                                                                                                                      <  2:602358/251‑1 (MQ=255)                                                                                                                gggCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATg                                                                                                                           <  1:564634/251‑1 (MQ=255)                                                                                                                         gAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGGTGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGa                                                                                                                                 >  1:421402/1‑236 (MQ=255)                                                                                                                         gAGCTGTACCGGGGGGCGGAAGACAGCGTCAATGTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGGTGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCGGAATTGCAACGCGTCATTCGTATGCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGa                                                                                                                                 <  2:421402/236‑1 (MQ=255)                                                                                                                                                          ttCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGCTGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCCACAGGAACGAAAAATGGAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCgctg                                                                                 >  2:68997/1‑248 (MQ=255)                                                                                                                                                                     aGTAAAAATTGATGCGTCGATTGCTGATGACCAAATCGATGAAGTGATCGATATCGTCAGTAAGGCGGCATACACCGGGAAAATTGGCGCCGGCAAAATCTTCGTCGCTGCATTGCAACGCGTCATTCGTATTCGAACCGGCAAAGCCGCCGAAGCGGCGCTGTAATCTCTGGCACACAGCAACCGGAACGAAAAATg                                                                                                                           >  2:640236/1‑198 (MQ=255)                                                                                                                                                                     aGTAAAAATTGAGGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGTCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCAGTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATg                                                                                                                           <  1:640236/198‑1 (MQ=255)                                                                                                                                                                         aaaaTTGATGTGGCGATTGCTGATGACCAACTCGATGCAGTGATCGATATCGTCAGGAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACAGGCGATGCtagatagggagggg                                                                  >  2:215904/1‑238 (MQ=255)                                                                                                                                                                         aaaaTTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCt                                                                               <  1:215904/238‑1 (MQ=255)                                                                                                                                                                              tGATGTGGCGATTGCTTATGACCAACTCGATGAAGTGATCGATATCGGCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTTAAGTGCAACGCTTCATTCGTATTCGTACCGGCGAAGCCGACGAAgc                                                                                                                                                                     <  2:60304/147‑1 (MQ=255)                                                                                                                                                                              tGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAgc                                                                                                                                                                     >  1:60304/1‑147 (MQ=255)                                                                                                                                                                              tGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGc                                                             >  1:164768/1‑251 (MQ=255)                                                                                                                                                                              tGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTAACGGCGAAGCCGACGCAGCGGCGCTGTAATCTCTGGCCCACAGCAACAGGTACGAAAAATGAAGATAGCGACGATAAAAAATGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGc                                                             >  2:482611/1‑251 (MQ=255)                                                                                                                                                                                 tgtgGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGc                                                          >  1:308630/1‑251 (MQ=255)                                                                                                                                                                                    ggctaTTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAATGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCagc                                                       <  2:145572/247‑3 (MQ=255)                                                                                                                                                                                      cGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTg                                                     >  1:562545/1‑251 (MQ=255)                                                                                                                                                                                               atgaCCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCg                                            >  1:354166/1‑251 (MQ=255)                                                                                                                                                                                                              aaGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATg                             <  1:214906/251‑1 (MQ=255)                                                                                                                                                                                                              aaGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGGCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATg                             <  1:44662/251‑1 (MQ=255)                                                                                                                                                                                                                            gTCAGTAAGGCGGCTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGAtt               <  1:37856/251‑1 (MQ=255)                                                                                                                                                                                                                                     gcggcTTACACCGGGAAAATTGGTGACGGCAAAAGCTTCGTCGCGGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTgcg      <  2:491326/251‑1 (MQ=255)                                                                                                                                                                                                                                       ggcTTACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCt    <  1:220646/251‑1 (MQ=255)                                                                                                                                                                                                                                         cTTACACCGGGACAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATACGAACAGGCGACACCGACGAAGCGGCGCTGTAATATCTGGCACACAGCACAAGGAACGAAACATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGCTGCTTCCGGGCCTGCAAATGACTGCACTTGCGGTGTCGGATAAAGCCGACCATGCGTTAATGGTGCGTTGTACttccctgg  >  2:105666/1‑244 (MQ=255)                                                                                                                                                                                                                                           tACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAGCGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACTAAAAATGAAGATAGCGACGATAAAAACTGGgcttgctt                                                                                            <  2:434078/159‑1 (MQ=255)                                                                                                                                                                                                                                           tACACCGGGAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGgcttgctt                                                                                            >  1:434078/1‑159 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 GTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGATCGATATCGTCAGTAAGGCGGCTTACACCGGAAAAATTGGCGACGGCAAAATCTTCGTCGCTGAATTGCAACGCGTCATTCGTATTCGTACCGGCGAAGCCGACGAAGCGGCGCTGTAATCTCTGGCACACAGCAACAGGAACGAAAAATGAAGATAGCGACGATAAAAACTGGGCTTGCTTCACTGGCGATGCTTCCGGGACTGGTAATGGCTGCACCTGCGGTGGCCGATAAAGCCGACAATGCGTTTATGATGATTTGTACTGCGCTGG  >  NC_000913/472607‑473089

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