Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	89,239	G→A	intergenic (+207/‑395)	cra → / → mraZ	DNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	89,239	0	G	A	97.6%	123.3 / ‑4.2	41	intergenic (+207/‑395)	cra/mraZ	DNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (21/19);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

CCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGG‑ACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTTCGCACGTTGG  >  NC_000913/89008‑89481                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    ccgtggcgTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTTACTATTATTGATAAGTTCTTTTAAACCCTCCCGTTAATTAACTCACCAGCTTAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTa                                                                                                                                                                                                                                  <  1:388329/247‑1 (MQ=255)      gcgTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTaacaag                                                                                                                                                                                                                             <  1:577591/251‑2 (MQ=255)            ccAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGTGCGTAGAGTACGGTAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACttt                                                                                                                                                                                                                       <  2:107407/250‑1 (MQ=255)              aGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTc                                                                                                                                                                                                                     >  1:418095/1‑251 (MQ=255)                  gTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGTATTTCCCGCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTg                                                                                                                                                                                                                 <  2:376287/251‑1 (MQ=255)                     gCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCtt                                                                                                                                                                                                              >  1:299418/1‑251 (MQ=255)                                      aaTACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGTCTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTgg                                                                                                                                                                                             <  1:373756/251‑1 (MQ=255)                                         aCGCGCCAGGTGAATTCCCCTCTTGCTCGTAGAGTACGGGAAGTGACATCAATATCCTTAAAGTAAAAAAATCTATTCCTTACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCAt                                                                                                                                                                                          <  2:138848/251‑1 (MQ=255)                                                    gAATTTCCCTCTGGCGGGTAGAGTGCGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTAGTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTGTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCcaca                                                                                                                                                                               <  2:139863/251‑1 (MQ=255)                                                            ctctGGCGCGTAGAGTACGGGAATGGACCTCCATATCCTTAAAGTAAATAAGACTATTCCTGACTCTTATTAATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCACGCGCCACACAGACgc                                                                                                                                                                        >  2:332361/1‑249 (MQ=255)                                                                     gTAGAGTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGG‑ACt                                                                                                                                                             >  1:464507/1‑251 (MQ=255)                                                                         agTACGGGAATGGACATCAATATCCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGG‑ACTGtt                                                                                                                                                          >  1:363228/1‑250 (MQ=255)                                                                          gTACGGGACGTGACAACAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATGAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGG‑ACTGTTaa                                                                                                                                                        <  2:422246/251‑1 (MQ=255)                                                                                      acctCAATATCCTTAAAGTAAATAAGACTATTCCTGACTCTTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGCTTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCTCTGGCATCAAGCGCCACACAGACGCAACAGGG‑ACTGTTAACCGGGGa                                                                                                                                                 >  2:78999/4‑246 (MQ=255)                                                                                               tccTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAAccc                                                                                                                                                                                                  <  2:655908/187‑1 (MQ=255)                                                                                               tccTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAAccc                                                                                                                                                                                                  >  1:655908/3‑189 (MQ=255)                                                                                                 cTTAAAGTAAATAAGACTATTCCTGAATATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGCATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGAAGCAACAGGG‑ACTGTTAACCGGGGAAGATATGTCCTAAAAt                                                                                                                                 >  2:82690/1‑251 (MQ=255)                                                                                                           aTAAGACTATTCCTGTATATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGGTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCGTTGCTTTTGCGCTAACCCTCTGGCATCAAGCGCCACACAGACGCAACAGGG‑ACTGTTAACCGGGGAAGATATGTCCTAAAAGGCCGCTcgcg                                                                                                                       <  2:661336/251‑1 (MQ=255)                                                                                                              aGACTATTCCTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTACCAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACCGGG‑AATGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCg                                                                                                                    >  2:669736/1‑251 (MQ=255)                                                                                                                      ccTGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTACTTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTCCATGCAGAGATGAATTGCCGATTTAACAAGAACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACG‑CACAAGGGACTGTTAACCGGGGAAGATATGTCCTAAAATGACGCTCGCGTCgaagactg                                                                                                             >  2:359670/1‑243 (MQ=255)                                                                                                                      ccTGACTATTATTGATAAGTTCTTTTAAACACGCCCGTTAATTAACTCCCCAGCTGAAATTCACAATAATTAAGTGATATCGACAGAGAGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGG‑ACTGTTAACCGGGGAAGATATGACCTAAAATGCCGCTCGCGTCGCAAACTg                                                                                                             >  1:596835/1‑250 (MQ=255)                                                                                                                        tGACTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGACACACAGACGAAACAGGG‑ACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGac                                                                                                           >  1:605733/1‑250 (MQ=255)                                                                                                                          aCTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGG‑ACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACt                                                                                                        <  1:134990/251‑1 (MQ=255)                                                                                                                          aCTATTATTGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCACCAGCGGAAATTCACAATAATTAAGTGATATCGACAGCTCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGG‑ACTGTTAACCGGGGAAGATATGTCCGAAAATGCCGCTCGCGTCGCAAACTGACACt                                                                                                        <  2:462260/251‑1 (MQ=255)                                                                                                                                  tGATAAGTTCTTTTAAACCCGCCCGTTAATTAACTCAACAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGG‑ACTGTTAACCGGGGAAGATATGTCCTAAACTGCCGCTCGCGTCGCAAACTGACACTTTAAAttt                                                                                                >  2:328901/1‑251 (MQ=255)                                                                                                                                                 aaCCCGGCCGTTAATTAACCAACCAGCTGAAATTCACAATAATTAAGTGATATCGCCAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGG‑ACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGt                                                                                 <  2:194244/251‑1 (MQ=255)                                                                                                                                                    ccgcccgTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGG‑ACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAg                                                                              <  1:461080/251‑1 (MQ=255)                                                                                                                                                                           cTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGG‑ACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACg                                                        >  1:161124/1‑250 (MQ=255)                                                                                                                                                                               aaTTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGG‑ACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAg                                                   <  1:61732/251‑1 (MQ=255)                                                                                                                                                                                   cacaATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGG‑ACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGAt                                               >  1:282088/1‑251 (MQ=255)                                                                                                                                                                                      aataatTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGG‑ACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGTTGAGGGAtttt                                            >  2:360482/1‑251 (MQ=255)                                                                                                                                                                                       ataatTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGG‑ACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGAttttt                                           <  2:161124/251‑1 (MQ=255)                                                                                                                                                                                        taatTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACCAGCACTTTTCTTTGCTTTTCCGCAAACCCGCTGGCATCAAGCGACACACAGACGCAACAGGG‑ACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTCAAACGGAGATGACGATGAGGGAtttttt                                          >  1:677646/1‑251 (MQ=255)                                                                                                                                                                                               gTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGG‑ACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTg                                                                                                             <  2:431518/177‑1 (MQ=255)                                                                                                                                                                                               gTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGG‑ACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTg                                                                                                             >  1:431518/1‑177 (MQ=255)                                                                                                                                                                                                    atCGACAGCTCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCCTTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGG‑ACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTAt                               <  2:418095/250‑1 (MQ=255)                                                                                                                                                                                                       gACAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGG‑ACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCAt                           >  1:367687/1‑251 (MQ=255)                                                                                                                                                                                                         cAGCGCGTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGG‑ACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATaa                         <  2:510380/251‑1 (MQ=255)                                                                                                                                                                                                              cgTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGG‑ACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGTTGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGtt                     >  2:432902/1‑250 (MQ=255)                                                                                                                                                                                                               gTTTTTGCATTATTTTGTTACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGAAACAGGG‑ACTGTTAACCGGGGAAGATATGTCCTAAAATGCAGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTa                    >  1:336939/1‑250 (MQ=255)                                                                                                                                                                                                                                 tACATGCAGAGATGAATTGCCGATTTAACAAGCACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGCAACAGGG‑ACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTgg  <  2:683494/250‑1 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    CCGCGGCGTGCTCAGCCGTAGCTAAGCCGCGAACAAAAATACGCGCCAGGTGAATTTCCCTCTGGCGCGTAGAGTACGGGACTGGACATCAATATGCTTAAAGTAAATAAGACTATTCCTGACTATTATTGATAAATGCTTTTAAACCCGCCCGTTAATTAACTCACCAGCTGAAATTCACAATAATTAAGTGATATCGACAGCGCGTTTTTGCATTATTTTGTTACATGCGGCGATGAATTGCCGATTTAACAAACACTTTTCTTTGCTTTTGCGCAAACCCGCTGGCATCAAGCGCCACACAGACGTAACAAGG‑ACTGTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTTCGCACGTTGG  >  NC_000913/89008‑89481

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: ‑