Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	155,463	A→T	*247K (TAA→AAA)	yadV ←	putative fimbrial chaperone YadV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	155,463	0	A	T	98.5%	245.1 / NA	68	*247K (TAA→AAA)	yadV	putative fimbrial chaperone YadV
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (29/38);  minor base G (0/1);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCTTCTGAACCTGACCAG  >  NC_000913/155217‑155707                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     tACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGc                                                                                                                                                                                                                                                  >  1:727226/1‑251 (MQ=255)  aCATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGATGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGCTGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGcc                                                                                                                                                                                                                                                 >  2:724489/1‑251 (MQ=255)    aTAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTg                                                                                                                                                                                                                                               >  1:568882/1‑251 (MQ=255)      aaaCACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                                                                                                                                                                                                              <  2:122719/250‑1 (MQ=255)        acacTGACGTCGTAAACACCCGGTATAGCGGNGTTACCGTCACTGGACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAt                                                                                                                                                                                                                                           <  2:340024/251‑1 (MQ=255)          acTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAtt                                                                                                                                                                                                                                          >  1:464421/1‑250 (MQ=255)           cTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTAc                                                                                                                                                                                                                                        >  2:223371/1‑251 (MQ=255)                   gtAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAAtt                                                                                                                                                                                                                                >  1:725573/1‑251 (MQ=255)                     aaaaacCCGGTATAGACGGGTTTGCTTCACCGTCCCTGCTGAGATCAATATTACATGCATTCTGCCCCATCAGAAAGGTATGTTCATATTCAACGAGTTCAGCACTGACAGCAGTATTGTAATACATCAGCGCGCAAAACGTGACGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATATTCACGCTATTCTTACTGCGTAATGATGTCCTTTCCTGATTTCTTCCTTGCATTACCTTCAATTgc                                                                                                                                                                                                                              <  2:418468/247‑1 (MQ=255)                      aaCACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCTCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTgcg                                                                                                                                                                                                                             <  1:724489/251‑1 (MQ=255)                               tataGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCaa                                                                                                                                                                                                                      >  1:558671/1‑249 (MQ=255)                                  aGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTc                                                                                                                                                                                                                  >  1:113462/1‑250 (MQ=255)                                    cGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCAtt                                                                                                                                                                                                               >  2:620635/1‑251 (MQ=255)                                     ggtgTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCCGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATTATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTCCCTTCAATTGCGCCACCAAAGTCATTa                                                                                                                                                                                                              <  2:555299/248‑1 (MQ=255)                                       ggTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAAt                                                                                                                                                                                                            <  2:23485/251‑1 (MQ=255)                                        gTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATg                                                                                                                                                                                                           >  1:177326/1‑251 (MQ=255)                                                cACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTaaa                                                                                                                                                                                                    >  2:39579/1‑250 (MQ=255)                                                   tGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATg                                                                                                                                                                                                >  1:28965/1‑251 (MQ=255)                                                                         aGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTg                                                                                                                                                                          <  2:584386/251‑1 (MQ=255)                                                                          gATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGc                                                                                                                                                                         <  1:620635/251‑1 (MQ=255)                                                                                    gCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAgg                                                                                                                                                                 <  1:711567/249‑1 (MQ=255)                                                                                         aTCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                                                                                                                                                                                                              >  1:662905/1‑167 (MQ=255)                                                                                         aTCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTgg                                                                                                                                                                                                                                              <  2:662905/167‑1 (MQ=255)                                                                                           cAGGAAGGTATGTTCATATTCAAAGAGTTCAGCACTGAAAGCAGGATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTTCCTTCAATTGCGCCACCAAAGTCATTAAGGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACt                                                                                                                                                        <  2:727226/251‑1 (MQ=255)                                                                                           cAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGAc                                                                                                                                                         <  1:554593/250‑1 (MQ=255)                                                                                            aGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACtt                                                                                                                                                       <  1:560028/251‑1 (MQ=255)                                                                                                aGGTATGTTCATATTCAACGAGTTCAGCACTGAAAGCATTATTGCAATACAGCATCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCCTTGCCATTAAGGCCATTGACTTTCAt                                                                                                                                                   <  2:526720/251‑1 (MQ=255)                                                                                                aGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACGTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCAt                                                                                                                                                   <  2:725573/251‑1 (MQ=255)                                                                                                 ggTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATNCGCTTTGCCATTAAGGCCATTGACTTTCATg                                                                                                                                                  >  2:390573/1‑251 (MQ=255)                                                                                                    aTGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTTCGGAATTATATCCTTGCCTGATTGCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATGCGCTTTGCCATTAAGGCCATTGACTTTCATGAcc                                                                                                                                               <  2:217478/251‑1 (MQ=255)                                                                                                      ggTCATATTCAACGAGTTCAGCACTGAAAGCAATATTGCACTACATCAGCGCGAACAACGTGGCGATAAGGGACAGATGATGATCATTTTTCGTATCTTATCTAGTCACGCTATACATCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTTAAATTACCTTCAATTGCGCCAACTCAGTCATTAATGGCGTACAACTGCACTTTTGCAGAAATCGCTTTGACATTAAAGCCATATCCTTTCATGACCtc                                                                                                                                             >  2:718503/1‑251 (MQ=255)                                                                                                         cATATTCAACGAGTTCAGCGCTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGTTCATATGATGATAAATTTTAGTATATTCTATAGTCCCGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAAATCGCTTTGCCATTAAGGCCATTGACTTTCATGACCtcatc                                                                                                                                          <  2:219321/251‑1 (MQ=255)                                                                                                             ttCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATAACt                                                                                                                                       >  2:607549/1‑250 (MQ=255)                                                                                                                  cGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATgg                                                                                                                                 <  1:598442/251‑1 (MQ=255)                                                                                                                       tCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATAGTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCTTTAATGGCGTAAAAATGCACTTTTGCATAATGCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCa                                                                                                                             <  1:726055/250‑1 (MQ=255)                                                                                                                           cACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCAt                                                                                                                        <  1:629265/251‑1 (MQ=255)                                                                                                                              ggaaagcattatggcgttaaagcagcGTGCATAACGTGGCTATTCGGATCAGATGATGATAATTGTTATGATATTCTATAGTCACGCTATGCTTCCTTCGGAATTATATCCGTGGCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAATTCCTTAATGGCGTAAAAATTCACTTTTGCAGTATTCGCTTTGCCATTAAGGCCATTGACTGTCATAACCTCATCATTAAATGGTGCcatcatgtt                                                                                                                     <  2:509807/231‑4 (MQ=255)                                                                                                                                aaaGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGCCCtcat                                                                                                                                           >  2:279159/1‑227 (MQ=255)                                                                                                                                aaaGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCtcat                                                                                                                                           <  1:279159/227‑1 (MQ=255)                                                                                                                                 aaGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCATATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTcac                                                                                                                  <  2:406468/251‑1 (MQ=255)                                                                                                                                        aTTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATc                                                                                                           <  1:693129/251‑1 (MQ=255)                                                                                                                                          tGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTGTATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCg                                                                                                          <  1:747528/250‑1 (MQ=255)                                                                                                                                              aTACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAg                                                                                                     <  1:438368/251‑1 (MQ=255)                                                                                                                                                    cagcGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACGAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTa                                                                                               <  2:464421/251‑1 (MQ=255)                                                                                                                                                       cgcgcacgaCGTCGCGATACGGGTCAGATGATGATAATTTTTAGTAGATTCTATAGTCACGATATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCAGTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACGTTTGCAGACTTCGCCTTGCCATTAAGGCCATTGACTTTCCTGACCTCATCACTAAATGGTGCAATCATTTTCACATCACTCGGATAGCGTTTAcct                                                                                            >  2:178555/9‑250 (MQ=255)                                                                                                                                                       cgcgcAAAACGTGGCGATCCGGGTCAGATGATGATAATTTTTAGGATATTCTATAGTCACGCTATTCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTGTGCCATTAAGGCCATTGACTTTCAGGACCTCATCACGAAATGGTGCAATCAGTTTCACATCAATCGGATATCGTTTACCg                                                                                            <  2:77572/251‑1 (MQ=255)                                                                                                                                                              aaCGTGGCGATACGGGTCAGATGATTATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCtt                                                                                     <  2:408778/251‑1 (MQ=255)                                                                                                                                                                       aTACGGGTCAGATGATGATAATTTGTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATcac                                                                            <  2:276158/251‑1 (MQ=255)                                                                                                                                                                           gggTAAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCATGCGGAATTATATCCTTGCATGATTTCAGCCTGGCATTACCTTAAATTGCGCAACCAAAGTCATTAATGGAGTAAAACTGCACTTTTGCAGACTACGATATGCCCTCAAGACGAATGACTTTCTGGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCTTTTaacacgatcctacaattaaccaat                                                                        >  2:471226/1‑228 (MQ=255)                                                                                                                                                                             gTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTa                                                                                  <  2:415137/239‑1 (MQ=255)                                                                                                                                                                             gTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTa                                                                                  >  1:415137/1‑239 (MQ=255)                                                                                                                                                                               cAGATGATGAAAATTTTTAGTATATTCTATAGTCACGCCATGCTTCCTGCGGATTTATATCCTTACCTGTTTTCAGCCTGGCATTACCTTCAATTGCGCCACAAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCTTTGACTTTAATGACCTCATAACTAAATGGTGCAATCCTTTTCACATCAATCGGATAGCGGTTACCGCTAGCTTCTAACTCACCACTGCt                                                                     >  2:244495/1‑250 (MQ=255)                                                                                                                                                                                          aaTTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACgt                                                          >  1:198205/1‑250 (MQ=255)                                                                                                                                                                                            tttttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACgtagt                                                       >  1:615620/1‑251 (MQ=255)                                                                                                                                                                                              tttAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTCCTAAATCACCACTGCTAAAAGAGACgtagta                                                      >  1:196062/1‑250 (MQ=255)                                                                                                                                                                                                         cTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGAtt                                          >  1:678134/1‑251 (MQ=255)                                                                                                                                                                                                          tataGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCACTCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGCGACGTAGTAAGGGGTTGGATTg                                         >  2:18229/1‑251 (MQ=255)                                                                                                                                                                                                                cACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACt                                   <  2:568882/251‑1 (MQ=255)                                                                                                                                                                                                                  cGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCg                                 <  1:39579/251‑1 (MQ=255)                                                                                                                                                                                                                   gCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGt                                >  1:401345/1‑251 (MQ=255)                                                                                                                                                                                                                      aTGCTTCCTGCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACGGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAAt                             <  1:216592/251‑1 (MQ=255)                                                                                                                                                                                                                               gCGGAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCtt                     <  2:139112/250‑1 (MQ=255)                                                                                                                                                                                                                                  gAATTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCt                 <  2:678134/251‑1 (MQ=255)                                                                                                                                                                                                                                   aaTTATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCtt                <  1:562880/251‑1 (MQ=255)                                                                                                                                                                                                                                     ttATATCCTTGCCTGATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCTTCt              >  1:92349/1‑251 (MQ=255)                                                                                                                                                                                                                                             ttGCCTGATTTCAGCCTGGCATTACCTTCAATTGAGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCTTCTGAACCTg       >  1:169217/1‑250 (MQ=255)                                                                                                                                                                                                                                               gCCTGATTTCAGCCTGGCATGACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCTTCTGAACCTGAc     <  2:432610/250‑1 (MQ=255)                                                                                                                                                                                                                                                   gATTTCAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCTTCTGAACCTGACCAg  <  2:401345/249‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     TACATAAACACTGACGTCGTAAACACCCGGTATAGCGGGGTTACCTTCACTGTACCGGCTGAGATCAATATTAGATGCATTCTGCCCCATCAGGAAGGTATGGTCATATTCAACGAGTTCAGCACTGAAAGCAGTATTGCAATACAGCAGCGCGCAAAACGTGGCGATACGGGTCAGATGATGATAATTTTTAGTATATTCTATAGTCACGCTATGCTTCCTGCGGAATTATATCCTTGCCTGATTACAGCCTGGCATTACCTTCAATTGCGCCACCAAAGTCATTAATGGCGTAAAAATGCACTTTTGCAGAATTCGCTTTGCCATTAAGGCCATTGACTTTCATGACCTCATCACTAAATGGTGCAATCATTTTCACATCAATCGGATAGCGTTTACCGCTAGCTTCTAAATCACCACTGCTAAAAGAGACGTAGTAAGGGGTTGGATTGGTCACTCGTAATGACGCCTTACCTTCTGAACCTGACCAG  >  NC_000913/155217‑155707

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