Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,069,065	G→A	19.4%	L259L (CTG→TTG)	rutG ←	pyrimidine permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,069,065	0	G	A	19.4%	135.2 / 4.9	66	L259L (CTG→TTG)	rutG	pyrimidine permease
Reads supporting (aligned to +/- strand):  ref base G (38/13);  new base A (5/8);  total (43/23)
Fisher's exact test for biased strand distribution p-value = 2.08e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTAAACACCGCCACCAGCCCAATACACAG  >  NC_000913/1068838‑1069299                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            cAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGcc                                                                                                                                                                                                                                  >  1:727651/1‑238 (MQ=255)      tGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGG‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTCAGTTTTCCGCCCC‑CCAAA‑TCACCCCCCCggg                                                                                                                                                                                                                             >  2:137815/1‑238 (MQ=255)      tGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGG‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACggg                                                                                                                                                                                                                             >  1:636815/1‑238 (MQ=255)                   agcagcaACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCCACGAATGCCCGCCCCCTATACGGGT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAAATGCCCTAGGTTTTCCCCCAC‑CAAAA‑TCCCCGCCACGGGCGCacaaaacacaa                                                                                                                                                                                                               >  2:60373/1‑229 (MQ=255)                          acaaacaCCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTg                                                                                                                                                                                                        >  1:415125/1‑239 (MQ=255)                                acCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCAGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATAACCTAAGCTTCACGCCAC‑CCAAA‑TCACCCCCACGGGCCCAATCACCACAACTGCCCGACCact                                                                                                                                                                                                  >  2:166159/1‑237 (MQ=255)                                  cAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCacaaa                                                                                                                                                                                                >  2:717393/1‑235 (MQ=255)                                  cAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCACCATTGCCTGACCATTaa                                                                                                                                                                                                >  1:758741/1‑239 (MQ=255)                                  cAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAAAA‑TCACCGCCCCGGGCGCAATCAGCACCATTGCCTGACCATaaa                                                                                                                                                                                                >  2:20878/1‑239 (MQ=255)                                     cGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCACC‑CCACA‑TCACAGCCCCGGTCGCAATCAGCATCATTGCCGGACCATTaaaaa                                                                                                                                                                                             >  1:1015465/1‑238 (MQ=255)                                          aGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTCAGTTTTCCGCCCC‑CAGAA‑TCACCGCCCCGGGGGCAAGCAGACTAATTGCCGGACCATTaaaaaaaggg                                                                                                                                                                                        >  2:382254/1‑233 (MQ=255)                                          aGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGC‑CCCTATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCCCCAC‑CAGAA‑TCACCGCCACGGGCGCAACCAGCATCATTGCCTGACCCTTAAAAACAggg                                                                                                                                                                                        >  2:1033884/1‑238 (MQ=255)                                                                     cGATATTTTCAGCATAGGTGGTCACGCCGCTACTGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCACCATCATTGCCTGCCCCTTAAAGACAGGAGGGGAAAAAGTGGGCAAACCGAACCAg                                                                                                                                                             >  1:466130/1‑239 (MQ=255)                                                                     cGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCGAACCGAACCAg                                                                                                                                                             >  1:420925/1‑239 (MQ=255)                                                                            ttCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCAt                                                                                                                                                      >  1:47546/1‑239 (MQ=255)                                                                                            aCGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAGACAGAAGTGGAAAAATGGGGCAGACCGAACCAGGCAGCATGGCGGACGAGAGTaaa                                                                                                                                      >  1:158322/1‑239 (MQ=255)                                                                                                             aaggccgaaagccttCGTGGCCAGCCCTCCGCCGACGAAGCCCCGCCCCATAAAGGGGT‑CCATTTTG‑CGTCGGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTGCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCa                                                                                                                     <  1:20878/225‑1 (MQ=255)                                                                                                               agCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCACCGGCTTTTCCCat                                                                                                                    >  1:712751/1‑237 (MQ=255)                                                                                                                               tGCCAGACCA‑CCCCTAAGGATGCCCCCCCCAAAAACCGGG‑CCAATTTGCCGCCCGGTCACCCCGGCTACTGCCTTGAGATGACCTAAGTTTCCCGCCCC‑CAAAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGtca                                                                                                   <  1:382254/239‑1 (MQ=255)                                                                                                                                  cagcccaccccttccaaaTGCCCGCCCCATACACGGGT‑CCATATTG‑CGCCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTCCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGtcatca                                                                                                <  1:137815/223‑1 (MQ=255)                                                                                                                                   acaccctcccctccccgatccccccccATTTTCGGGTCCCATTTGG‑CGCCCGGTCCACCCGGCGTCTGCCTTGAGATGACCTAAGGTTTCCCCCCA‑CAAAA‑GCGCCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGtcatca                                                                                                <  2:712751/220‑1 (MQ=255)                                                                                                                                    gccaccccccatcgcacgcccgcccGATATCCGGTC‑CCATTTGG‑GCTCCGGCCATCCCGGCGCCTGCCTTGGACGGCCCAAAGTTTTCCCCCCC‑AACAA‑TCCCCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACg                                                                                              <  2:658788/222‑1 (MQ=255)                                                                                                                                         tCGCCTACGAATGCCCGCCCCATATACGGGT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACTTTGGTCATCACGCCATa                                                                                         >  1:207454/1‑239 (MQ=255)                                                                                                                                                       ccgcccCATATACGGGT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGCCATCACGCAATACAGAAAACAgcgca                                                                           >  2:318339/1‑235 (MQ=255)                                                                                                                                                       ccgcccCATATACGGGT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCACCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGGCCTAACCCCAAAAAGCCGAAAGGGGa                                                                           >  2:425191/1‑239 (MQ=255)                                                                                                                                                                     ggT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGGGTAAAATCGACGGCTTTTCCCAACCCCAAAAAGTTGGGCATAACCCCATAAAGCAAACAGGCCACAACAAATCCAAcc                                                             >  2:868733/1‑238 (MQ=255)                                                                                                                                                                     ggT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGCCAATTAATCCTACg                                                             >  1:143098/1‑239 (MQ=255)                                                                                                                                                                     ggT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAAAAAAGCGAAACTTAAACCTAca                                                             >  2:222497/1‑238 (MQ=255)                                                                                                                                                                      gT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCCCGCCATACCGCAGACAgggcaaaaataaacccacgg                                                            >  2:1025648/1‑221 (MQ=255)                                                                                                                                                                       t‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACgag                                                           >  1:201944/1‑239 (MQ=255)                                                                                                                                                                          cATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGCAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGGAAAATCGACGCCTTTTCCCAGCCCTAAAACGTTGGGCATCACCCCAAAAAGCAAAAAGGCGAAAAATAAACCCACggaaa                                                         >  2:504762/1‑235 (MQ=255)                                                                                                                                                                          cATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCCGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATAAATCCAAAGAGaa                                                         >  2:762367/1‑239 (MQ=255)                                                                                                                                                                          cATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACgaga                                                          >  1:330233/1‑238 (MQ=255)                                                                                                                                                                          cATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCAAGCCATACAGCAGACAGGCGACACTTAATCCTACGCGaa                                                         >  1:764674/1‑239 (MQ=255)                                                                                                                                                                          cATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCAAGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACCTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGaa                                                         >  1:576452/1‑239 (MQ=255)                                                                                                                                                                          cATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGCAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGAAAAATCGCCGGCCTTTCCCAGCCCTAAAACGCTGGTCATCACGCCATAAAGCAGACAGGCGACGCTTAATCCTACGAGaa                                                         >  1:1051990/1‑239 (MQ=255)                                                                                                                                                                                ttctgcgggtctcccgggcgctggCCTGGAGATGCCTGATGTTTCCCGCACC‑AAGAA‑TCACCGGCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAAGCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAAt                                                    <  2:633324/216‑1 (MQ=255)                                                                                                                                                                                tgcccccGGTCCCCCCGGGGACTTCCTTGAGTTGCCCTATGTTTGCCGCCCC‑CAAAA‑CCCCCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAAt                                                    <  2:143098/234‑1 (MQ=255)                                                                                                                                                                                gtcccccGGTCACCCCGGCGCCTCCCTGGAGATGACTTGATTTTCCCCCCAC‑AATGA‑TCCCCGACACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTCTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAAt                                                    <  2:415125/234‑1 (MQ=255)                                                                                                                                                                                gccgtccctcccccccGGCGCCTCCCTTGAGATGACC‑AAGTTTTCCCCCCCACAGAA‑CCCCCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAAt                                                    <  2:47023/226‑1 (MQ=255)                                                                                                                                                                                gccGTCGGGTCACCCGGGCGATTGCCTTGAGAGGACCTAAGTTTTCCGCCAC‑CAAAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAAt                                                    <  1:955186/237‑1 (MQ=255)                                                                                                                                                                                    tCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGCCAGGCGACAATTAATCCTACGAGAATCAATAAc                                                 >  1:273999/1‑238 (MQ=255)                                                                                                                                                                                    tCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAAcc                                                >  1:47023/1‑239 (MQ=255)                                                                                                                                                                                        gTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAAGCAGCCTGGCTGACGCGAGGAAAATCGACGGCTTTTCCCGGCCCTAAAACGTTGGTCATCACGCCAAACAGCGAAAAGGCGACAATTCATCCTACGAGAAACAATAAcccccg                                            >  2:590637/1‑235 (MQ=255)                                                                                                                                                                                                  cctcttccttgggtttccttatttttcCCCCCAC‑CAGAA‑CCCCCGCCCCGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTCGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCa                                  <  2:420925/216‑1 (MQ=255)                                                                                                                                                                                                   cacctgtggggagccaCCTAATTTTCCCGCCCA‑CAGAA‑CCCCCGCCAAGGGCCCAAGCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTTGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCa                                  <  2:676366/223‑1 (MQ=255)                                                                                                                                                                                                     cTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGCCAGGGGAAAATTAATCCTACGAGAATCAATAACCCCTGAATCATGCCCCgg                               >  1:188282/1‑239 (MQ=255)                                                                                                                                                                                                      tGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCACTAACCGCGGGATCATCCCACggg                              >  1:724085/1‑239 (MQ=255)                                                                                                                                                                                                       gCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGt                             >  1:656639/1‑239 (MQ=255)                                                                                                                                                                                                        ccTGGAGATCACTTAATTTTCCCGCCAC‑CAGAC‑TCCCCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTa                            <  1:717393/239‑1 (MQ=255)                                                                                                                                                                                                        catggagaTGACATAATTTTTCCGCCAC‑TAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGa                                                                                                                                 <  2:268438/134‑1 (MQ=255)                                                                                                                                                                                                                cttccaacttttcCCCCCCC‑CAAAA‑CCCCCGCCCCGGGCGAAATCAGCACCATTGCCTGACCATTAAAAGCAGGAGGGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTAAACAcc                      <  2:194016/228‑1 (MQ=255)                                                                                                                                                                                                                  aCCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGACCCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTcc                                                                                                                       >  1:268438/1‑138 (MQ=255)                                                                                                                                                                                                                  aCCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGCCAGGCGACCATTAATCCTACAAGAATCAAAAACCGCTCGATCATGCCACAGGAAAAAAACGcc                   >  1:501964/1‑238 (MQ=255)                                                                                                                                                                                                                   cttAATTTTCCCGCCCC‑AAAAA‑CCACCGCCACGGGCCCAATCAGCATCATTTCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTAAACACCGccac                 <  2:743167/237‑1 (MQ=255)                                                                                                                                                                                                                   ccTAAGTTTCCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTAAACACCGccac                 <  1:60373/239‑1 (MQ=255)                                                                                                                                                                                                                       aGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGCCAGGCGCCAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTAAACACCGCCACCAg              >  1:743167/1‑238 (MQ=255)                                                                                                                                                                                                                          ttcccGCCAC‑AAGAA‑CCCCCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTAAACACCGCCACCAGCCCa          <  1:222497/236‑1 (MQ=255)                                                                                                                                                                                                                          tatCCGCCCA‑CAGAA‑TAACCGGCACGGGCGCAATCAGCAGCACTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGGGGCAGACCGAACCAGGCAGCATGGCTGACTAGAGTAAAATCGACGGCTTTTCCCATCCCGAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTAAACCCCGCCACCAGCCCa          <  1:590637/237‑1 (MQ=255)                                                                                                                                                                                                                           ttCCCCCAC‑CAAAAGTCCCCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAAGGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTAAACACCGCCACCAGCCCa          <  1:166159/239‑1 (MQ=255)                                                                                                                                                                                                                            tCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGgcagca                                                                                                                                                       <  2:1037609/96‑1 (MQ=255)                                                                                                                                                                                                                            tCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGgcagca                                                                                                                                                       >  1:1037609/1‑96 (MQ=255)                                                                                                                                                                                                                              cGCCAC‑CAAAA‑CCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTAAACACCGCCACCAGCCCAATac      <  1:425191/239‑1 (MQ=255)                                                                                                                                                                                                                                 cac‑caGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTa                                                                                                               >  1:788697/1‑131 (MQ=255)                                                                                                                                                                                                                                 cac‑caGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAA‑TAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTAAACACCGCCACCAGCCCAATACACAg  >  1:223544/1‑239 (MQ=255)                                                                                                                                                                                                                                 cac‑aaGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTa                                                                                                               <  2:788697/131‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            CAGCATGGCGATGACGGCAGCAGCAACAAACACCAGCGTCGAGTAGACTTTGGTCACGGCCATCACACCGATATTTTCAGCATAGGTGGTCACGCCGCTACCGCCGACAGAGCCGGAAAGCATCGTTGCCAGACCATCGCCTACGAATGCCCGCCCCATATACGGGT‑CCATATTG‑CGTCCGGTCATCCCGGCGACTGCCTTGAGATGACCTAAGTTTTCCGCCAC‑CAGAA‑TCACCGCCACGGGCGCAATCAGCATCATTGCCTGACCATTAAAAGCAGGAGTGGAAAAATGTGGCAGACCGAACCAGGCAGCATGGCTGACGAGAGTAAAATCGACGGCTTTTCCCAGCCCTAAAACGTTGGTCATCACGCCATACAGCAGACAGGCGACAATTAATCCTACGAGAATCAATAACCGCTGGATCATGCCACGGGTAAACACCGCCACCAGCCCAATACACAG  >  NC_000913/1068838‑1069299

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