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.1%	163.1 / ‑4.3	53	intergenic (+5/‑16)	radA/nadR	DNA recombination protein/DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (24/28);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01
*	NC_000913	4,627,300	0	T	A	98.1%	163.2 / NA	53	intergenic (+6/‑15)	radA/nadR	DNA recombination protein/DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base A (23/28);  minor base C (1/0);  total (24/29)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTAGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACC  >  NC_000913/4627061‑4627537                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                              ggCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGGGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTTTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACggaggga                                                                                                                                                                                                                                    <  2:205904/250‑1 (MQ=255)     tGCGTGACAGACCGCTGCCACAGGATCTGGTGGCGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCCGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGCTTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGAGGCAGATTTTTGGCGTTAAAAAACTCTCAGCCGCGCTTAGAGTGTTCGACGACTTATACTGAGGAATACGGAGGGGGat                                                                                                                                                                                                                                 >  1:28630/1‑250 (MQ=255)      gcgtgccagcccGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACTGTCTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATCCGGTGGGAGat                                                                                                                                                                                                                                 <  2:59294/239‑1 (MQ=255)          gACAGACCGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATgtcg                                                                                                                                                                                                                           <  1:256445/251‑1 (MQ=255)                ccGCTGCCACAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCAttt                                                                                                                                                                                                                     >  1:265195/1‑251 (MQ=255)                       cacAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTCACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAAATCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGGCATTTGAttccc                                                                                                                                                                                                              >  1:273004/3‑248 (MQ=255)                                 ggtggtGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAAc                                                                                                                                                                                                       <  2:86351/248‑1 (MQ=255)                                 ggtggcGTTTGGTGAAGTCGGGCTGGCAGGGGCGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGCGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGc                                                                                                                                                                                                    >  2:244186/1‑251 (MQ=255)                                         ttGGGGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAg                                                                                                                                                                                             <  1:288655/250‑1 (MQ=255)                                              gAAGTCGGGCTGGCAGTGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGACTTTCAGAAGCGGCGAAACACGGTTTTCGCCTGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTGAGCGTGTTCGACGACTTATAATGAGGAATACGGAGTGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACa                                                                                                                                                                                        <  1:313954/250‑1 (MQ=255)                                               aaGTCGGGCGGGCAGGGGAGTTCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAgg                                                                                                                                                                                      <  1:18085/251‑1 (MQ=255)                                                        gggCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCTAGCAACAGGGCTGCACGc                                                                                                                                                                             <  2:303867/250‑1 (MQ=255)                                                          gCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTa                                                                                                                                                                           <  1:4673/251‑1 (MQ=255)                                                            gggggAGATCCGCCCGGGTCCCAGGGGTCAGGAACGAATTTAAGAAGCGGCGAAACACTGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTAca                                                                                                                                                                         <  2:396638/250‑1 (MQ=255)                                                                          cGGTGCCCAGCGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATg                                                                                                                                                           <  1:325931/251‑1 (MQ=255)                                                                                 cAGCGGTCAGGAACGAATTTCAGAAGCAGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCgg                                                                                                                                                    <  2:348727/251‑1 (MQ=255)                                                                                    cGGTCAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTAt                                                                                                                                                 <  1:438479/251‑1 (MQ=255)                                                                                        cAGGAACGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGAc                                                                                                                                              <  1:1891/250‑1 (MQ=255)                                                                                             aCGAATTTCAGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAgg                                                                                                                                        <  1:490622/251‑1 (MQ=255)                                                                                                    tcagAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTAttt                                                                                                                                  <  1:424515/247‑1 (MQ=255)                                                                                                      agAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGCATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTa                                                                                                                                 >  1:42100/2‑249 (MQ=255)                                                                                                      agAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAg                                                                                                                               <  2:235353/250‑1 (MQ=255)                                                                                                       gAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGc                                                                                                                              >  1:384237/1‑251 (MQ=255)                                                                                                            ggcgAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGtt                                                                                                                         >  1:103859/1‑251 (MQ=255)                                                                                                              cgAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGcc                                                                                                                             >  1:360337/1‑245 (MQ=255)                                                                                                              cgAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGcc                                                                                                                             <  2:360337/245‑1 (MQ=255)                                                                                                                   cactgTTTTCGTCGGGAGATTGTTCCGGAGGCTCACGTACCGAAAACAGCGCCGGAAGTGATGCAGAATTGTGGCGTTAAAAAACTCTCCGACGCTCGTAGCGTGTTCGACGACTTATAAAGAGGAATCCGGAGGGAGATATGTCGTCATTTGATTACCGGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGATATTTAAGCCAGTTCCTGAAt                                                                                                                  <  2:109062/247‑1 (MQ=255)                                                                                                                    ccgggtttCGCCGGGGGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGTTATTTAAGCCAGTTACTGAATg                                                                                                                 <  2:277760/246‑1 (MQ=255)                                                                                                                    acGGTTTTCGCTGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGAGTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATg                                                                                                                 <  2:103859/251‑1 (MQ=255)                                                                                                                           tCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGTAAGGGATGCAGATTTTTGGCGGTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCaaaa                                                                                                           <  2:504742/250‑1 (MQ=255)                                                                                                                           tCGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAAt                                                                                                          >  1:90311/1‑251 (MQ=255)                                                                                                                            cGCCGGGCGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGCCCAAAGGGTATTTAAGCCAGTTACTGAATGCCACAATc                                                                                                         >  1:263884/1‑251 (MQ=255)                                                                                                                                   ggATTGTTCCGGCGGCTAACGTACAGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACGGAATGCCAAAATCAAAAGcc                                                                                                  <  2:395397/250‑1 (MQ=255)                                                                                                                                   cGATTGTTCCGGCGGCTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACGCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGcc                                                                                                  <  1:280862/251‑1 (MQ=255)                                                                                                                                                ggcTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCaaa                                                                                      >  1:450667/1‑250 (MQ=255)                                                                                                                                                  cTAACGTACCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCCAAATCAAAAGCCCAAGCGCGCAAAAGc                                                                                   >  1:288554/1‑251 (MQ=255)                                                                                                                                                      cgaaacGAAAAAAGCGCCGGAAGGGATCCAGTATTTTGACGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTGGACGACTTATAATGAGGAATACGGAGGTAGATATGTCGGCATTAGATTACCTGAAAACGGCCATCAAGCAACATGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGGATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCTCAGCGCGCAAAAGCTGGa                                                                               <  2:40284/246‑1 (MQ=255)                                                                                                                                                        taccGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGGATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAgg                                                                             >  2:154487/3‑251 (MQ=255)                                                                                                                                                            gAAAAAAGCGCCGTAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGGCGGCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGtt                                                                         <  2:384237/251‑1 (MQ=255)                                                                                                                                                            gAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGCCTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTAAAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGtt                                                                         >  1:17014/1‑251 (MQ=255)                                                                                                                                                                          aaGGGATGCAGATTTTTGGCGTTAAAAAACGCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGGCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTgg                                                           <  2:304876/251‑1 (MQ=255)                                                                                                                                                                           aGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTggg                                                          >  1:356291/1‑251 (MQ=255)                                                                                                                                                                                   aGATTTTTGGCGTTAAACAACTCTCCGCAGCGCTAAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTCCCTGAAAACTGCCCTCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGCTGCCAGCGGTATGAGCAAAGGGTCTTTAAGCAAGTTACTGACTGCCAACATAAAAAGCCCCAGCGCGCAACCGCTGGAGGAGTTGCACCGTTTTTTGGGGCTTGAGt                                                  >  2:38958/1‑251 (MQ=255)                                                                                                                                                                                           ggCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCgg                                          >  2:379004/1‑251 (MQ=255)                                                                                                                                                                                               ttAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGCGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGAATTTACGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGc                                         >  2:192139/1‑248 (MQ=255)                                                                                                                                                                                                       ctctCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAAcc                              >  1:94924/1‑250 (MQ=255)                                                                                                                                                                                                                  gcgtAGCGTGTGCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCTCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGgtgtggt                   <  1:186178/248‑4 (MQ=255)                                                                                                                                                                                                                  gcTTAGCGTGTTCGACGCCTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGgtggggt                   >  2:437141/1‑247 (MQ=255)                                                                                                                                                                                                                     tAGCGTGGTCGACGACTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCATAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTAt                 <  2:265195/250‑1 (MQ=255)                                                                                                                                                                                                                          tgtTCGACGACTTATAAGGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGAAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACAATTGGTGTGGTATTCGGt            >  2:25249/1‑250 (MQ=255)                                                                                                                                                                                                                              cgacgaCTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGtt       >  2:273660/1‑251 (MQ=255)                                                                                                                                                                                                                                  gacTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCCTTGGTGTGGTATTCGGTAAGTTCTAc   >  1:328638/1‑251 (MQ=255)                                                                                                                                                                                                                                   acTTATAATGAGGAATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTGGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTAcc  >  2:235003/1‑251 (MQ=255)                                                                                                                                                                                                                                               ||                                                                                                                                                                                                                                              AGCCTGCGCGACAGACCGCTGCCGCAGGATCTGGTGGTGTTTGGTGAAGTCGGGCTGGCAGGGGAGATCCGCCCGGTGCCCAGCGGTCAGGAACGAATCTCTGAAGCGGCGAAACACGGTTTTCGCCGGGCGATTGTTCCGGCGGCTAACGTGCCGAAAAAAGCGCCGGAAGGGATGCAGATTTTTGGCGTTAAAAAACTCTCCGACGCGCTTAGCGTGTTCGACGACTTATAATGAGATATACGGAGGGAGATATGTCGTCATTTGATTACCTGAAAACTGCCATCAAGCAACAGGGCTGCACGCTACAGCAGGTAGCTGATGCCAGCGGTATGACCAAAGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACC  >  NC_000913/4627061‑4627537

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