Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,816,148	G→A	C134C (TGC→TGT)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,816,148	0	G	A	97.2%	225.4 / ‑3.8	72	C134C (TGC→TGT)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (42/27);  total (44/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

GAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACC  >  NC_000913/3815905‑3816393                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      gAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCAc                                                                                                                                                                                                                                                <  1:262726/251‑1 (MQ=255)   gTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCAccc                                                                                                                                                                                                                                              <  2:306623/251‑1 (MQ=255)   gTAGCCCTTCATGGGTGAACCGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGCTGCGCACGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCAccc                                                                                                                                                                                                                                              >  2:222902/1‑251 (MQ=255)     aGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGt                                                                                                                                                                                                                                            <  1:211593/251‑1 (MQ=255)        tctTCATGGGTGAACGGCTCGCCTTCTGCCGTACCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCTCCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAAt                                                                                                                                                                                                                                         >  1:136764/1‑251 (MQ=255)         ctTCATGGGTGACCGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTaa                                                                                                                                                                                                                                          <  2:358316/249‑1 (MQ=255)          ttCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCGACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCg                                                                                                                                                                                                                                       >  1:164877/1‑251 (MQ=255)          ttCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAAt                                                                                                                                                                                                                                         <  1:50134/249‑1 (MQ=255)          ttCATGGGTGAACGGCTCGCCTTCTGCAGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCCCACAGGCACCCggcat                                                                                                                                                                                                                                         >  2:50134/1‑245 (MQ=255)             aTGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGaag                                                                                                                                                                                                                                    <  1:213632/251‑3 (MQ=255)                ggTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCgg                                                                                                                                                                                                                                 >  2:354385/1‑251 (MQ=255)                             ccTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACCATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGaagc                                                                                                                                                                                                                                   >  1:41527/1‑233 (MQ=255)                             ccTTCTGCCGTCCCCGTCACTGCAATGAGGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCCCCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCGGTAGCGCATCTGCCATCGCCACACAGGCACCCGTatgcgaagc                                                                                                                                                                                                                                   <  2:41527/236‑9 (MQ=255)                                                cTTCAATGATGCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCACATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCAAAGCGGTACGCGTGCCACCATCAGCCTGACGCACGTCg                                                                                                                                                                                                 >  2:207745/1‑251 (MQ=255)                                                       gatgCGCCCGTCTTCGGTCATCCCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTGTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCAGCAGCCTGAAGCACGTCGCAGTCCa                                                                                                                                                                                          <  2:438059/251‑1 (MQ=255)                                                         tgCGCCCTTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCTTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAg                                                                                                                                                                                         <  1:29982/250‑1 (MQ=255)                                                            gcCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGt                                                                                                                                                                                      >  2:90046/1‑250 (MQ=255)                                                            gcCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTa                                                                                                                                                                                     <  1:180941/251‑1 (MQ=255)                                                                  cTTCGGTCATCCCTAAGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTAACATCCCTTTCATTGGATTGGTTTTCCGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCCCCATCAGCCCGACGCACGTCCCCGTCCAGCGTAATGGTg                                                                                                                                                                               >  1:359794/1‑251 (MQ=255)                                                                     cGGTCATCACTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAAc                                                                                                                                                                            <  1:32142/251‑1 (MQ=255)                                                                           tcaCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGCTCGCAAATTGCTTCACCGTTTACAATAACGACAGAAACCGCGGCTACCATCCCTTTCCTTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCAccc                                                                                                                                                                      >  1:431867/1‑251 (MQ=255)                                                                             aCTACGTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCAccc                                                                                                                                                                      >  1:98349/1‑249 (MQ=255)                                                                                  gTTCATGTCGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCTCGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCtt                                                                                                                                                               <  2:39345/251‑1 (MQ=255)                                                                                          cGGTCTCTGCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCTACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTTCCACCATCAGCCTGAAGCACGTCGCAGGCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATc                                                                                                                                                        <  2:261527/250‑1 (MQ=255)                                                                                                  gCGGCTGAGTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACGCACCCAGCGCTTTCAAATCTACTGCcg                                                                                                                                                <  2:43264/250‑1 (MQ=255)                                                                                                       tgagTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGCAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTCATGGTGAACTCACCCAGCGCTTTCAAATCTCCTGCCGCGCGa                                                                                                                                           >  2:75994/2‑250 (MQ=255)                                                                                                       tgagTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGcca                                                                                                                                                                                                                      <  2:41135/174‑1 (MQ=255)                                                                                                       tgagTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGcca                                                                                                                                                                                                                      >  1:41135/2‑175 (MQ=255)                                                                                                       tgagTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAACGCTTTCAAATCTACTGCCGCGCGa                                                                                                                                           >  1:261153/2‑250 (MQ=255)                                                                                                          gTCTTCAACGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAagtg                                                                                                                                       >  2:85361/1‑249 (MQ=255)                                                                                                                 aCGTATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGCAGTGCACGGGc                                                                                                                                >  2:200923/1‑251 (MQ=255)                                                                                                                    tATTCCAGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGAt                                                                                                                             >  1:243691/1‑251 (MQ=255)                                                                                                                          aGATCGCAAATTGCTTCACCGTTTACAATACCGACAGAACCCGCGGCTACCATCCCTTTCATTGGATTGGTTGTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGatg                                                                                                                       <  2:302192/251‑3 (MQ=255)                                                                                                                            aTCGCAAATTGCTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCt                                                                                                                     <  1:396126/251‑1 (MQ=255)                                                                                                                                      gCTTCACCGTTTACAATACAGACAGAAACCCCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCAGTTTTCCAACAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGCAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCAt                                                                                                           >  1:159969/1‑251 (MQ=255)                                                                                                                                       cTTCACCGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCACATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCAtt                                                                                                          >  1:151756/1‑251 (MQ=255)                                                                                                                                          caccGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATtgtg                                                                                                       <  1:184819/249‑1 (MQ=255)                                                                                                                                           accGTTTACAATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTATAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGAGGTACGCGTGCCACCATCGGCCTGACGCACGTCGCAGTCCAGAGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCCAAGTGCACGGTCGAACAGCCGCTGGTTTTCCAttacgc                                                                                                      >  1:430867/2‑247 (MQ=255)                                                                                                                                                 tacaATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCCGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCAc                                                                                                >  2:352252/2‑251 (MQ=255)                                                                                                                                                  acaATACCGACAGAAACCGCGGTTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCATACGCTGGATTTCCATTGTGCGTCCAcc                                                                                               >  1:139371/1‑251 (MQ=255)                                                                                                                                                  acaATACCGACAGAAACCGAGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGCGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCAcc                                                                                               >  2:49218/1‑251 (MQ=255)                                                                                                                                                   caATACCGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCCCGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCACATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCAccc                                                                                              >  2:193784/1‑251 (MQ=255)                                                                                                                                                     ataccGACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGGTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTg                                                                                            <  2:90908/248‑1 (MQ=255)                                                                                                                                                          gACAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTAAGCTTGCCGTTTTCCACCAGCTTATGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAAATCACACAGCGCTTTCAAATCTACTGCCGCGAGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCttc                                                                                        >  1:279953/1‑249 (MQ=255)                                                                                                                                                           aCAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTAcc                                                                                      <  1:85361/251‑1 (MQ=255)                                                                                                                                                           aCAGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTAcc                                                                                      >  1:51288/1‑251 (MQ=255)                                                                                                                                                             aGAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCTCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCtt                                                                                    <  2:167034/251‑1 (MQ=255)                                                                                                                                                              gAAACCGCGGCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCttt                                                                                   >  2:301654/1‑251 (MQ=255)                                                                                                                                                                      ggCTACCATCCCTTTCATTGGATGGGTTTTCAGCTTTCCGTTTTCCACCAGCATCTGTAGCCCATATGCCAGCGCCACACAGGCACCCGTAATCGCAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCTAGCGTAATGGAGAACTCACCCAGCGCTTTCAAATCTACCGCTGCACGAGGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTACACCCTGCTTACCTTTCGCTGCtt                                                                           >  2:130871/1‑251 (MQ=255)                                                                                                                                                                       gCTACCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCACGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTc                                                                          >  1:399401/1‑251 (MQ=255)                                                                                                                                                                          aCCATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTTTGTAGCGCATCTGCCAGCGCCACTCAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGCAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACg                                                                       >  2:59546/1‑251 (MQ=255)                                                                                                                                                                          aCCATCCCTTTCATTGGATTGGTGTTCAGCGTGCCGTTTTCCTCCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGGGCCACCATCAGCCTGAAGCACGTCGCAGTCCGGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCCCGAAGGGCACGTGCGATCAGACGCTGGATGGCCATTGTGCGTCCACCCTGCTTTCCTTTCGCTGCATCACg                                                                       <  2:431867/251‑1 (MQ=255)                                                                                                                                                                            cATCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACg                                                                       >  2:134325/1‑249 (MQ=255)                                                                                                                                                                             aTCCCTTTCATTGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGc                                                                    >  1:208247/1‑251 (MQ=255)                                                                                                                                                                                   atcattgcaTTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACCCGCGCCACCATAAGCCTGAAGCACGTCGCAGTCCCGCGTATTGGTGAACTCAACCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCGTTCTGCGTCCACACTGCTTACCTTCCGCTGCTCCACGCGCGTTACg                                                              >  2:275140/9‑251 (MQ=255)                                                                                                                                                                                     cattggATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACgg                                                             <  1:356388/246‑1 (MQ=255)                                                                                                                                                                                     cattggATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGCATCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACTTCGCAGTCCCGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCTCACGTAGTGCACGGGCGCTCAGACGCTGGCTTTCCATTGGGCGTCCACCCTGCTTACCTTTCGCTGCTTCCCGATCGTTACggg                                                            >  2:232991/5‑251 (MQ=255)                                                                                                                                                                                      attggATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGAATCTGCCAGCGCCAAACATGCACACGTAATCGAAGCAGTACGCGTGCCACCATCCGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACATTTCGCTGCTACACGAGCGTTACGGgt                                                           >  2:182279/4‑251 (MQ=255)                                                                                                                                                                                         ggATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGtgtg                                                        <  1:67286/251‑1 (MQ=255)                                                                                                                                                                                          gATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTgg                                                       >  1:289358/1‑251 (MQ=255)                                                                                                                                                                                            ttGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGt                                                     >  1:42037/1‑251 (MQ=255)                                                                                                                                                                                              ggTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGCAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCCGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGt                                                     >  2:278577/1‑249 (MQ=255)                                                                                                                                                                                                  ttCAGCTTGCCGTTTTCCACCAGCTGCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTACTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAAc                                                <  2:419077/250‑1 (MQ=255)                                                                                                                                                                                                    cAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGAt                                                                                                                             >  1:188296/1‑171 (MQ=255)                                                                                                                                                                                                    cAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGAt                                                                                                                             <  2:188296/171‑1 (MQ=255)                                                                                                                                                                                                    cAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACgcg                                             >  1:176700/1‑249 (MQ=255)                                                                                                                                                                                                     aGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACgcgg                                            >  1:127309/1‑248 (MQ=255)                                                                                                                                                                                                      gcgtGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGCAGCTGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTTAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCGCCCTGCTTACCGTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGatcgcggc                                           <  2:210769/248‑8 (MQ=255)                                                                                                                                                                                                         tGCCGTTTTCCACCATCTTCGGTAGCGCATCTGCCAGCGCCACGCAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCAcc                                                                                                                                                                       <  1:261137/124‑1 (MQ=255)                                                                                                                                                                                                         tGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGgcagc                                        >  1:190173/1‑251 (MQ=255)                                                                                                                                                                                                         tGCCGTTTTCCACCAGATTCTGTAGCGCATCTGCCAGCGCCACGACGGCACCCGTCATCGCAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCAcc                                                                                                                                                                       >  2:261137/1‑124 (MQ=255)                                                                                                                                                                                                                                               gcCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGGCAGCATACCGTACTCTGCGGTGATCCAGCCCTGGCCCTGAcc  >  2:428132/1‑251 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      GAGTAGCTCTTCATGGGTGAACGGCTCGCCTTCTGCCGTCCCCTGCACTTCAATGATGCGCCCGTCTTCGGTCATCACTACGTTCATGTCGGTCTCTGCGGCAGAGTCTTCAACGTATTCCAGATCGCAAACCGCTTCGCCGTTCACAATTCCGACAGAAACTGCGGCTACCATCCCTTTCATCGGATTGGTTTTCAGCTTGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACC  >  NC_000913/3815905‑3816393

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: ‑