Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,816,349	T→C	P67P (CCA→CCG)	rph ←	truncated RNase PH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,816,349	0	T	C	98.3%	203.1 / NA	60	P67P (CCA→CCG)	rph	truncated RNase PH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (22/37);  minor base A (1/0);  total (23/37)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCGGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑AAACAAATGTGGCTGCGCATTATACGGA  >  NC_000913/3816105‑3816595                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        tGCCGTTTTCCACCAGCTTCTGTAGCGCATCTGCCAGCGCCACACAGGCACCCGTAATCGAAGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGTGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCTGCTTCACGAGCGTTACGGGTGTGGGTAGAACGCGgcagc                                                                                                                                                                                                                                                   >  1:190173/1‑251 (MQ=255)                                       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)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        TGCCGTTTTCCACCAGCTTCTGTAGCGCATCTACCAGCGCCACGCAGGCACCCGTAATCGACGCGGTACGCGTGCCACCATCAGCCTGAAGCACGTCGCAGTCCAGCGTAATGGTGAACTCACCCAGCGCTTTCAAATCTACTGCCGCGCGAAGAGCACGGGCGATCAGACGCTGGATTTCCATTGTGCGTCCACCCTGCTTACCTTTCGCCGCTTCACGAGCGTTACGGGTGTGGGTAGAACGTGGCAGCATGCCGTACTCTGCGGTGATCCAGCCCTGGCCCTGACCTTTCAGGAAGCGCGGCACGCCTTCTTCAATAGAGGCGGTACACAACACTTTGGTATCGCCAAATTCGACCAGCACCGAGCCTTCTGCATGTTTTGTATAGTTACGAGTCAGGGTAACGGGACGCACCTGATTATTGCTACGGCCTGCTGGACGCATATTGAAATCTCCGGCTTG‑AAACAAATGTGGCTGCGCATTATACGGA  >  NC_000913/3816105‑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: ‑