Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,628,805	G→T	intergenic (+258/+50)	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,805	0	G	T	92.3%	110.7 / ‑1.8	39	intergenic (+258/+50)	nadR/ettA	DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR/energy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base T (23/13);  total (25/14)
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.83e-01

AAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCT‑GGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATC  >  NC_000913/4628565‑4629041                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             aaaGGCGTGGTAAAAAATAAGTTTTTT‑GATGAGTTGAAACCCTCTCGGGAAGAGTCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACAAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCTTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGc                                                                                                                                                                                                                                   <  1:128432/251‑1 (MQ=255)    ggCGAGGTATAAAATAAGTTTTTT‑GATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGTTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCgg                                                                                                                                                                                                                                <  1:412455/251‑1 (MQ=255)           tataAAATCAGTTATTT‑GATGAGTTGAAAACCTCTCTGGCAGAGGCTGTCAAGCTTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTACGGCTGATGTAGGCCGGAAAAGCCGTTCACGCCGTATCCGGCATTAGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGCGTCTCACAACTTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGc                                                                                                                                                                                                                         >  2:679700/1‑251 (MQ=255)             aaaaaagaagttttttGAGGGGTTGAAAACCGCTCTGGAAGAGGTGGTCGAGATTAACAAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATGGCGGATATTAAGGCTTTTGTGGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGc                                                                                                                                                                                                                         <  2:410576/240‑1 (MQ=255)                            ggatgatTTGAAACCCTCTCTGGAAGAGGTTGACGAGATAAAAAAAGGTTTAAAAAACCCTGCACGGGTGACCCGCTACGAAATGGCGGATACAAAGGTTGTTGTAGGCCGAAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                                                                                                                         <  1:619549/250‑1 (MQ=255)                                     gAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAAGGTATTGAAATTTCTTGGTTTT‑TAGGCCGGATAAAGCGTTCACGCCGCATCCGGCATTTTAc                                                                                                                                                                                                >  2:461615/1‑251 (MQ=255)                                          ccTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCAtt                                                                                                                                                                                           >  1:274487/1‑251 (MQ=255)                                                             cgataatAAACAATGTTTGAAAAAACCGTCACGGGTGCCCCTTTACGAATTTGCGGATACCAAGGCTGTTGTAGCCCGGAAAGGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACACCGTATTGAAATTTCT‑TGGTTTGTAGGCCGTATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGt                                                                                                                                                                         <  1:379647/244‑1 (MQ=255)                                                               agaTTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACt                                                                                                                                                                      >  2:33897/1‑251 (MQ=255)                                                                      aCAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCT‑GGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGAAATACGCTTGTACTTGATAcc                                                                                                                                                               >  1:600911/1‑250 (MQ=255)                                                                        aagttttGAAAAACCCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTc                                                                                                                                                                                                                     <  1:250383/191‑1 (MQ=255)                                                                                 aaaaaCCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTcc                                                                                                                                                     >  1:672002/1‑250 (MQ=255)                                                                                  aaaaCCAGCAAGGGTGACCCGCTACGAAATTGCGGATACTACGGCTGATGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTAAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCAAGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACACTTCGGCTcc                                                                                                                                                     >  2:449584/1‑249 (MQ=255)                                                                                   aaaCCTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAgt                                                                                                                                                  >  1:533766/1‑250 (MQ=255)                                                                                      ccTGCACGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGtgcg                                                                                                                                               >  1:412526/1‑251 (MQ=255)                                                                                            cGGGTGACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCCTTTTACGCATTACTTCCCAATACGCTTGTACTTGATAAGCTTCGGCTCCAGTTCGTCTgcc                                                                                                                                         >  2:498114/1‑250 (MQ=255)                                                                                                 gACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTTCGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCCGTGCGTCTGTGCCCac                                                                                                                                    >  2:212432/1‑250 (MQ=255)                                                                                                 gACCCGCTACGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGccc                                                                                                                                      >  1:177115/1‑249 (MQ=255)                                                                                                  acccgataCGAAATTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCCCGACGTATCCGGCATTCGGTGCTCAATGCCTGATGCGCCGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTCTTGAAATTTCT‑TGGTTTGTATGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTCCTTCGCAATACGCTTGTACTTTATACGCTTCGGCTCCCGTGCGTCTGCGCCCAgc                                                                                                                                   >  2:675468/8‑251 (MQ=255)                                                                                                         aCGAAATTGCGGATACTAAGGCTTTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAgcgtgcg                                                                                                                              <  1:33897/249‑1 (MQ=255)                                                                                                          cgcaaTTGCGGATACTAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCCGCTCCAGTGCGTCTGCGCCCAGCGTGGGttt                                                                                                                           >  2:589034/4‑251 (MQ=255)                                                                                                                         tAAGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGt                                                                                                             >  1:568799/1‑250 (MQ=255)                                                                                                                          aaGGCTGTTGTAGGCCGGAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGtcc                                                                                                           >  2:113171/1‑249 (MQ=255)                                                                                                                                        ccgaaaaaggcgtTCACGCCGTTTCCGGCATCCGGTGCTAAATGCCTGATGCGACGTTTACCGCGGCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTAcc                                                                                             <  2:572039/240‑1 (MQ=255)                                                                                                                                          ggAAAAGGCGTTCACGCCGTATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCtt                                                                                           >  1:115462/1‑251 (MQ=255)                                                                                                                                                               gccGGCATTCTGTGCTCAATGCCTGATGCGACGCTTACCTCGTCTTATCAGGCCTACAGTTTTACAACTTATTGAAATTGCT‑GGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTAcc                                                                        <  2:600911/248‑1 (MQ=255)                                                                                                                                                                     aTTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGAGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCt                                                                >  2:380122/1‑251 (MQ=255)                                                                                                                                                                       gcGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTg                                                               <  1:618194/249‑1 (MQ=255)                                                                                                                                                                                 aTGCCTGATGCGACGCTGACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAgga                                                    <  1:559680/251‑3 (MQ=255)                                                                                                                                                                                      tGATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGc                                               >  1:136786/1‑251 (MQ=255)                                                                                                                                                                                       gATGCGACGCTTACCGCGTCTTATCAGGCCTGGAGATTCACAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCg                                              >  1:182525/1‑251 (MQ=255)                                                                                                                                                                                                           ttATCAGGCCTGGAGATTCATAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCCCGCCGCATCCGGCATTTTTCGCATTACCTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTAACTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGaa                           <  2:177115/250‑1 (MQ=255)                                                                                                                                                                                                               caggcctggagattcaCAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAAc                      >  2:265543/16‑251 (MQ=255)                                                                                                                                                                                                                 ggcctggagattcaCAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACgg                    >  1:474189/14‑251 (MQ=255)                                                                                                                                                                                                                      ggagattcaCAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTg               >  1:315263/9‑251 (MQ=255)                                                                                                                                                                                                                               cAACGTATTGAAATTTCT‑TTGTTTGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTTTACTTGATACGCTTCGGCTCCAg                                                                                                                                                   <  2:101603/110‑1 (MQ=255)                                                                                                                                                                                                                               cAACGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAg                                                                                                                                                   >  1:101603/1‑110 (MQ=37)                                                                                                                                                                                                                                aaCGTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGAATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAAcc     <  1:652233/251‑1 (MQ=255)                                                                                                                                                                                                                                   gTATTGAAATTTCT‑TGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGTGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCAGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAGTCCAGGATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATc  >  1:581496/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                             AAAGGCGAGGTATAAAATGAGTTTTTTTGATGAGTTGAAAACCTCTCTGGAAGAGGCTGTCGAGATTAAACAAGGTTTGAAAAAACCTGCACGGGTGACCCGCCACGAAATTGAGGATGCTAAGGCTGTTGTAGACCGGAAACGGTGTTCACGCCGCATCCGGCATTCGGTGCTCAATGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACAGTTTTACAACGTATTGAAATTGCT‑GGTTTTGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCATTTTACGCATTACTTCGCAATACGCTTGTACTTGATACGCTTCGGCTCCAGCGCGTCTGCGCCCAGCGTGCGTTTCTTGTACTCTTCGTACTCGGTAAAGTTACCTTCGAAGAACTCAACTTTACCTTCATCCTGGTAATCCAGAATGTGCGTGGCGATACGGTCGAGGAACCAACGGTCGTGCGAGATAACCATC  >  NC_000913/4628565‑4629041

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