Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,068,537	G→A	7.3%	A91V (GCA→GTA)	argO ←	arginine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,068,537	0	G	A	7.3%	193.4 / 3.1	64	A91V (GCA→GTA)	argO	arginine transporter
Reads supporting (aligned to +/- strand):  ref base G (24/33);  new base A (2/2);  total (28/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

ATAATGCGCTGTGCTTTTGCCGTGCGCAGACGCGGTGCCAGCCAGGCTGCGAGAAGAGCCAGACCAAAGAACCACAGGAAAGAGGCGCTAATTGTCCCGAGTGCAAACCAGCGTTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTT  >  NC_000913/3068261‑3068811                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                   gtaagtcgcTGTGCTTTTGCCGTGCGGAGACGCGGTGCCAGCCAGGCTGCGAGAAGAGCCAGACTACAGAACCACGGGAAAGAGGCGCTAATTGGCCCGAGTGCAAACCAGCGTGTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAATAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTAGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTa                                                                                                                                                                                                                                                                            <  2:5052/279‑1 (MQ=255)  tAATGCGCTGTGCTTTTGCCGTGCGCAGACGCGGTGCCAGCCAGGCTGCGAGAAGAGCCAGACCAAAGAACCACAGGAAAGAGGCGCTAATTGTCCCGAGTGCAAACCAGCGTTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCa                                                                                                                                                                                                                                                             >  1:327565/1‑299 (MQ=255)   aaTGCGCTGTGCTTTTGCCGTGCGCAGACGCGGTGCCAGCCAGGCTGCGAGAAGAGCCAGACCAAAGAACCACAGGAAAGAGGCGCTAATTGTCCCGAGTGCAAACCAGCGTTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTTCTGTTTTaaa                                                                                                                                                                                                                                                                          >  2:461734/1‑285 (MQ=255)   aaTGCGCTGTGCTTTTGCCGTGCGCAGACGCGGTGCCAGCCAGGCTGCGAGAAGAGCCAGACCAAAGAACCACAGGAAAGAGGCGCTAATTGTCCCGAGTGCAAACCAGCGTTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTaaa                                                                                                                                                                                                                                                                          >  2:242356/1‑285 (MQ=255)   aaTGCGCTGTGCTTTTGCCGTGCGCAGACGCGGTGCCAGCCAGGCTGCGAGAAGAGCCAGACCAAAGAACCACAGGAAAGAGGCGCTAATTGTCCCGAGTGCAAACCAGCGCTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAGCAAAAGTATCCAGGTAAACATGCGGCTTCAGCCAGGTCACTGCCAACACGATGGCGATAATTTTCCATCTGCTTTCCTTCATGACTTCGCCGCTGGCTCACTCACTATTACTGCTCATTCCTTTTTTaa                                                                                                                                                                                                                                                                           >  2:126370/1‑284 (MQ=255)    aTGCGCTGTGCTTTTGCCGTGCGCAGACGCGGTGCCAGCCAGGCTGCGAGAAGAGCCAGACCAAAGAACCACAGGAAAGAGGCGCTAATTGTCCCGAGTGCAAACCAGCGTTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTaaaa                                                                                                                                                                                                                                                                         <  2:57763/285‑1 (MQ=255)                ttttCCGTGCGCAGACGCGGGGCCAGCCAGGCTGCGAGAAGAGCCAGACCAAAGAACCACAGGAAAGAGGCGCTAATTGTCCCGAGTGCAAACCAGCGTTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAAcc                                                                                                                                                                                                                                                              <  2:304074/284‑1 (MQ=255)                      gTGCGCAGACGCGGTGCCAGCCAGGCTGCGAGAAGAGCCAGACCAAAGAACCACAGGAAAGAGGCGCTAATTGTCCCGAGTGCAAACCAGCGTTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTAc                                                                                                                                                                                                                                         >  1:430026/1‑299 (MQ=255)                       atcgcAGACGCGGTGCCAGCCAGGCTGCGAGAAGAGCCAGACCAAAGAACCACAGGAAAGAGGCGCTAATTGTCCCGAGTGCAAACCAGCGTTTTGATTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACcaca                                                                                                                                                                                                                                                      <  2:288617/283‑1 (MQ=255)                         cgcAGACGCGGTGCCAGCCAGGCTGCGAGAAGAGCCAGACCAAAGAACCACAGGAAAGAGGCGCTAATTGTCCCGAGTGCAAACCAGCGTTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTAcgc                                                                                                                                                                                                                                       <  1:564228/298‑1 (MQ=255)                            aGACGCGGTGCCAGTCAGGCTGCGAGAAGAGCCAGACCAAAGAACCACAGGAAAGAGGCGCTAATTGTCCCGAGTGCAAACCAGCGTTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTACTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTAcgccgc                                                                                                                                                                                                                                    >  1:330064/1‑298 (MQ=255)                              aCGCGGGGCCAGCCAGGCTGCGAGAAGAGCCAGACCAAAGAACCACAGGAAAGAGGCGCTAATTGTCCCGAGTGCAAACCAGCGTTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCaagaa                                                                                                                                                                                                                                               <  2:302141/285‑1 (MQ=255)                                       cagccagGCTGCGAGAAGAGCCAGACCAAAGAACCACAGGAAAGAGGCGCTAATTGTCCCGAGTGCAAACCAGCGTTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCa                                                                                                                                                                                                                        <  1:274749/299‑1 (MQ=255)                                               cTGCGAGAAGAGCCAGACCAAAGAACCACAGGAAAGAGGCGCTAATTGTCCCGAGTGCAAACCAGCGTTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGc                                                                                                                                                                                                                >  1:123676/1‑299 (MQ=255)                                                           ccAGACCAAAGAACCACAGGAAAGAGGCGCTAATTGTCCCGAGTGCAAACCAGCGTTTTGGTTCCACATCAAGTTGCACGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACCTGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACAACAGCAAGAAGGCTACTCCGCCCCAGGTGACCAgcgc                                                                                                                                                                                                                    >  2:421785/1‑283 (MQ=255)                                                                    aGAACCACAGGAAAGAGGCGCTAATTGTCCCGAGTGCAAACCAGCGTTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTACCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAAt                                                                                                                                                                                           <  1:129095/299‑1 (MQ=255)                                                                         cacaGGAAAGAGGCGCTAATTGTCCCGAGTGCAAACCAGCGTTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGAc                                                                                                                                                                                                    >  2:338272/1‑285 (MQ=255)                                                                                 atagtcgcTAATTGTCCCGAGTGCAAACCAGCGTTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCaa                                                                                                                                                                                            <  2:432707/280‑1 (MQ=255)                                                                                     gcgcTAATTGTCCCGAGTGCAAACCAGCGTTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAAc                                                                                                                                                                                        >  2:157675/1‑285 (MQ=255)                                                                                           aTTGTCCCGAGTGCAAACCAGCGTTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATccc                                                                                                                                                                    >  1:169966/1‑299 (MQ=255)                                                                                             tGTCCCGAGTGCAAACCAGCGTTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCt                                                                                                                                                                                   <  2:310140/282‑1 (MQ=255)                                                                                                                    ttGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCCCACATCa                                                                                                                                                         >  2:196298/1‑285 (MQ=255)                                                                                                                     tGGTTCCACATCAAGTTGCCCGCCAAGGCAGCCCAGTACAACAAAAGGATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATc                                                                                                                                                          <  2:314145/283‑1 (MQ=255)                                                                                                                       gTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAAAAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGa                                                                                                                                                      <  2:287456/285‑1 (MQ=255)                                                                                                                       gTTCCACATCAAGTTGACCGCCAAGGCTGCCCAGTACAACAAAAGTATGCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTGCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAgg                                                                                                                                                       <  2:108462/284‑1 (MQ=255)                                                                                                                          ccACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGca                                                                                                                                     <  1:568587/299‑1 (MQ=255)                                                                                                                            acaTCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCACCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATACCGGCGCAAATCAGGACCaac                                                                                                                                                 >  2:474390/1‑284 (MQ=255)                                                                                                                             caTCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCCAAAATCCCGGCGCAAATCAGGACCAAAt                                                                                                                                                >  2:22344/1‑285 (MQ=255)                                                                                                                               tCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATc                                                                                                                                               <  2:77510/284‑1 (MQ=255)                                                                                                                                    ttGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAg                                                                                                                            >  1:562804/1‑298 (MQ=255)                                                                                                                                                   gCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACCGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAgg                                                                                                                           <  2:22325/284‑1 (MQ=255)                                                                                                                                                    cccAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGt                                                                                                            >  1:427515/1‑298 (MQ=255)                                                                                                                                                                     aTCCATGTAAACATGCGGATTCAGCCAGGTCACTGCAAACATGGTGGGGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGAACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGcc                                                                                            <  1:555587/297‑1 (MQ=255)                                                                                                                                                                     aTCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGcc                                                                                            <  1:167555/297‑1 (MQ=255)                                                                                                                                                                        cAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTACTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTga                                                                                                       <  2:289758/283‑1 (MQ=255)                                                                                                                                                                         aGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTgacg                                                                                                     >  2:203377/1‑284 (MQ=255)                                                                                                                                                                         aGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTAATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGAGCAGCGCCAGCAACCACGGCGGCTGCATCAATAACGCGCTGCCACCAACAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGGAAAGCAATCATAATGTGGTACTgacg                                                                                                     >  2:312370/1‑284 (MQ=255)                                                                                                                                                                             aaaCATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATtcatc                                                                                  <  1:223251/299‑1 (MQ=255)                                                                                                                                                                              aaCATGCGGAGTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTGCTGACGACGTAt                                                                                               <  2:376422/285‑1 (MQ=255)                                                                                                                                                                                 aTGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACaa                                                                              >  1:483638/1‑299 (MQ=255)                                                                                                                                                                                 aTGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCACCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACaa                                                                              >  1:327868/1‑299 (MQ=255)                                                                                                                                                                                   gCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGAATCAATAACGCGCTGCCACCACAAATCCCGGCGCAACTCAGGCCCAAATCGCTGATAGCACAAAGCAAGTCAATCATTATGTGGTAATTACGACGTAAGCCCTGATTCATTACaaa                                                                             >  1:335398/1‑298 (MQ=255)                                                                                                                                                                                   gCGGATTAAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGATTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCAGTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCAAGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCt                                                                                          <  2:514980/285‑1 (MQ=255)                                                                                                                                                                                          cAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCGGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACC‑AAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATtcat                                                                                   <  2:521387/284‑1 (MQ=255)                                                                                                                                                                                            gacAGGTCACTGCCCCCATCCTTGCTCTAATTTTCAATCTCCCTTGCTTCATCTCTTCGGCGCAGGCTACCTCACTATTACTGCTCATTGCTGTTTTAAAAGCGCCAACACCATCCCCCAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATAAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATtcatca                                                                                 <  2:452494/283‑1 (MQ=255)                                                                                                                                                                                                 gTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAACACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGCCCAAATCGCTGATAGCACAAAGTCAGGCAATCATAATGTGGTACTGCCTACGTATGCCCTGATTCATCACaaaa                                                                            >  2:14335/1‑285 (MQ=255)                                                                                                                                                                                                  tCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAg                                                                           >  2:45150/1‑285 (MQ=255)                                                                                                                                                                                                    aCTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAg                                                           <  1:231522/299‑1 (MQ=255)                                                                                                                                                                                                         cAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGt                                                       <  1:52597/298‑1 (MQ=255)                                                                                                                                                                                                         cAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTa                                                      >  1:264738/1‑299 (MQ=255)                                                                                                                                                                                                         cAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTa                                                      <  1:144003/299‑1 (MQ=255)                                                                                                                                                                                                          aaCATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTACTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAg                                                     >  1:504796/1‑299 (MQ=255)                                                                                                                                                                                                            cATGGTGGCTATCATTTTCCATCTGCCTTGCTTCATGACTTCGTCGCTGGCTAACTCAATATTACTGCTCATTTCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCTACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGCCCGAGCGGTAGGa                                                   <  1:459929/299‑1 (MQ=255)                                                                                                                                                                                                                gtggCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAATGCTTTtgtgat                                                                <  2:149286/282‑3 (MQ=255)                                                                                                                                                                                                                                 tCTGCCTTGCTTCATGACTTCGGCGCTGGATAACTCAATATTACTGCACATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCGCGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCaa                              <  1:483001/299‑1 (MQ=255)                                                                                                                                                                                                                                     ccTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAgccgc                                         <  2:518279/284‑1 (MQ=255)                                                                                                                                                                                                                                            tCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAg                   <  1:529320/299‑1 (MQ=255)                                                                                                                                                                                                                                                  cTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGAACATAGCCGCCCCAAGGGCAAGAc                           >  2:130366/1‑285 (MQ=255)                                                                                                                                                                                                                                                         gcTGGCTAAATCAAAATTACTGCGCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGaa                     <  2:233103/284‑1 (MQ=255)                                                                                                                                                                                                                                                          cTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAg                   <  2:374194/285‑1 (MQ=255)                                                                                                                                                                                                                                                             gCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGtt  >  1:477942/1‑299 (MQ=255)                                                                                                                                                                                                                                                                aaCTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGtaata              <  2:133788/284‑1 (MQ=255)                                                                                                                                                                                                                                                                aaCTCAATATTACTGCTCATTACTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGtaataa             <  2:440834/285‑1 (MQ=255)                                                                                                                                                                                                                                                                  cTCAATATTACTGCTCATTGCTGTTTTCAACGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACACAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTCCAAGACCTTGAAAGTAATAAGa           >  2:278948/1‑285 (MQ=255)                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                   ATAATGCGCTGTGCTTTTGCCGTGCGCAGACGCGGTGCCAGCCAGGCTGCGAGAAGAGCCAGACCAAAGAACCACAGGAAAGAGGCGCTAATTGTCCCGAGTGCAAACCAGCGTTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAACGCGCTGCCACCAAAAATCCCGGCGCAAATCAGGACCAAATCGCTGATAGCACAAAGTAAGGCAATCATAATGTGGTACTGACGACGTATGCCCTGATTCATCACAAAAGCATTTTGTGGACCGAGCGGTAGGATCATAGCCGCCCCAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTT  >  NC_000913/3068261‑3068811

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