Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,445,981	T→G	27.2%	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	27.2%	70.8 / 3.8	51	G34G (GGT→GGG)	ydbL	DUF1318 family protein
Reads supporting (aligned to +/- strand):  ref base T (12/22);  new base G (14/2);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 7.17e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCACCTAAGGAACCGATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTG‑GTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑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‑TTTTACGGTTATCTG‑GTTGCGCTGAAAACg                                                                                                                                                                     <  1:209378/226‑1 (MQ=255)                                  aaCGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGTGGAAGCCAGAACTCAGGGGCGGGAAGGGGAAAAA‑TTTTAAGGTTAATTG‑GTTTGGCTGAAAACg                                                                                                                                                                     >  1:296512/1‑229 (MQ=255)                                  aaCGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGGGAAACA‑TTTTACGGTTAATTG‑GTTGGGCTGAAAACg                                                                                                                                                                     >  1:411200/1‑229 (MQ=255)                                  aaCGTTAAAATTGAGCATGAGATCATCACCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATGACTTCTTTGTGCGTTTCTTGGTGGGCTGGTAAGCAGAAATGTAATGGAATTGACTCTGGATGAAGCCAGAACTCAGGGTCGGGTATGTCAACAA‑TTTTACGGATATCTG‑GGTgctggcaaatcg                                                                                                                                                                     >  1:349735/1‑219 (MQ=255)                                    cgctAAAATTGAGCATGATCTCATCATCAAGGCAGACAAATATGTCGAAGAGATGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTTCGTTTCTTGTTGGGCTTGTAAGAAGCAATGTAATCGCATTGCCTCTGCAGGATGCCAGAATTCAGGGACGGCTTGGTTAAACC‑TTTTTCGTGTTTCTg‑gcatccatgcacacgac                                                                                                                                                                   >  1:344196/4‑213 (MQ=255)                                    cGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGAGGGGAAAAA‑TTTTTTGGTGATTCT‑GTGTCGCTGGAAACggg                                                                                                                                                                   >  2:413740/1‑228 (MQ=255)                                      ttAAAATTGATCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTG‑GTTGCGCTGAAAACGGATg                                                                                                                                                                 <  1:694791/229‑1 (MQ=255)                                        aaaaTTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGAGG‑GTAACACTTTTACGGTTATTTG‑GTTGGCGTGGAAACGGACGCt                                                                                                                                                               >  1:290143/1‑229 (MQ=255)                                                     agaTCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTG‑GTTGCGCTGAAAACGGATGCTGAAACAGAGAAAt                                                                                                                                                  <  1:448573/229‑1 (MQ=255)                                                             atcaAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTG‑GTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGcc                                                                                                                                          <  1:459956/229‑1 (MQ=255)                                                                 aGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCCGGGGCGGGGAGGGGAAACA‑TTTTACGGGTATCTT‑GTTGCGCTTGAAACGGATGCTGAAACAGAGAAATTATTAGCCGACa                                                                                                                                      >  1:469448/1‑229 (MQ=255)                                                                 aGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGGGGAACA‑TTTTACGGGTATCTG‑GTTGCGCTGGAAACGGGGGCGGGAAAAGGGAAATTATTAGCCGCCa                                                                                                                                      >  1:306763/1‑229 (MQ=255)                                                                          aggATGTCGAATGGCTGTTTGAATCTCGTAGCGATTTTTTCTGAGGTGATGAGGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTG‑GTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCg                                                                                                                             <  2:469448/227‑1 (MQ=255)                                                                          aaGATGGCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTG‑GTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCg                                                                                                                             <  1:568444/229‑1 (MQ=255)                                                                                  gAAGAGCTGCTTGAAACTCGTATCGATCTTTTCTGAGGAGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGGAATTGCTTTGACTCTGGATGAAGCCAGAACCCAGGGGCGGGGAAGTGTAACC‑TTTTACGGATAATTGAGTGACCCAGAAAAGGGA‑ACTGAAACAGAGAAATAAGATGGCGACATTaaatcacgacgtaaa                                                                                                                     >  2:262897/1‑216 (MQ=255)                                                                                    agagGTGATTGAAACTCGTGGCGATCTGTTCTGGGGGGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTG‑GTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGc                                                                                                                   <  2:270388/229‑1 (MQ=255)                                                                                    agagCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTG‑GTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGc                                                                                                                   <  1:89352/229‑1 (MQ=255)                                                                                          gctTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGACGAAAACATTACTTCTTTGTGCGTTTCTTGTTGTGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGCACTCAGGGGCGGGTGAGTGCAAAA‑TTTTACGGTTATCTT‑GTGTCGCTGAAAAAGGGTCGTTGAACTGAGAAACTCATCGCCCACTTTAAAGCCTAACATaaaacagagga                                                                                                             >  2:342836/1‑221 (MQ=255)                                                                                             tGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTTCGGTGATCTG‑GATGCGCTGGAAAAGGTGGCTGAAACAGAGCAATTATTGACCGACATTAATACCGAAAGAAAAGAGAGTTAcaa                                                                                                          >  2:520783/1‑227 (MQ=255)                                                                                                       tgcggttttttttGGGGTGATGTAGAAGAAAACATTACTTGTTTGTGCGTTTCTTGGTGGGAGGGTAAGGGGAAATGCAAGGGCATTGACTTTGGTTGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTG‑GTTGCGATGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGgc‑a                                                                                               <  2:290143/222‑1 (MQ=255)                                                                                                       tgcggttttttttGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTGGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTG‑GTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGgc‑a                                                                                               <  2:306763/222‑1 (MQ=255)                                                                                                                    gcggTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGCAAGTGAAGAA‑TTTTCCAGTGATCTG‑GTTGCGCTGAAAACGGAGATCGCACCAGAGAAATTAGCAACCGCCATTATTGCCGAAACAAAAGCGCAGTACCAACAACTGGA‑CACGCAAaataaaa                                                                                  >  2:41751/3‑227 (MQ=255)                                                                                                                       gTGATTATGAAGAAAACATTCCTTCGTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTGTACGGTTATCTG‑GTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATgtgt                                                                               <  1:302905/229‑1 (MQ=255)                                                                                                                         gatgatgaAGAAAGTATTACTTCTTTGTGCGGTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTATGGATGAAGCCAGAACTCAGGGGCAGGTAGGTGAAACA‑TTTTACGGTTATCTG‑GTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCg                                                                             <  1:345058/229‑1 (MQ=255)                                                                                                                          atgatgaAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTG‑GTTGCGCTGAAAACg                                                                                                                                                                     <  1:717270/141‑1 (MQ=255)                                                                                                                          atgatgaAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGGGGAACA‑TTTTTCGGGTATCTG‑GTTTCGCGGGAAACGGATGCGTAAACAAATCAATGAGGAGCCGACATTAATGCCCAACGGAAAGAGAGTGAACAAAAACAGGcgaaagcaacaaaagttgtc                                                                              >  2:522893/1‑217 (MQ=255)                                                                                                                          atgatgaAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAAGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTG‑GTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGCGTTACAAACAACTGGC‑AAAGCAAAATACTGTGACg                                                                             >  1:405016/1‑228 (MQ=255)                                                                                                                                  gaaAACTTTTCTTTTTTGTGGGTTTGTTGTTGGTCTGGTAAGGGGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTG‑GTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGTAgatga                                                                     <  2:422979/228‑1 (MQ=255)                                                                                                                                          tACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGATAGGGGAAAAA‑TTTTAAGGTTATTTG‑GATGGGCTGCAAACTGATTCTGAAACCGCGGAATTTGTAGCCGGCATTAATGTCGAACGTACAGCGAGTTACCCACAACAGGG‑CAAGCCAAAACATGAGTCGGTGGCTGatatatcgca                                                            >  1:543917/1‑223 (MQ=255)                                                                                                                                                         ttCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTG‑GGTGCGCTGAAAACGGATGCTGGAACAGAGAAATTTGTAGCCGACATTAATGCCGCACGTAAAGCGGGTTTCCAAAAACAGGG‑ACAGCCAAAAAAAAGGGCGGTAGAGGGGATCCGCAAACACCCCGGGAaaaa                                             >  2:454557/1‑229 (MQ=255)                                                                                                                                                                   gCTGGTAAGGAGCAATGTAATGGCATTGACTTTGGATGAAGCCAGAACTGAGGGGCGGGTAGGTTGAACA‑TTTTACGGTTATCTG‑GTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAGCGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACGCGCCGGGCAAAAGCTCGTAGCt                                   <  1:41751/229‑1 (MQ=255)                                                                                                                                                                            cagcaATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTTCGGTTATCTG‑GTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTATTAGCCGACATTAAAGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGGGGCGGTAAATGGAATCGCGAAACTCGCCGGGCCAAAGATCGGAGCTCGGGCCaaa                          >  1:523691/1‑229 (MQ=255)                                                                                                                                                                             agcaATGTAATGGCATTGACTTTGGATGAAGCCAGAACTCAGGGCCGGGTAGGTGAAACA‑TTTTACGGTTATCTG‑GATGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTCGGGCCaaa                          <  1:520783/228‑1 (MQ=255)                                                                                                                                                                              tgaaTGTAAGGGCATGGACTCTGGATTAAGGCAGAACTCGGGGGCGGGGAGGGGAAACA‑TTTTACGTTTAACTG‑GTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAACCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACCCGCCGGGCAAAAGCTCGTAGCTCGGGCCAAAc                         <  2:744833/226‑1 (MQ=255)                                                                                                                                                                                 gtGTAATGGAGTTGACTTTGGATGAACACAGAACTAAGGGGCGGGTTGGGGAAACA‑TTTTACGGTTATCTG‑GTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGGCGACATTAATGCCGAACGTAAAGCGAGTTACCAAAAACTGGC‑AAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTCGGGCCAAACCggg                     <  1:338942/228‑1 (MQ=255)                                                                                                                                                                                 aTGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGAGGGGAAACA‑TTTTACGGTTATCTG‑GTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGGAGATGATATCGCGAAACTCACCGGGCACAAAATCGTAGCTCGGGCCAAAcacgg                     >  1:423455/1‑225 (MQ=255)                                                                                                                                                                                    tAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGAGGTGAAACA‑TTTTACGGTTATCTG‑GGTGCGCTGAAAACGGATGCTGAAACAAAGAAATTATTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGGAGATGATATCGCGGAACTCGCCGGGCAAAAGACCGTAGCTCGGGCCAAACCGCGGGa                  >  1:119037/1‑229 (MQ=255)                                                                                                                                                                                            tttatttgggATGAAGCCAAACCGCAGTGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTG‑GTTGCGGTGAAAACGGATGCGGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTCGGGCCAAACCGGGGGAATATGTAc          <  2:184659/221‑1 (MQ=255)                                                                                                                                                                                            ttGCCTCTGGATGAAGCAAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTG‑GTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTCGGGCCAAACCGGGGGAATATGTAc          <  2:622613/229‑1 (MQ=255)                                                                                                                                                                                            ttGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGGGAAACA‑TTTTACGGTTATCTG‑GGTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGCACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGGAGATGGTATCGCGGAACACGCCGGGCAAAAGCTCGAAGCTCGGGCCAAACCGAGGGAATATGtca          >  1:184659/1‑227 (MQ=255)                                                                                                                                                                                              gactttggataaatgtagacaccAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTG‑GTTGCGCTGAAAACGGATGCGGAAACAGAGAAATTAGTAGCCGGCATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTCGGGCCAAACCGGGGGAATATGTACaa        <  2:618513/207‑1 (MQ=255)                                                                                                                                                                                              gACTTTGGATGAAGCCTGAACCGAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTAGCTG‑GTTGCGCTGAAAACGGATGCTGAAACGGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTCGGGCCAAACCGGGGGAATATGTACaa        <  2:296512/229‑1 (MQ=255)                                                                                                                                                                                                     gATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTG‑GTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTCGGGCCAAACCGGGGGAATATGTACAAGGGATa  <  1:413740/228‑1 (MQ=255)                                                                                                                                                                                                                            ggTAGGTGAAACA‑TTTTACGGTTATCTG‑GTTGCGCTGAAGACGGATCCTGAAACAGAGAAATTAGTGGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑AAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTCGGGCCAAACCggg                     <  1:405/186‑1 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                          GCTGCACCTAAGGAACCGATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACA‑TTTTACGGTTATCTG‑GTTGCGCTGAAAACGGATGCTGAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGC‑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: ‑