Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,816,358	G→A	G64G (GGC→GGT)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,816,358	0	G	A	100.0%	202.0 / NA	59	G64G (GGC→GGT)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/37);  total (22/37)

GCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCGGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑AAACAAATGTGGCTGCGCATTATACGGA  >  NC_000913/3816143‑3816595                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               gcCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGAcc                                                                                                                                                                                                             >  2:428132/1‑251 (MQ=255)   cacacAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCtt                                                                                                                                                                                                           >  1:79537/5‑251 (MQ=255)    acacAGGCACCCGTAATCGAAGCGGTACGCGGTCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCGGGCCCTGACCttt                                                                                                                                                                                                          <  2:208247/248‑1 (MQ=255)        aGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAgg                                                                                                                                                                                                      >  1:140711/1‑251 (MQ=255)         ggCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGa                                                                                                                                                                                                     <  1:354385/251‑1 (MQ=255)          gCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGaa                                                                                                                                                                                                    <  1:360591/251‑1 (MQ=255)             cccGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCATTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCg                                                                                                                                                                                                                                    <  2:251066/216‑1 (MQ=255)             cccGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCg                                                                                                                                                                                                                                    >  1:251066/1‑216 (MQ=255)             cccGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAgcg                                                                                                                                                                                                 <  2:151756/251‑1 (MQ=255)                  aaTCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCa                                                                                                                                                                                            <  1:86489/251‑1 (MQ=255)                      gaagcagTACGCGTGCCACCATCCGCCTGACGCACGTCGCAGTCCAGCGTACTGGTGCACTCACCCAGCGCTTTCAAATATACTGCCGCGCGAATTGCACGGGCGATCAGACGCTGGATTTACCTTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGAGGCAGCATACCGGCCTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGCAGCGAGGCCCGcc                                                                                                                                                                                        >  2:186850/7‑251 (MQ=255)                          cgGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGCGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCcttct                                                                                                                                                                                    >  1:29974/1‑251 (MQ=255)                                    gCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATagag                                                                                                                                                                           <  1:292008/250‑1 (MQ=255)                                      caccaTCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCg                                                                                                                                                                        <  1:428132/251‑1 (MQ=255)                                                  gAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACacaatac                                                                                                                                                             <  2:254615/250‑4 (MQ=255)                                                           cGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCg                                                                                                                                                   <  1:403903/251‑1 (MQ=255)                                                              aGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGcc                                                                                                                                                 >  1:306125/1‑250 (MQ=255)                                                                  cAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTAAACAATACTTTGTTATCGCCAAAtt                                                                                                                                            >  2:240380/1‑251 (MQ=255)                                                                      gTAATGGTGAAATCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTATTCAATAGAGGCGGTACACAATACTGTGGGATCGCCAAATTctac                                                                                                                                        >  1:88851/1‑248 (MQ=255)                                                                        aaTGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGAcc                                                                                                                                       <  1:200923/250‑1 (MQ=255)                                                                            gTGACCTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTAAACAATACTTTGGTATCGCCAAAttttcccagccc                                                                                                                                  >  2:35267/1‑241 (MQ=255)                                                                                aCTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGTCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTGTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAg                                                                                                                              <  2:136764/251‑1 (MQ=255)                                                                                              ttCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGtt                                                                                                                 <  1:222902/250‑1 (MQ=255)                                                                                                cAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATGTCCATTGTGCGGCCACCCGGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGtttt                                                                                                               <  2:139371/250‑1 (MQ=255)                                                                                                       aCTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTTCCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTCCTGCATGTTTTGTATAg                                                                                                         <  2:325160/249‑1 (MQ=255)                                                                                                       aCTGCCGAGCGAAGTGCACGGGCGACCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACAGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATTCTTTGGTATCGCCAAATTCGACCATCACCGAGCCTTCTGCATGTTTTGTATAg                                                                                                         <  2:7044/249‑1 (MQ=255)                                                                                                         tGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACCCAATACTTTGGTATCGCCAAATTCGACCCACCCCGCGCCTTCTGCCTGTTTTTTATAGttcc                                                                                                     >  2:30754/1‑249 (MQ=255)                                                                                                          gCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTGCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGTTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTGGTATAGTTAc                                                                                                     <  2:164877/250‑1 (MQ=255)                                                                                                                 gAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTAACCAATACTTTGTTATCTCCAAAATCGACCCGCACCCAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAgg                                                                                             >  2:80745/1‑251 (MQ=255)                                                                                                                       cACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCTCCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACg                                                                                       <  2:243691/251‑1 (MQ=255)                                                                                                                        aCGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTTGACCCGCACCGAGCCCTCTGCATGTTTTGGATTGTTTCGAGTAAGGGTAACgg                                                                                      >  2:120068/1‑251 (MQ=255)                                                                                                                         cGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACggg                                                                                     >  1:299716/1‑251 (MQ=255)                                                                                                                                    aCGCTGGATTTCCATTGTGCGTCCACCCTGCTTCCCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGat                                                                          <  2:289358/251‑1 (MQ=255)                                                                                                                                    aCGCTGGATTTCCATTGTGCGTCCACCCTGCTGACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGat                                                                          <  1:402223/251‑1 (MQ=255)                                                                                                                                            tttCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCATGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGGAACGGGACGCACCTGATTATTGCt                                                                   >  1:182864/1‑250 (MQ=255)                                                                                                                                             ttCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTGTGTATAGTTACGAGTCAGGGTAACGGGACGCGCCTGATTAGTGCTa                                                                  <  2:190173/250‑1 (MQ=255)                                                                                                                                             ttCCATTGTGCGGCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACGCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCATCACCGAGCCTTCTTCATGTTTTGTATAGTTACGAGTCAGGGTAACGGTACGCATCTGATGAGTGCTa                                                                  <  2:399401/250‑1 (MQ=255)                                                                                                                                              tCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAggg                                                                                            <  1:133278/223‑1 (MQ=255)                                                                                                                                               ccATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACgg                                                               <  1:90046/251‑1 (MQ=255)                                                                                                                                                   tgtgCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCct                                                            <  1:49218/250‑1 (MQ=255)                                                                                                                                                   tgtgCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCct                                                            <  1:301654/250‑1 (MQ=255)                                                                                                                                                   tgtgCGTCCACCCGGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCGACGGCct                                                            <  2:42037/250‑1 (MQ=255)                                                                                                                                                        gTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGGATCGCCAAATTCGGCCAGCACCGGGCCTTCTGCATGTTTTGTATAGTTACGAGTAAGGGTAAAGGGAAGCACCTGATTATTGCTACGGCCTGCTggg                                                      >  1:129834/1‑250 (MQ=255)                                                                                                                                                          ccACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGGATCGCCAAATTTGACCAGCACCCAGCCTTCTGCCTGTTTTGTATAGTTACGGGTCCGGGTAAAGGGCCGCACCTGATTATTGCTACGGCctgcctggcg                                                    >  1:17819/1‑245 (MQ=255)                                                                                                                                                              ccTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCata                                                <  2:176700/251‑1 (MQ=255)                                                                                                                                                              ccTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCata                                                <  1:75994/251‑1 (MQ=255)                                                                                                                                                                 gCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCTCCTGATTATTGCTACGGCCTGCTGGACGCATAtt                                              <  1:372273/250‑1 (MQ=255)                                                                                                                                                                 gCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCatat                                               <  1:35267/249‑1 (MQ=255)                                                                                                                                                                     aCCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGaaa                                          <  1:134325/250‑1 (MQ=255)                                                                                                                                                                             ctgcTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACATTGGTATCGACAAATTCGACCAGCACCGAGCCTTCTGGATGTTTTGTATAGTTACGAGTCAGGGTAACGGGGCGGACCTGATTATTGCTACGGCCTGCTGGACGCAAATTGAAATCTCCGGCt                                 >  1:20466/3‑251 (MQ=255)                                                                                                                                                                                 ttCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑a                             >  1:331027/1‑250 (MQ=255)                                                                                                                                                                                    aCGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTCCCCAAAACTTTGGTATTGCCCAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGGGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGATGCATATTGAAATCTCCGGCTTG‑AAAc                          >  1:148260/1‑250 (MQ=255)                                                                                                                                                                                                ggTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATGATTGCTACGGCCTGCTGGACGCAGATTGAAATCTCCGGCTTG‑AAACAAATGTGGCTgcg             <  2:29974/251‑1 (MQ=255)                                                                                                                                                                                                      ggTAGAACGCGGCAGCATACCGGACTCTGCGGTGATCCAGCCATGGCCCTGTCCTTTCAGGAAGCGAGGCTCGCCTTCTTCAATAGAGGCTGTACACACTACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTCCGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGAATATTGAAATCTCCGGCTTG‑AAACAAATGTGGCGGCGCATtat       <  2:235132/251‑1 (MQ=255)                                                                                                                                                                                                          gAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑AAACAAATGTGGCTGCGCATTATACgg   <  2:51288/251‑1 (MQ=255)                                                                                                                                                                                                          gAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCCAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCCCCTGATTATTGCTACGGCCTGCTGGACGCATATTGGAATCTCCGGCTTG‑AAACAAATGTGGCTGCGCATTCTACgg   >  1:184155/1‑251 (MQ=255)                                                                                                                                                                                                           aaCGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGGGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCCCCTGATTATTGCTACGGCCTGCTGGGCGCATATTGAAATCTCCGGCTTGTAAACaat                       >  2:72180/1‑230 (MQ=255)                                                                                                                                                                                                           aaCGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑AAACAAAt                      <  1:72180/231‑1 (MQ=255)                                                                                                                                                                                                           aaCGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGAGGCACGCCTTCTTCAATAGAGGCGGTACACAATACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑AAACAAATGTGGCTGCGCATTATACGGa  <  1:278577/251‑1 (MQ=255)                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               GCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCGGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑AAACAAATGTGGCTGCGCATTATACGGA  >  NC_000913/3816143‑3816595

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