Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,627,299	2 bp→GA	intergenic (+5/‑15)	radA → / → nadR	DNA recombination protein/DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,627,299	0	A	G	98.5%	220.4 / ‑4.8	68	intergenic (+5/‑16)	radA/nadR	DNA recombination protein/DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (30/37);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01
*	NC_000913	4,627,300	0	T	A	98.5%	215.6 / NA	68	intergenic (+6/‑15)	radA/nadR	DNA recombination protein/DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base A (30/37);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑ATATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTAGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCA  >  NC_000913/4627056‑4627539                                                                                                                                                                                                                                                     ||                                                                                                                                                                                                                                                tttCCAGACTGCGTGACAGACCGCTGCCACCGGACCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGCACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTAACGAAAAACGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACgg                                                                                                                                                                                                                                           >  2:60270/1‑251 (MQ=255)   tCCAGCCTGCGGGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACgga                                                                                                                                                                                                                                          <  1:356892/250‑1 (MQ=255)    ccAGCCTGCGTGACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACggagg                                                                                                                                                                                                                                        >  2:173334/1‑251 (MQ=255)           gcgtgACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGTGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGat                                                                                                                                                                                                                                   <  1:364362/245‑1 (MQ=255)           gcgtgACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGatat                                                                                                                                                                                                                                 >  1:298425/5‑251 (MQ=255)             gtgACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATgt                                                                                                                                                                                                                               <  1:492208/249‑1 (MQ=255)             gtgACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATgt                                                                                                                                                                                                                               >  1:197608/3‑251 (MQ=255)               gACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGAAGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATgtcg                                                                                                                                                                                                                             >  2:345190/1‑251 (MQ=255)                      cGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTGTTGGCGTTAAAAACCTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTg                                                                                                                                                                                                                      <  2:19933/251‑1 (MQ=255)                            cacAGGATCTGGTGGTGTTTGGTGAAGTCGGGATGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGCTTTTTGGCGTTAAAAAACTCTCCGCCGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTAcc                                                                                                                                                                                                                >  1:156094/3‑251 (MQ=255)                             acAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCt                                                                                                                                                                                                               <  2:298425/250‑1 (MQ=255)                              cAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAAATCAGCCAGGGGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGATTAGCGGGTTCGCAGACTTATAATGGG‑GAATACGGAGAGAGAAATGTAGTCATTTGCTTACCTg                                                                                                                                                                                                              >  2:356633/1‑251 (MQ=255)                                 gATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGaa                                                                                                                                                                                                            >  1:342694/1‑250 (MQ=255)                                 gATCTGGAAGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCAAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAACACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGGGGGAGCTATGTAGTCATTTGATTACCTGaa                                                                                                                                                                                                            >  2:421016/1‑250 (MQ=255)                                           ggtTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCa                                                                                                                                                                                                 <  2:156094/250‑1 (MQ=255)                                            gtTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCa                                                                                                                                                                                                 <  1:189968/250‑1 (MQ=255)                                              ttGGTGAAGTCGGGCTGGCAGGGGAGATCCTGCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCGGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGGTAAAAAACGGTCGGACGCGCGTAGCGTGTTCGGCGGCTTATAAGGAG‑GAATGCGGAGGGGGATATGTCGTCAGTTGATTGCCTGAAATCTGcctgcaag                                                                                                                                                                                               <  2:260565/250‑7 (MQ=255)                                               tGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGc                                                                                                                                                                                              <  1:4575/250‑1 (MQ=255)                                                  tGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCACCGGGCGATTGTTCCGGCGGCTAACGCACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGACGCCATTTGATTACCTGAAAACTGCCATCAAGCAAc                                                                                                                                                                                           >  1:597657/1‑250 (MQ=255)                                                  tGAAGTCGCACTGGCAGGGGGGATCCGCCCGGTGCCCAGCGGACCGGAACGAATTTCAGCAGCGGCGACACCAGGTTTTCGCCGGGCGCTTGTTCCGGCGGCTCACGCACCGCAACAAACGCCGGAAGGGACGCCGATTTTTTGCGTTAAAAGACTCTCGGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAT‑TACGGCGGGGGATCTGTCGTCCTTTGATTCACTGATAACTTCCCTCAAGCAAc                                                                                                                                                                                           >  2:234432/1‑250 (MQ=255)                                                    aaGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAgg                                                                                                                                                                                        <  1:162447/251‑1 (MQ=255)                                                      gTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGc                                                                                                                                                                                      >  1:603184/1‑251 (MQ=255)                                                           gCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCAc                                                                                                                                                                                 <  1:280320/251‑1 (MQ=255)                                                                 aGGGGAGAGCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACCCGGTTTTCGCCGGGCGAGTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGCAAACTGCCATCAAGCAACATGGCTGCACGCTAca                                                                                                                                                                           <  2:277706/251‑1 (MQ=255)                                                                   gggAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCAGTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTAcagc                                                                                                                                                                         <  1:37388/251‑1 (MQ=255)                                                                       gaTCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGt                                                                                                                                                                     <  1:616319/251‑1 (MQ=255)                                                                         tCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGCGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGtgg                                                                                                                                                                   >  2:69139/1‑249 (MQ=255)                                                                             cccGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTg                                                                                                                                                                >  1:238350/1‑250 (MQ=255)                                                                                ggTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGc                                                                                                                                                            <  2:12899/251‑1 (MQ=255)                                                                                 gTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGcc                                                                                                                                                           >  1:268103/1‑251 (MQ=255)                                                                                     ccAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCTTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCGCGCTGCAGCAGGTGGCTGATGCTAGCg                                                                                                                                                       <  2:396074/251‑1 (MQ=255)                                                                                       aGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGt                                                                                                                                                     >  2:101371/1‑251 (MQ=255)                                                                                          ggTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATg                                                                                                                                                  <  1:110846/251‑1 (MQ=255)                                                                                              aGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGAcc                                                                                                                                               >  1:623805/1‑250 (MQ=255)                                                                                                   cGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAggg                                                                                                                                         <  1:270827/251‑1 (MQ=255)                                                                                                         tcagAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTAttt                                                                                                                                    >  1:380432/4‑250 (MQ=255)                                                                                                         tcagAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTAttt                                                                                                                                    <  2:487810/247‑1 (MQ=255)                                                                                                                    gAAACACGGTTTGCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGGCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCCGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACt                                                                                                                        <  2:55617/251‑1 (MQ=255)                                                                                                                        cacGGTTTTCGCCGGGCGATTGGTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAAt                                                                                                                    <  1:253306/251‑1 (MQ=255)                                                                                                                            gTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCa                                                                                                                <  1:173334/251‑1 (MQ=255)                                                                                                                            gTTTGCGCCGGGCGATTGTTCCGGCGGCGAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCAACAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCa                                                                                                                <  2:251392/251‑1 (MQ=255)                                                                                                                                 cGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATc                                                                                                           >  1:198959/1‑251 (MQ=255)                                                                                                                                              ttCCGGCGGCTAACGTACCGCAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAgcg                                                                                              >  1:461850/1‑251 (MQ=255)                                                                                                                                                  ggcggcTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTAGTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCa                                                                                          <  2:51972/251‑1 (MQ=255)                                                                                                                                                    cggcTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAACATCAAAAGCCCCCgcgcgcca                                                                                         >  1:617525/1‑248 (MQ=255)                                                                                                                                                           cgtaccGAAAAAAGCGACGGAAGGGATGCAGAGTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGa                                                                                 <  2:53969/247‑1 (MQ=255)                                                                                                                                                            gtaccGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCa                                                                                                                                                                                  >  2:319468/4‑153 (MQ=255)                                                                                                                                                            gtaccGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCa                                                                                                                                                                                  <  1:319468/150‑1 (MQ=255)                                                                                                                                                              accGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGc                                                                              >  1:382390/2‑251 (MQ=255)                                                                                                                                                                          gcCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGa                                                                                                                                                 <  1:388405/172‑1 (MQ=255)                                                                                                                                                                          gcCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGa                                                                                                                                                 >  2:388405/1‑172 (MQ=255)                                                                                                                                                                          gcCGGAAGGGATGCAGATTTTTGGCGGTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACGGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCGCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGttt                                                                  <  1:234432/251‑1 (MQ=255)                                                                                                                                                                                  ggATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGAGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGc                                                          >  1:533780/1‑251 (MQ=255)                                                                                                                                                                                   gATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCt                                                         <  1:110918/251‑1 (MQ=255)                                                                                                                                                                                                   gTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGc                                           >  2:620646/1‑249 (MQ=255)                                                                                                                                                                                                   gTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCag                                         <  2:603184/251‑1 (MQ=255)                                                                                                                                                                                                     tAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGCGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTCCAGCAGGTGGCTGATGCCAGCGGAATGACCAAAGGGGATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCaga                                        >  2:458885/1‑250 (MQ=255)                                                                                                                                                                                                       aaaaaCTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAACGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCaga                                        >  2:603763/1‑248 (MQ=255)                                                                                                                                                                                                           aCTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAAc                                 <  1:48942/251‑1 (MQ=255)                                                                                                                                                                                                            ctctCCGACGCGCTTAGAGTGTTCGACGACTTATAATGAG‑GAATACGGCGGGCGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCACCCGGGCTGAACGCTACAGCAGGCGGCTGATGCCAGCGGAATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAACAGCCCCAGCGCGCAAAAGATGAAGACGTTGCACCGTTTTTTGGGGCTTGAGTTTCCACGGCAGAAGAAAAcc                                >  1:304990/1‑250 (MQ=255)                                                                                                                                                                                                                ccGACGCGCTTAGCGTGTGCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTGGATGACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGGGGCTGATGCCAGCGGTATGACCAAAGGGTATGTAAGCCAGTTACGGAATGCCAAAATTAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTGCCCGGCAGAAGAAAAccat                              <  2:274052/249‑4 (MQ=255)                                                                                                                                                                                                                    cgcgCTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATAGGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGgtgt                        <  1:614623/251‑1 (MQ=255)                                                                                                                                                                                                                      cgcTTAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTGGATTACCGGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTTCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGgtgtgg                      <  1:313269/251‑3 (MQ=255)                                                                                                                                                                                                                          tAGCGTGTTCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTTGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAATAAAACCATTGGTGTGGTAt                   <  2:377975/250‑1 (MQ=255)                                                                                                                                                                                                                                 ttCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTAGTCGGTaa            <  1:83395/250‑1 (MQ=255)                                                                                                                                                                                                                                 ttCGACGACTTATAATGAG‑GAATACGGAGGGAGATATGCCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGAACAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTaa            >  2:83395/1‑250 (MQ=255)                                                                                                                                                                                                                                   cgacgaCTTATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGtt         <  1:92442/251‑1 (MQ=255)                                                                                                                                                                                                                                          ttATAATGAG‑GAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTAcccg  >  1:336201/1‑250 (MQ=255)                                                                                                                                                                                                                                                     ||                                                                                                                                                                                                                                                TTTCCAGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAG‑ATATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTAGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCA  >  NC_000913/4627056‑4627539

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