Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,445,981	T→G	29.3%	G34G (GGT→GGG)	ydbL →	DUF1318 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,445,981	0	T	G	29.3%	67.4 / 4.3	48	G34G (GGT→GGG)	ydbL	DUF1318 family protein
Reads supporting (aligned to +/- strand):  ref base T (10/21);  new base G (14/2);  total (25/23)
Fisher's exact test for biased strand distribution p-value = 5.09e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCACCTAAGGAACCGATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTCGGGCCAAACCGGGGGAATATGTACAAGGGATA  >  NC_000913/1445756‑1446179                                                                                                                                                                                                                                  |                                                                                                                                                                                                         gCTGCACCTAAGGAACCGATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGTTGAAGCCAGAACTCCGGGGCGGGAGGGTgga                                                                                                                                                                                                       >  1:618514/1‑227 (MQ=255)               aCCGATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGGGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGATTCTGGTTGAAGCCAGAACTCCGGGGCGGGGAGGGGAAAAA‑TTTTACGGttt                                                                                                                                                                                        >  2:694790/1‑228 (MQ=255)                ccGATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCATGGGCGGGTAGGTGAAAAA‑TTTTTCGGttgt                                                                                                                                                                                       >  1:744834/1‑227 (MQ=255)                 cGATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCGGCAATGTAAGGGAATTGACTTTGGTTGAAGCCAGAACTCAGGGGCGGGAAGGGGAAAAA‑tttttagggaaat                                                                                                                                                                                      >  2:338942/1‑220 (MQ=255)                  gATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGGGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTATGGCTGAAGCCATAAATCAGGGGGGGGGAGGGGAAACA‑ATTTAAGGTTaaag                                                                                                                                                                                     >  2:89352/1‑226 (MQ=255)                   aTCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGCTCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGCCTCTGGATGAAGCCCGAACTCCGGGGCGGGAAGGGGAAAAA‑TTTATAGGTTTTCTg                                                                                                                                                                                    >  1:270388/1‑229 (MQ=255)                   aTCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGGTGAACA‑TTTTCGGGTTATCTg                                                                                                                                                                                    >  1:422979/1‑229 (MQ=255)                                  aatgTTAAAATTGAGCATGAAGTCATCAACAAGGCAGGCAAAGATGGCGAAGAGCTGTTTGAAACTCGTAGCGATTTTTTCTGAGGTTATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTGGTTGCGCTGAAAACg                                                                                                                                                                     <  1:209378/226‑1 (MQ=255)                                  aaCGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGTGGAAGCCAGAACTCAGGGGCGGGAAGGGGAAAAA‑TTTTAAGGTTAATTGGTTTGGCTGAAAACg                                                                                                                                                                     >  1:296512/1‑229 (MQ=255)                                  aaCGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGGGAAACA‑TTTTACGGTTAATTGGTTGGGCTGAAAACg                                                                                                                                                                     >  1:411200/1‑229 (MQ=255)                                  aaCGTTAAAATTGAGCATGAGATCATCACCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATGACTTCTTTGTGCGTTTCTTGGTGGGCTGGTAAGCAGAAATGTAATGGAATTGACTCTGGATGAAGCCAGAACTCAGGGTCGGGTATGTCAACAA‑TTTTACGGATATCTGGGTgctggcaaatcg                                                                                                                                                                     >  1:349735/1‑219 (MQ=255)                                    cGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGAGGGGAAAAA‑TTTTTTGGTGATTCTGTGTCGCTGGAAACggg                                                                                                                                                                   >  2:413740/1‑228 (MQ=255)                                      ttAAAATTGATCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTGGTTGCGCTGAAAACGGATg                                                                                                                                                                 <  1:694791/229‑1 (MQ=255)                                        aaaaTTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGAGG‑GTAACACTTTTACGGTTATTTGGTTGGCGTGGAAACGGACGCt                                                                                                                                                               >  1:290143/1‑229 (MQ=255)                                                     agaTCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAAt                                                                                                                                                  <  1:448573/229‑1 (MQ=255)                                                             atcaAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGcc                                                                                                                                          <  1:459956/229‑1 (MQ=255)                                                                 aGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCCGGGGCGGGGAGGGGAAACA‑TTTTACGGGTATCTTGTTGCGCTTGAAACGGATGCTGAAACAGAGAAATTATTAGCCGACa                                                                                                                                      >  1:469448/1‑229 (MQ=255)                                                                 aGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGGGGAACA‑TTTTACGGGTATCTGGTTGCGCTGGAAACGGGGGCGGGAAAAGGGAAATTATTAGCCGCCa                                                                                                                                      >  1:306763/1‑229 (MQ=255)                                                                          aggATGTCGAATGGCTGTTTGAATCTCGTAGCGATTTTTTCTGAGGTGATGAGGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCg                                                                                                                             <  2:469448/227‑1 (MQ=255)                                                                          aaGATGGCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCg                                                                                                                             <  1:568444/229‑1 (MQ=255)                                                                                    agagGTGATTGAAACTCGTGGCGATCTGTTCTGGGGGGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGc                                                                                                                   <  2:270388/229‑1 (MQ=255)                                                                                    agagCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGc                                                                                                                   <  1:89352/229‑1 (MQ=255)                                                                                          gctTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGACGAAAACATTACTTCTTTGTGCGTTTCTTGTTGTGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGCACTCAGGGGCGGGTGAGTGCAAAA‑TTTTACGGTTATCTTGTGTCGCTGAAAAAGGGTCGTTGAACTGAGAAACTCATCGCCCACTTTAAAGCCTAACATaaaacagagga                                                                                                             >  2:342836/1‑221 (MQ=255)                                                                                             tGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTTCGGTGATCTGGATGCGCTGGAAAAGGTGGCTGAAACAGAGCAATTATTGACCGACATTAATACCGAAAGAAAAGAGAGTTAcaa                                                                                                          >  2:520783/1‑227 (MQ=255)                                                                                                       tgcggttttttttGGGGTGATGTAGAAGAAAACATTACTTGTTTGTGCGTTTCTTGGTGGGAGGGTAAGGGGAAATGCAAGGGCATTGACTTTGGTTGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTGGTTGCGATGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGgc‑a                                                                                               <  2:290143/222‑1 (MQ=255)                                                                                                       tgcggttttttttGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTGGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGgc‑a                                                                                               <  2:306763/222‑1 (MQ=255)                                                                                                                    gcggTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGCAAGTGAAGAA‑TTTTCCAGTGATCTGGTTGCGCTGAAAACGGAGATCGCACCAGAGAAATTAGCAACCGCCATTATTGCCGAAACAAAAGCGCAGTACCAACAACTGGA‑CACGCAAaataaaa                                                                                  >  2:41751/3‑227 (MQ=255)                                                                                                                       gTGATTATGAAGAAAACATTCCTTCGTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTGTACGGTTATCTGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATgtgt                                                                               <  1:302905/229‑1 (MQ=255)                                                                                                                         gatgatgaAGAAAGTATTACTTCTTTGTGCGGTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTATGGATGAAGCCAGAACTCAGGGGCAGGTAGGTGAAACA‑TTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCg                                                                             <  1:345058/229‑1 (MQ=255)                                                                                                                          atgatgaAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTGGTTGCGCTGAAAACg                                                                                                                                                                     <  1:717270/141‑1 (MQ=255)                                                                                                                          atgatgaAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGGGGAACA‑TTTTTCGGGTATCTGGTTTCGCGGGAAACGGATGCGTAAACAAATCAATGAGGAGCCGACATTAATGCCCAACGGAAAGAGAGTGAACAAAAACAGGcgaaagcaacaaaagttgtc                                                                              >  2:522893/1‑217 (MQ=255)                                                                                                                          atgatgaAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAAGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGCGTTACAAACAACTGGC‑AAAGCAAAATACTGTGACg                                                                             >  1:405016/1‑228 (MQ=255)                                                                                                                                  gaaAACTTTTCTTTTTTGTGGGTTTGTTGTTGGTCTGGTAAGGGGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGTAgatga                                                                     <  2:422979/228‑1 (MQ=255)                                                                                                                                          tACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGATAGGGGAAAAA‑TTTTAAGGTTATTTGGATGGGCTGCAAACTGATTCTGAAACCGCGGAATTTGTAGCCGGCATTAATGTCGAACGTACAGCGAGTTACCCACAACAGGG‑CAAGCCAAAACATGAGTCGGTGGCTGatatatcgca                                                            >  1:543917/1‑223 (MQ=255)                                                                                                                                                         ttCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTGGGTGCGCTGAAAACGGATGCTGGAACAGAGAAATTTGTAGCCGACATTAATGCCGCACGTAAAGCGGGTTTCCAAAAACAGGG‑ACAGCCAAAAAAAAGGGCGGTAGAGGGGATCCGCAAACACCCCGGGAaaaa                                             >  2:454557/1‑229 (MQ=255)                                                                                                                                                                   gCTGGTAAGGAGCAATGTAATGGCATTGACTTTGGATGAAGCCAGAACTGAGGGGCGGGTAGGTTGAACA‑TTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAGCGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACGCGCCGGGCAAAAGCTCGTAGCt                                   <  1:41751/229‑1 (MQ=255)                                                                                                                                                                            cagcaATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTTCGGTTATCTGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTATTAGCCGACATTAAAGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGGGGCGGTAAATGGAATCGCGAAACTCGCCGGGCCAAAGATCGGAGCTCGGGCCaaa                          >  1:523691/1‑229 (MQ=255)                                                                                                                                                                             agcaATGTAATGGCATTGACTTTGGATGAAGCCAGAACTCAGGGCCGGGTAGGTGAAACA‑TTTTACGGTTATCTGGATGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTCGGGCCaaa                          <  1:520783/228‑1 (MQ=255)                                                                                                                                                                              tgaaTGTAAGGGCATGGACTCTGGATTAAGGCAGAACTCGGGGGCGGGGAGGGGAAACA‑TTTTACGTTTAACTGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAACCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACCCGCCGGGCAAAAGCTCGTAGCTCGGGCCAAAc                         <  2:744833/226‑1 (MQ=255)                                                                                                                                                                                 gtGTAATGGAGTTGACTTTGGATGAACACAGAACTAAGGGGCGGGTTGGGGAAACA‑TTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGGCGACATTAATGCCGAACGTAAAGCGAGTTACCAAAAACTGGC‑AAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTCGGGCCAAACCggg                     <  1:338942/228‑1 (MQ=255)                                                                                                                                                                                 aTGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGAGGGGAAACA‑TTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGGAGATGATATCGCGAAACTCACCGGGCACAAAATCGTAGCTCGGGCCAAAcacgg                     >  1:423455/1‑225 (MQ=255)                                                                                                                                                                                    tAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGAGGTGAAACA‑TTTTACGGTTATCTGGGTGCGCTGAAAACGGATGCTGAAACAAAGAAATTATTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGGAGATGATATCGCGGAACTCGCCGGGCAAAAGACCGTAGCTCGGGCCAAACCGCGGGa                  >  1:119037/1‑229 (MQ=255)                                                                                                                                                                                            tttatttgggATGAAGCCAAACCGCAGTGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTGGTTGCGGTGAAAACGGATGCGGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTCGGGCCAAACCGGGGGAATATGTAc          <  2:184659/221‑1 (MQ=255)                                                                                                                                                                                            ttGCCTCTGGATGAAGCAAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTCGGGCCAAACCGGGGGAATATGTAc          <  2:622613/229‑1 (MQ=255)                                                                                                                                                                                            ttGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGGGAAACA‑TTTTACGGTTATCTGGGTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGCACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGGAGATGGTATCGCGGAACACGCCGGGCAAAAGCTCGAAGCTCGGGCCAAACCGAGGGAATATGtca          >  1:184659/1‑227 (MQ=255)                                                                                                                                                                                              gACTTTGGATGAAGCCTGAACCGAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTAGCTGGTTGCGCTGAAAACGGATGCTGAAACGGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTCGGGCCAAACCGGGGGAATATGTACaa        <  2:296512/229‑1 (MQ=255)                                                                                                                                                                                                     gATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTCGGGCCAAACCGGGGGAATATGTACAAGGGATa  <  1:413740/228‑1 (MQ=255)                                                                                                                                                                                                                            ggTAGGTGAAACA‑TTTTACGGTTATCTGGTTGCGCTGAAGACGGATCCTGAAACAGAGAAATTAGTGGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTCGGGCCAAACCggg                     <  1:405/186‑1 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                         GCTGCACCTAAGGAACCGATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTCGGGCCAAACCGGGGGAATATGTACAAGGGATA  >  NC_000913/1445756‑1446179

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