Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,628,560	T→G	intergenic (+13/+295)	nadR → / ← ettA	DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR/energy‑dependent translational throttle protein EttA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,628,560	0	T	G	89.7%	114.7 / ‑2.7	39	intergenic (+13/+295)	nadR/ettA	DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR/energy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base G (18/17);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGT  >  NC_000913/4628318‑4628807                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        atgaATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAg                                                                                                                                                                                                                                                 <  1:213959/251‑1 (MQ=255)  tgaATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGCGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAgg                                                                                                                                                                                                                                                >  1:485447/1‑251 (MQ=255)   gaATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGTTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATACCCGCGATGAAACGGCTCAAAGGc                                                                                                                                                                                                                                               <  1:424286/251‑1 (MQ=255)    aaTACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCg                                                                                                                                                                                                                                              >  1:193221/1‑251 (MQ=255)     aTACCGTTTCGATCTGGTGATCCGGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGAGGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGa                                                                                                                                                                                                                                             <  2:483017/251‑1 (MQ=255)          gTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGGGGtat                                                                                                                                                                                                                                        >  1:268221/1‑251 (MQ=255)          gTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGtat                                                                                                                                                                                                                                        <  2:202623/251‑1 (MQ=255)             tCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATaa                                                                                                                                                                                                                                      <  1:164482/250‑1 (MQ=255)              cGATCTGATGATCCTGCTGGAGACCAGCACGCCGTGGGAGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGAAAAGAGTTCCCGAACTTGCTGGTGGAGATGCTGGAAGAGAACCATATCGAATTTGTGCGGGATGAAGAGGCCGATTATGACAGCCGTTTCCTGCGCTTCGTGGAGCTGGTGCGGGCGATGCTGGGGGAGCCGAGCTAACCGCGATGACACGGCTCAAAGGCGCGGGATaaa                                                                                                                                                                                                                                     >  2:455330/1‑250 (MQ=255)                    ggTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGCGGATCGCAAAGAGGTCCCGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAAATAAACGCGATGAAACGGCTCAAGGGCGAGGTAAAAAATAAGtt                                                                                                                                                                                                                              >  2:462134/1‑251 (MQ=255)                    ggTGATCCTGCTGGAGAACAACACGCCGAGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTAGATCGCAAAGAGGTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGCGGAAGATTAAGACAGTCGTTTCCTGCGCAGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCCGGGGTAACCGTGATGAACCTGCTCAAAGGCggtggctaaaatgaggt                                                                                                                                                                                                                              >  2:201310/1‑235 (MQ=255)                    ggTGATCCTGCGGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGtt                                                                                                                                                                                                                              <  2:28320/251‑1 (MQ=255)                     gTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGttt                                                                                                                                                                                                                             <  1:109192/251‑1 (MQ=255)                        aTCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCTCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGtttttt                                                                                                                                                                                                                          <  1:180784/251‑1 (MQ=255)                              ctgGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGAAAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATTGGGGAGCAGAGATAACCGTGATGAAACTGCTCACAGGCGTGGTATCAAATGAGTTTTTTtgatg                                                                                                                                                                                                                     >  2:169341/1‑250 (MQ=255)                                            gCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGATGGAAGAGAACAATATCGCATTCGTGCGGGTTGCAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACctctct                                                                                                                                                                                                     >  1:457920/1‑251 (MQ=255)                                                  ggtggCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACGTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAaga                                                                                                                                                                                               <  2:206325/251‑1 (MQ=255)                                                                  cGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTaaa                                                                                                                                                                               <  1:290271/251‑1 (MQ=255)                                                                  cGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTaa                                                                                                                                                                                >  1:289420/1‑250 (MQ=255)                                                                    cAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAAc                                                                                                                                                                              >  1:83770/1‑250 (MQ=255)                                                                                  cGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGTGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGaaaaaa                                                                                                                                                                <  1:242479/250‑1 (MQ=255)                                                                                     tGGATCGCAAAGAGTTCCAGAACTTGATGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGGTGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAAACTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAACCCTg                                                                                                                                                            >  2:352352/1‑251 (MQ=255)                                                                                                       aGAACTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGCGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGCAAAAACCTGCACGGGTGGCCCGCTACg                                                                                                                                          >  2:339120/1‑251 (MQ=255)                                                                                                          aCTTGCTGGTGGAGATGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGaaa                                                                                                                                       <  1:295243/251‑1 (MQ=255)                                                                                                                     gagaTGCTGGAAGAGAAAAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGAtact                                                                                                                            <  2:257814/251‑4 (MQ=255)                                                                                                                        aTGCTGGAAGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTa                                                                                                                           >  1:34424/1‑249 (MQ=255)                                                                                                                               aaGAGAACAATATCGAATTCGTGCGGGTTGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCtgttg                                                                                                                  >  1:1834/1‑251 (MQ=255)                                                                                                                                                           tGAAGAGGACGATTATGACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCAcgccgtat                                                                                      >  1:497523/1‑248 (MQ=255)                                                                                                                                                                       ttatgACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATACCCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAg                                                                                                                         >  2:213344/5‑204 (MQ=255)                                                                                                                                                                       ttatgACAGCCGTTTCCTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAg                                                                                                                         <  1:213344/200‑1 (MQ=255)                                                                                                                                                                                      ccTGCGCTGCGTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATg                                                           >  1:47133/1‑251 (MQ=255)                                                                                                                                                                                        tgcgctgcgTGGAGCTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCg                                                         >  2:104251/1‑251 (MQ=255)                                                                                                                                                                                                                          ggtgCAGATATAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATTAGTTGAAATCCTCTCTGGAAGCGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGGGTCCCGCTACGAAATTGCGGATACTAATGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCtggagattca                       <  1:485711/248‑10 (MQ=255)                                                                                                                                                                                                                                 agataccctcgataaaaCTGCTCAAAGGCGGGGTATAAAAGAAGTTTTTTAG‑TGAGTGGAAACCCTCTCTGGAAGAGGCTGCCGAGATTAAACAAGGTTGGAAAAAACCGGCACGGGTGACCCGCTACTAAATTGCGGATACGAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTa                <  2:269064/237‑1 (MQ=255)                                                                                                                                                                                                                                 agaTAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCtggagattcacaacgt                 >  1:94665/1‑235 (MQ=255)                                                                                                                                                                                                                                   aTAACCGCGATGAAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGAAAGAGCCTGTCGATATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTAtt              <  1:294025/251‑1 (MQ=255)                                                                                                                                                                                                                                            atgaAACGGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCg                                                                                                           <  1:408308/149‑1 (MQ=255)                                                                                                                                                                                                                                            atgaAACGGCTCAAAGGCGAGGGATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGCAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCg                                                                                                           >  2:408308/1‑149 (MQ=255)                                                                                                                                                                                                                                               aaaCTGCTCAAAGGCGAGGTATAAAATAAGTTTTTT‑GAGGGGATGAAACCCTCTCTGGAAGAGGCTGTCGAAATTAAACAAGGTTTGAAAAACCCTGCAGGGGTGACCCGCTCCGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAAtttcttgg  <  2:190454/251‑7 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                        ATGAATACCGTTTCGATCTGGTGATCCTGCTGGAGAACAACACGCCGTGGGTGGCGGATGGTTTACGCAGCCTCGGCAGTTCGGTGGATCGCAAAGAGTTCCAGAACTTGCTGGTGGAGATGCTCGAAGAGAACAATATCGAATTTGTGCGGGTTGAAGAGGAAGATTACGACAGTCGTTTCCTGCGCTGCGTGGAACTGGTGCGGGAGATGATGGGGGAGCAGAGATAACCGTGATGAAACTGCTCAAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCTGGT  >  NC_000913/4628318‑4628807

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