Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,624,587	G→A	S68L (TCA→TTA)	ytjB ←	putative inner membrane protein Smp

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,624,587	0	G	A	97.6%	106.4 / NA	43	S68L (TCA→TTA)	ytjB	putative inner membrane protein Smp
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (25/16);  minor base C (1/0);  total (27/16)
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.82e-01

TATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAG‑TCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAAT  >  NC_000913/4624342‑4624819                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          tATGGGTGTAGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGCTAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAG‑TCg                                                                                                                                                                                                                                     >  2:346331/1‑251 (MQ=255)  aTGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTCACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAG‑TCgg                                                                                                                                                                                                                                    >  1:376787/1‑251 (MQ=255)    gggTGTCAAGTGTCAGGGACAGATAGCCGAACGGTCCTTTTTTACCCGCAATTGGCAAGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGGCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTGCATCCGTTAACTGATCGAGGATCGACTGAATGCGTTTTTCATCCGGTAAG‑TCGGTa                                                                                                                                                                                                                                  <  2:431482/251‑1 (MQ=255)          cGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATGGGCTCGACAATCTGCTGGTTAAAATAGCCGCCAGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTAGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATCCACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGATTTTCATCCGGTAAG‑TCGGTACGCAt                                                                                                                                                                                                                             >  2:170096/1‑250 (MQ=255)                        aGATAGCCGGGCGCTCCGTTTTTCCCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAG‑TCGGTACGCATCAGCGTTGCAACGtg                                                                                                                                                                                                              <  2:230589/251‑2 (MQ=255)                         gATAGTCGGGCGTTCCGTTTTCCCCCGCAATGGGCCCGACAATCTGCTGTTGAAAATAGCCGCCTGCTTTTTTACCGGCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACACGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAG‑TCGGTACGCATCAGCGGTGCAACGTTc                                                                                                                                                                                                             <  2:235085/251‑1 (MQ=255)                               cGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGACGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTAGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGAACGCCTGAATGCGTTTTTCATCCGGTAAG‑TCGGTACGCATCAGCGGTGCAACGTTCAGCGTc                                                                                                                                                                                                       >  1:180480/1‑251 (MQ=255)                                gAGCGGGCCGTTTTTACCCGCAATTGGCTGGACATGCTGCGGGGTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTGCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGATGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAAGGCGTTTTTGATCCGGTAAG‑TCGGTAGGCGTCAGCGTTTCAACGTGCAGCGTCa                                                                                                                                                                                                      <  2:453015/251‑1 (MQ=255)                                  gCGGTCCGTTTTTACCCGCAATTGGCTCGGCCATCTGCTGGGTAAAATAGCCGACTGCTTTATTACCGTCCTGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGGTACGGCTTTCCTCCGTTAACTGCTCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAG‑TCGGTAGGCATCAGCGGTGCCACGTTCTGCGTCAcc                                                                                                                                                                                                    >  2:179380/1‑251 (MQ=255)                                      tCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAG‑ACGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGAc                                                                                                                                                                                                >  1:20059/1‑251 (MQ=255)                                             ttACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGAGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTAATACACACCCGCGTCGAGGATACGGCTTTCATCCTTTAACTCATCGAGGATCGCCTCAATGCGTTTCTCATCCGGCAAG‑TCGGAACGCACCAGCGGTGACACGTTCAGCGTCCCCttgcgcacc                                                                                                                                                                                           >  2:236717/1‑241 (MQ=255)                                                 ccGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAG‑TCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTg                                                                                                                                                                                     <  1:346331/251‑1 (MQ=255)                                                       ttGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAG‑TCGGTACGCATCAGCGGTGCAACGTTCAGCGTCCCCTGACGCGCCAGGGTGCGGGcc                                                                                                                                                                               >  1:97866/1‑251 (MQ=255)                                                                  atcgtttgtgtAAAATCGCCTCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGTCGCTTGCGCCAGAACGTGCTATAAGATCGTCGTGTTCTTCATACACACCCGGGTCGAGTATACGGCTTGCATCCGTTAACTGAGCGAGGATCGCCGGAATTCGTTTGTCATCCGGTAAG‑TCGGTACGCATCAGCGGTGCATCGTTCAGCGTCACCGGACTCGCCAGGGTGCTGGCCATTTCTTCCAg                                                                                                                                                                    <  2:294397/241‑1 (MQ=255)                                                                         ggTTAAAATCGCCGCCTGTTTTTTTACCGTCGACCGCCAGCCGGGCGCGCACTTCGAAGCGTTCGCCAGACCGTGCGATAAGATCGCCTTGTTCGTCATACACACCCTCGTCGAGGATACTGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAG‑TCGGTACGCATCAGCGGTGCAACGTTCAGCGTCATCTGACGCGCCAGGGGGCGTGCGAGTTCTTCCAGCtgt                                                                                                                                                                <  2:295078/248‑1 (MQ=255)                                                                          gTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAG‑TCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGAt                                                                                                                                                            >  1:382535/1‑251 (MQ=255)                                                                                   gccgccTGCTTTTTTACCGTCGAGCGCCAGCAGATCGCGCACTTCGACGCTTTCGCCAAAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAG‑TCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTc                                                                                                                                                   >  1:410251/1‑251 (MQ=255)                                                                                                ttACCGTAGAGCACCAGCCGCACGCGCAATTCGACTATTTCGCCAGAACCTGCGTTACACTAGCCTTGTTCGTCATACACACCCGAGTCTAGGATACGGCTTTTCTCCGTTAACTGATCGAGGATCGCCTGAATGCGTGTTTCATCCGGTAAG‑TCGGTACGCATCAGCGTCGCAACGGTGCGCGTCACCTGACGCGCCAGGGTGCGGGCCCGTTCTTCCAGCCGTGGATTATAATGTCGCTGGTGGtttgt                                                                                                                                      >  2:386099/1‑249 (MQ=255)                                                                                                    cGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGCTAAGCTCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAG‑TCGGTAAGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACt                                                                                                                                   >  1:172662/1‑250 (MQ=255)                                                                                                     gTCGAGCGCCAGCCGGTCGAGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCAGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGAGTTTTTCATCCGGTAAG‑TCGGTACGCATAAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTa                                                                                                                                  >  1:356583/1‑250 (MQ=255)                                                                                                          gcgcCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAG‑TCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCAt                                                                                                                            <  1:8234/251‑1 (MQ=255)                                                                                                                                     tttCGCCAGAACGTGCGATACGATCGCCTTGTTCGTCATACACTCACGCGTCGAGGATACGGCTTTCATCCGTTATCTGATAGAGGATCGACTGAATGCGTTTTTCATCCGGTAAG‑TCGGTACGCATCAGCGGTGCAACGTTCAGCGTCAACTGACGCGCCAGGGAGCGGGCCAGTTCTTCCAGCTGTGGATTACGATGTCGCTGGTGGTTTTGCCTAAACCATGACGCTCCCTGCATCCGCGCCACTAAc                                                                                                 >  2:443339/1‑251 (MQ=255)                                                                                                                                                  tGCGATAAGATCGCATTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCATGAATGCGTTTTTCATCCGGTAAG‑TCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGAc                                                                                       >  1:132504/1‑248 (MQ=255)                                                                                                                                                     gATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAG‑TCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAAc                                                                                  >  1:114007/1‑250 (MQ=255)                                                                                                                                                          gATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATACGTTCACTGCTCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAG‑TAGGTACGCATCAGCGGTGCCACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTAGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTACCACGGCAAGACAGAACAGGAc                                                                             >  1:226725/1‑250 (MQ=255)                                                                                                                                                                cTTGTTAGTCATACACACCCGCGTCGATAATACGGCTTTCATCCGTTAACTGATCGAGGAACACCTGACTGCGTTTTTCATCCGGTAAG‑TCGGTACGCCTCAGCGGAGCAACGCTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTATTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACAATGAAGCTCACTGCATCAGCGCCACTAACAAGGCAAGACAGAACCGGACAATCACt                                                                      >  2:264117/1‑251 (MQ=255)                                                                                                                                                                     tCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAG‑GCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCg                                                                 >  2:7856/1‑251 (MQ=255)                                                                                                                                                                          aTACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAG‑TCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGc                                                             >  2:346524/1‑250 (MQ=255)                                                                                                                                                                             cacatCCGCGTCGAGGATCCGGCTTTCATCCGTTAACTGATCGGGGATCGCCTGAATGCGTTTTTCAGCCGGTAAA‑TCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTGTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGcc                                                         <  2:355297/251‑1 (MQ=255)                                                                                                                                                                                acCCGCGTCGTGGTTACGGCTTTCTGCCGTTAGCGGATCGAGGATCGCCTGAATGCGTTTTGCATCCGTTAAG‑TCGGGACGCATCAGCGGTGCAACTTTCAGCGTCACCTGACGCGCCAGGGTGCGTGCCAGTTCTTCCAGCTGTTGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAATGCAAGACAGAACAGGACAATCACTGCCCGATTCAGCCGGa                                                      <  2:436025/251‑1 (MQ=255)                                                                                                                                                                                  ccGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAG‑TCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAAt                                                    <  1:115077/251‑1 (MQ=255)                                                                                                                                                                                         gAGGATACGGCTTTCAGCCGTTAACTGAGCGAGGATCGCCTGAATGCGTTTTTCAGCCGGTAAG‑TCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTTACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAGGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGtt                                             <  1:179380/251‑1 (MQ=255)                                                                                                                                                                                           ggATTCGGCTTTCATCCGTTAACTTATCGAGGATCGCCTGAATGCGTTTTGCATCCGGTAAG‑GCGGGACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTCCGCTGTCGCTGGTGGTTTTGACTAACCCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGtttt                                           <  2:317845/251‑1 (MQ=255)                                                                                                                                                                                                         tCCGTTAACTAATCGCGGATCGCCTGAATGCGTTTTTCATCCGGTCAG‑TCGGTACGCATCAGAGGTGCAACGTTCAGCGTCAACTGACGCGCCAGGGAGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGAATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATACCTGCCCGATGCAGCCGCAATTTCAGTTTTGTGCGTGCCATCtt                             >  2:273424/1‑251 (MQ=255)                                                                                                                                                                                                                 ctaaTCGAGAATCGCCTGAATGCGTTTTTCATCCGGT‑AGATCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCAGTGGATTACGCTGTCGCTGGTGGTTTTGAATAAACCATGACGCACCCTGCATCAGAGCCACTAACAAGGCAAGACAGAACAGAACAATAACTGCCCGAAGAAGCCGGAATTTCAGTGTTGTGCTAGCCGTCTTCAACCCtt                     >  2:316086/4‑251 (MQ=255)                                                                                                                                                                                                                   gATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTAAG‑TCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCtt                     <  1:23158/249‑1 (MQ=255)                                                                                                                                                                                                                   gATCGAGGATCGCCGGAATGCGTTTTTCATCCGGTAAG‑TCGGTACGCAGCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTGGTGCGAGCCATCTTCCACCCTTTg                   <  1:175255/251‑1 (MQ=255)                                                                                                                                                                                                                     tagAGGATCGCCTGAATGCGTTTTTCATCCGGTAAG‑TCGGTACGCATCAGCGGAGCAACGTTCAGCGTAACCTGACGCGGCAGGGAGCGGGCCAGTTCTTCCAGATGTCGACTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGACTCAGAGCAACTAACAAGGCAAGACAGAACAGGCCAATCACTGCCAGATGCAGCCGGAATTTCAGTTTTGTGCGCGACATCTTCCACCCttgga                  >  1:260652/3‑247 (MQ=255)                                                                                                                                                                                                                         ggATCGCCTGAATGCGTTTTTCATCCGGTAAG‑TCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAAtt             >  1:422826/1‑251 (MQ=255)                                                                                                                                                                                                                                 tGAATGCGTTTTTCATCCGGTAAG‑TCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACtt     >  2:4366/1‑251 (MQ=255)                                                                                                                                                                                                                                  gAATGCGTTTTTCATCCGGTAAG‑TCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACtt     >  1:6444/1‑250 (MQ=255)                                                                                                                                                                                                                                  gAATGCATTTTTCCTCCGGTAAG‑TCGGTACGCATCAGCGGTTCAACGTTCAGCTTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGTCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTa    <  2:422447/251‑1 (MQ=255)                                                                                                                                                                                                                                    aTGCGTTTTTCATCCGGTAAG‑TCGGTACGCATCAGAGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAAt  >  1:161703/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                          TATGGGTGTCGAGTGTCAGGCGCAGATAGCCGAGCGGTCCGTTTTTACCCGCAATTGGCTCGACAATCTGCTGGTTAAAATAGCCGCCTGCTTTTTTACCGTCGAGCGCCAGCCGGTCGCGCACTTCGACGCTTTCGCCAGAACGTGCGATAAGATCGCCTTGTTCGTCATACACACCCGCGTCGAGGATACGGCTTTCATCCGTTAACTGATCGAGGATCGCCTGAATGCGTTTTTCATCCGGTGAG‑TCGGTACGCATCAGCGGTGCAACGTTCAGCGTCACCTGACGCGCCAGGGTGCGGGCCAGTTCTTCCAGCTGTGGATTACGCTGTCGCTGGTGGTTTTGACTAAACCATGACGCTCCCTGCATCAGCGCCACTAACAAGGCAAGACAGAACAGGACAATCACTGCCCGATGCAGCCGGAATTTCAGTTTTGTGCGAGCCATCTTCCACCCTTTGAAAATTTGAGACTTAAT  >  NC_000913/4624342‑4624819

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