Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,581,316	G→A	G60E (GGA→GAA)	yhhY →	N‑acetyltransferase YhhY

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,581,316	0	G	A	98.4%	209.2 / NA	61	G60E (GGA→GAA)	yhhY	N‑acetyltransferase YhhY
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (27/33);  minor base C (1/0);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCATTTGACTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGACT  >  NC_000913/3581072‑3581560                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     gCATTTGACTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACATGTGCCTCAGCCTTCCGATCATATGTGGCAGTATCGACTCGCCGATCGGCCCGGCATCAAGCAACTCGTCGCCTGTATTGAGGAAGAcg                                                                                                                                                                                                                                                 <  2:305728/250‑1 (MQ=255)  cATTTGACTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGAcgtc                                                                                                                                                                                                                                               >  2:113575/1‑251 (MQ=255)     ttGACTCGCATTTGCAAGTCAATGCTGTTGAAAGGGCCATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCACGCAACTCGTCGCCTGTATTGATGAAGACGTCGTg                                                                                                                                                                                                                                            >  2:69415/1‑251 (MQ=255)     ttGACTCGCATTTGAAATTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGAcgtcgt                                                                                                                                                                                                                                             <  1:223742/250‑1 (MQ=255)      tGACTCGCATTTGAAAGTCAATTATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGAGTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTTATGAAGACGTCGTg                                                                                                                                                                                                                                            <  2:297843/250‑1 (MQ=255)       gACTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGACATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTggg                                                                                                                                                                                                                                          >  1:286003/1‑251 (MQ=255)           cGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCa                                                                                                                                                                                                                                       <  2:33968/250‑1 (MQ=255)            gCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGACGACGTCGTGGGCCAtc                                                                                                                                                                                                                                     >  2:8724/1‑251 (MQ=255)              aTTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCAtctc                                                                                                                                                                                                                                   <  1:18791/251‑1 (MQ=255)               tttGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCAtctcc                                                                                                                                                                                                                                  >  2:129786/1‑250 (MQ=255)                  gAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCa                                                                                                                                                                                                                               <  1:190170/251‑1 (MQ=255)                   taaGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACCAAGGAATGAGAGAGATAGTAATACGCCACGCAGCAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTAATGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGCTCGTCCCGGCATCAAGCAACTCCCACCCTGTATTGATGAAGCCGTCGTGGGCCATCTCACCAt                                                                                                                                                                                                                              >  2:294051/2‑251 (MQ=255)                      gTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGCGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCCAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTg                                                                                                                                                                                                                            >  2:222946/1‑250 (MQ=255)                        cAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACg                                                                                                                                                                                                                         <  1:8724/251‑1 (MQ=255)                          atgatgTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGTGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCTGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAACCAACGCGTCGCCTGTAGTGATGAAGACTGCGTGGGCCGTTTCTCCATTGACGTg                                                                                                                                                                                                                       <  2:433121/251‑1 (MQ=255)                            gatgTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGCAGACGTCGTGGGCCATCTCACCATTGACGTGcc                                                                                                                                                                                                                     >  2:240415/1‑250 (MQ=255)                               gTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAAc                                                                                                                                                                                                                   <  1:441022/250‑1 (MQ=255)                                            tttACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGCCGTGCAACAGCGCCCAcgccgc                                                                                                                                                                                                     >  2:89546/1‑251 (MQ=255)                                              tACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCAcgccgc                                                                                                                                                                                                     >  2:3238/1‑249 (MQ=255)                                               cccccAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGt                                                                                                                                                                                                  <  1:3238/250‑1 (MQ=255)                                                 cccAAAGAGGAAAAAGGAATGAGGGAGATAGTAATACGCAACGCAGAAACACGGGATTACGAGGCCCTAAGGCAGAGTCACGCCCCGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATTTGGCAGGAGCGACTCGCCTATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTTATGAAGACGTCGGGATCCATCTCACCATGGACGTGCAACAGCGCCCACGCCTCAGTCa                                                                                                                                                                                                <  2:202462/251‑1 (MQ=255)                                                  ccAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCAt                                                                                                                                                                                               <  2:3145/251‑1 (MQ=255)                                                    aaaGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGt                                                                                                                                                                                             >  2:81032/1‑251 (MQ=255)                                                         ggACAAATGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGGCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCg                                                                                                                                                                                        <  2:274732/251‑1 (MQ=255)                                                          gACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCg                                                                                                                                                                                        >  2:292605/1‑250 (MQ=255)                                                           aCAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCCACTCGTCGCCTGTATTGATGAAGACGTCGTGTGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGAt                                                                                                                                                                                      >  2:246015/1‑251 (MQ=255)                                                             aaaGGAAGGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGAGGGGCAACAGCGCCCACGCCGCAGTCATGTTGCCGCttt                                                                                                                                                                                    >  1:185885/1‑251 (MQ=255)                                                                      agtgagATAGTACTACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGAGCCTCATCCTTCCGATCATATGTGGCAGGAGCGCCTCACCGATAGTCCCGGCATCAAGCAACTCGTCGCATGTATTGCTGAAGACGTCGGTGGCAATCTCACCATTGACGTGCCACAGCTCCCACGCCGCAGTCATGTTGCCGATTTTGATATCtg                                                                                                                                                                           >  2:314227/1‑251 (MQ=255)                                                                           gaTAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCg                                                                                                                                                                      >  1:361442/1‑251 (MQ=255)                                                                            aTAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGa                                                                                                                                                                     <  2:318775/251‑1 (MQ=255)                                                                              aGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGAct                                                                                                                                                                   <  1:222946/251‑1 (MQ=255)                                                                                   tACGCCACGCAGAAACACGGGATTACGCGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCACCTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCgcc                                                                                                                                                                    >  2:194630/1‑243 (MQ=255)                                                                                   tACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGAc                                                                                                                                                                    <  1:194630/245‑1 (MQ=255)                                                                                      gCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTgg                                                                                                                                                           >  1:314467/1‑251 (MQ=255)                                                                                               aaaCACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACcgc                                                                                                                                                  <  1:37129/251‑1 (MQ=255)                                                                                                                  cATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACGGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGAATCGCCGATAGTCCCGGCATCAAGCAACTCGTAGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCAAGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCGCGTTGGAAGGAACGCGGCGTCGCCAGCGACCTg                                                                                                                                >  1:241869/1‑250 (MQ=255)                                                                                                                           aaTTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCgagaga                                                                                                                      <  2:295804/250‑1 (MQ=255)                                                                                                                              tCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCTCCAGCGCCCTGATGCGAGAgatg                                                                                                                    <  1:129786/250‑1 (MQ=255)                                                                                                                                   cccAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCACCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATATGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGa                                                                                                               >  1:315496/1‑250 (MQ=255)                                                                                                                                    ccAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCCCCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGaaa                                                                                                             <  1:294051/251‑1 (MQ=255)                                                                                                                                       gCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAAtgt                                                                                                          <  1:335624/251‑1 (MQ=255)                                                                                                                                          ggAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACCCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCGCTGATGCGAGAGATTATTGAAATGTGCg                                                                                                       <  2:440207/251‑1 (MQ=255)                                                                                                                                            aGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGCGGGCTGATTGAAATGAGCGAc                                                                                                     >  2:75659/1‑251 (MQ=255)                                                                                                                                            aGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGAc                                                                                                     <  1:113575/251‑1 (MQ=255)                                                                                                                                                gtATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACt                                                                                                 <  1:8926/251‑1 (MQ=255)                                                                                                                                                     gCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGtt                                                                                             <  1:320828/250‑1 (MQ=255)                                                                                                                                                               aCAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGt                                                                                       >  1:129027/1‑246 (MQ=255)                                                                                                                                                               aCAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCAGCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATc                                                                                  <  1:226648/251‑1 (MQ=255)                                                                                                                                                                   gTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGAACGTCCCGGCATCAAGCAACTCGTCGACTGTATTGATGAAGTCGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGGCAGCGCCCTGATGCGAGAGATGATTGACATGTGCGACCACTGGTTGCGGGTAGCTCGCAt                                                                              >  1:273150/1‑251 (MQ=255)                                                                                                                                                                                 cGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCgtgt                                                                <  1:389410/251‑1 (MQ=255)                                                                                                                                                                                  gATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGtt                                                               <  1:318500/251‑1 (MQ=255)                                                                                                                                                                                  gATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGtt                                                               >  1:444571/1‑251 (MQ=255)                                                                                                                                                                                       tatGTGGCAGGAGCGACTCGCCGATCGGCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGTCGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAAcc                                                                    <  2:346960/241‑1 (MQ=255)                                                                                                                                                                                               aGGAGCGACTCGCCGATCGGCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGGAGATCGCATTGAACTAACCGTGTGTGTCGATAACGcc                                                  <  2:452795/251‑1 (MQ=255)                                                                                                                                                                                                         cGCCGATCGTCCCGGACTCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCAGTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAAt                                          <  2:444571/249‑1 (MQ=255)                                                                                                                                                                                                                 gTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCAGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCtat                                >  1:417986/1‑251 (MQ=255)                                                                                                                                                                                                                   cccGGCATCAAGCAACTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCTATaa                              <  1:246015/251‑1 (MQ=255)                                                                                                                                                                                                                                  cTCGTCGCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCTATAAAAAATACGGCTTTg                >  1:3327/1‑250 (MQ=255)                                                                                                                                                                                                                                    cgtcggccgcaTTGATGAAGACTTCGTGGGCAATCTCACCATTGACGTGCGACTGCGCCCACGCCGCAGTGATGTTGGCGATTTTGGTATCTGTGACAACTCTCGTTGGACGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGGTTGAAATGTGCGACAACTGGATGCGGGTAGACCGCATTGAACTAACCGTGTTTGTCGATAACGACCCGGCAATTAAGGTCTATAAAAAATACGGCtttaca              >  2:250520/11‑247 (MQ=255)                                                                                                                                                                                                                                       cgCCTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAACGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCTATAAAAAATACGGCTTTGAAAtt           <  1:249158/250‑1 (MQ=255)                                                                                                                                                                                                                                          cTGTATTGATGAAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGACTATAAAAAATACGGATTTGAAATTGa         >  2:29650/1‑249 (MQ=255)                                                                                                                                                                                                                                                tgatgaAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGGCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGACt  <  1:261376/250‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     GCATTTGACTCGCATTTGAAAGTCAATGATGTTGAAAGGGACATTTACCCCAAAGAGGACAAAGGAATGAGTGAGATAGTAATACGCCACGCAGAAACACGGGATTACGAGGCCATCAGGCAGATTCACGCCCAGCCGGAGGTGTATTGCAACACACTACAGGTGCCTCATCCTTCCGATCATATGTGGCAGGAGCGACTCGCCGATCGTCCCGGCATCAAGCAACTCGTCGCCTGTATTGATGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGTGTTTGTCGATAACGCCCCGGCAATTAAGGTCTATAAAAAATACGGCTTTGAAATTGAAGGGACT  >  NC_000913/3581072‑3581560

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