Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,821,992	C→T	R189C (CGC→TGC)	gmk →	guanylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,821,992	0	C	T	100.0%	232.3 / NA	69	R189C (CGC→TGC)	gmk	guanylate kinase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (42/27);  total (42/27)

GGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCTACGCGGTCGCGGTCAGGACAGCGAAGAGGTCATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGCGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGTCGCGCTCGTCAGATGCAGGTAGGCGGAAAGGATCCGCTGGTAC  >  NC_000913/3821745‑3822223                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                        gggCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCTACGCGGTCGCGGTCAGGACAGCGAAGAGGTCATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCTCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCa                                                                                                                                                                                                                                      <  1:375555/251‑1 (MQ=255)       gcagcaAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCTACGCGGTCGCGGTCAGGACAGCGAAGAGGTCATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAgcc                                                                                                                                                                                                                                <  2:25359/251‑1 (MQ=255)            aaaTTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCTACGCGGTCGCGGTCAGGACAGCGAAGAGGTCATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCaa                                                                                                                                                                                                                           <  1:364573/251‑1 (MQ=255)                  gCCAGATGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAAATGAACTGGACCGCCGTCTACGCGGTCGCGGTCAGGACAGCGAAGAGGTTATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAagcag                                                                                                                                                                                                                     <  2:303200/251‑1 (MQ=255)                    aagaagaTGCCGCTCTCGTGGAGTATCTTTATGATACCGCCGTCCAAAAATGAACTGGACGTCCGTCTACGCTGTCGCGGTCAGGACAGCGAAGAGGTCATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTTCGCCGTATATGATTATCTTATTGTGAATGATGACTTCGATACCGCGTTGACCTATTTTAAGATCATTAGTCGTGCCGAACTTCTGTGCAGGAgccgcttaaagaagcg                                                                                                                                                                                                                   <  2:229057/250‑12 (MQ=255)                     agaagaTGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCTACGCGGTCGCGGTCAGGACAGCGAAGAGGTCATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAACAGCAGCGt                                                                                                                                                                                                                  >  2:356750/1‑251 (MQ=255)                      gaagaTGCCGCACGCGCGGAGTTTCATTATTTGACCGCCGTCCAAAATTGAACTGGACCGCCGGCTACGCGGTCGCGGGCAGGACAGCGAAGAGGTCATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTGGACCGATTTGAAGTCCAGTATTCGCGCCGATCGTCTGTGCATGAGCCGCCAAAAGCAGCGTc                                                                                                                                                                                                                 <  2:336503/251‑1 (MQ=255)                            gCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCTACGCGGTCGCGGTCAGGACAGCGAAGAGGTCATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACg                                                                                                                                                                                                           >  1:163753/1‑251 (MQ=255)                                      cgGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCTACGCGGTCGCGGTCAGGACAGCGAAGAGGTCATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATc                                                                                                                                                                                                   >  2:184967/1‑249 (MQ=255)                                              cTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCTACGCGGTCGCGGTCAGGACAGCGAAGAGGTCATTGCACAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAAttgt                                                                                                                                                                                         >  2:148798/1‑251 (MQ=255)                                              cTTTATTTTACCGCCGACCAAAATTGTACTGGACCGCCGTCTACGCGGTCGCGGTCAGGACAGCGAAGCGGTCATTGCAAAGCGAATGGCGCAAGCTGTTGCAGAAATGAGCCATTACTCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGAGCCGAACGTCTGTGCATAAGCCGACAAAAGCAGCGTCCTGACGCTTTAATCAGAAAAttgt                                                                                                                                                                                         >  1:417114/1‑251 (MQ=255)                                                ttATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCTACGCGGTCGCGGTCAGGACAGCGAAGAGGTCATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATGCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAAttgtt                                                                                                                                                                                        <  1:245674/250‑1 (MQ=255)                                                 tATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCTACGCGGTCTCGGTCAGGACAGAGAAGAGGTCATTGCAAAGCGTATGACGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATCCCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTTCATGTGCCGCCACAAGCAGCGTCATGACGCTTTAATAAGCAAATTGTTgg                                                                                                                                                                                      >  2:161095/1‑251 (MQ=255)                                                 tATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCTACGCGGTCGCGGTCAGGACAGCGAAGAGGTCATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTgg                                                                                                                                                                                      >  2:69768/1‑251 (MQ=255)                                                                 aaaaTTGAACTGGACCGCCGTCTACGCGGTCGCGGTCAGGACAGCGAAGAGGTCATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATGACGCCGAATATGAGTATCTGATTGTGAATGATGACTTCGATTCCGCGTTGAACGAGTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATTACGCTTTAATCAGCAAATTGTTAGCAGTCTGAACCTGAtt                                                                                                                                                                      <  2:201819/251‑1 (MQ=255)                                                                    aTTGAACTGGACCGCCGTCTACGCGGTCGCGGTCAGGACAGCGAAGAGGTCATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTc                                                                                                                                                                    >  1:169677/1‑250 (MQ=255)                                                                        aaCTAGACCGCAGTCTACGCGGTCGCGGCAAGGACAGCGAAGAGGGCATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTAt                                                                                                                                                               >  2:605/1‑251 (MQ=255)                                                                            ggACCGCCGTCTACGCGGTCGCGGTCAGGACAGCGAAGAGGTCATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATg                                                                                                                                                           >  1:72890/1‑251 (MQ=255)                                                                                    gTCTACGCGGTCGCGGTCAGGACAGCGAAGAGGTCATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACGTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCa                                                                                                                                                   <  2:334386/251‑1 (MQ=255)                                                                                      cTACGCGGTCGCGGTCAGGACAGCGAAGAGGTCATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCAtt                                                                                                                                                 >  1:400692/1‑251 (MQ=255)                                                                                               cgcggtcAGGACAGCGAAGAGGTCATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCGGATTTCAGTATCATGCCCAGTCATTTCTTCAcc                                                                                                                                         <  1:344609/250‑1 (MQ=255)                                                                                                  ggtcAGGACAGCGAAGAGGTCATCGAAAAGCGTCTGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATAACTTCGAGACCGCGTTGACCGATTTGAAGACCATTATTAGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGAGTCATGACGCTTTAATCAGCACATTGTTTGCAGACTGAACCTGATTTCAGTATCATGCCCAGCAATTTCTTCACCtgtg                                                                                                                                     >  2:34237/1‑251 (MQ=255)                                                                                                   gtcAGGAAAGCGCAGAGGTCACTGCAAAGCGTTTGAAGCAAGCTGTTGCAGAAATGAGCCATTACGCAGAATATGATTATCTGCTTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGACCCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTgg                                                                                                                                    >  2:244848/1‑251 (MQ=255)                                                                                                        gACAGCGAAGAGGTCATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCtt                                                                                                                               <  1:134866/251‑1 (MQ=255)                                                                                                         aCAGCGAAGAGGTCATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCttt                                                                                                                              <  1:278673/251‑1 (MQ=255)                                                                                                          cAGCGAAGAGGTCATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGGTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGGCAAAAGCAGCGTCATGACGCTTTAAGCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCtttt                                                                                                                             <  1:208643/251‑1 (MQ=255)                                                                                                                      cATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCAc                                                                                                                 >  2:308793/1‑251 (MQ=255)                                                                                                                       aTTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCAcg                                                                                                                >  1:114742/1‑251 (MQ=255)                                                                                                                          gCAAAGCGTATGGCGCAAGCGGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGt                                                                                                             <  1:320718/251‑1 (MQ=255)                                                                                                                                 gTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGtt                                                                                                      >  2:81425/1‑251 (MQ=255)                                                                                                                                   agggCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGATCGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCa                                                                                                    <  2:224692/249‑1 (MQ=255)                                                                                                                                      gcgcAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAgg                                                                                                  >  1:43455/1‑250 (MQ=255)                                                                                                                                        gcAAGCTGTTGCAGACATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGTACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACg                                                                                               >  1:223082/1‑251 (MQ=255)                                                                                                                                        gcAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACg                                                                                               >  1:187790/1‑251 (MQ=255)                                                                                                                                         cAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGc                                                                                              >  2:2716/1‑251 (MQ=255)                                                                                                                                             cTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGt                                                                                           >  1:181090/1‑250 (MQ=255)                                                                                                                                              gattGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCATACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTTTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTag                                                                                         <  2:256660/249‑1 (MQ=255)                                                                                                                                                ttGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTag                                                                                         >  2:105032/1‑249 (MQ=255)                                                                                                                                                 tGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTagag                                                                                       <  1:192787/250‑1 (MQ=255)                                                                                                                                                    aGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGaaa                                                                                    >  2:208091/1‑250 (MQ=255)                                                                                                                                                          gAGCCATTACGCCGAATATGAGTATCTGATTGTCAATGATGACTGCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGGCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACGGAACCTGATTTCAGTATCATGCCCAGTCAGTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACTCGGAACTGTTTAGGTCGCTGTAGAGAAAAtttgg                                                                              <  2:150157/250‑4 (MQ=255)                                                                                                                                                                ttACGCCGAATATGAGTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGCGAAAATTGGTAACCGtt                                                                       >  2:143627/1‑251 (MQ=255)                                                                                                                                                                 tACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGttt                                                                      >  1:319386/1‑251 (MQ=255)                                                                                                                                                                  aCGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGTCACGCGTAACTGGTCAGGACGCTGTAGAGAAAATTGGTAACCGtgtt                                                                     <  2:266827/251‑4 (MQ=255)                                                                                                                                                                     ccGAATATGATTATCTGATTGTGAATGCTGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGCCTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGAc                                                                  >  1:42755/1‑251 (MQ=255)                                                                                                                                                                         atatGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGCTTTGAAGACCATTATTCGCGACGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAACTTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTgg                                                              >  1:243017/1‑251 (MQ=255)                                                                                                                                                                            tGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTAc                                                           >  2:335175/1‑251 (MQ=255)                                                                                                                                                                             gATTATCTGATTGTGAATGATGCCTTCGATCCCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTATTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCCGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACt                                                          >  2:195781/1‑251 (MQ=255)                                                                                                                                                                               ttATCTGATTCTGAATGATGACTTCTATACCGCGTTGACCGATTTGTAGACCATTATTCGCGCCGGACGTCTGTGCATGAGCCGCCAAAAGCCGCGTCATGACGCTTTAATCCGCAAATTGTTGGCAGACTGACCCTGATTTCAGTATCATGCACAGTCATTTCTTCACCTGTGGAGCTTTTTATGTATGTCACGCGTAGCTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTgg                                                        >  2:331891/1‑251 (MQ=255)                                                                                                                                                                                        tgtgAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCcgcgc                                               >  2:147478/1‑251 (MQ=255)                                                                                                                                                                                        tgtgAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGAACTGGTACTGGTCGCcgcg                                                <  1:17484/250‑1 (MQ=255)                                                                                                                                                                                           gAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCTAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACATGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGTc                                            >  1:276273/1‑251 (MQ=255)                                                                                                                                                                                           gAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGt                                             <  1:308793/250‑1 (MQ=255)                                                                                                                                                                                              tgatgaCTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCAc                                                                                                                                          >  2:313082/1‑154 (MQ=255)                                                                                                                                                                                              tgatgaCTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCAc                                                                                                                                          <  1:313082/154‑1 (MQ=255)                                                                                                                                                                                                atgaCTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGTCGCGCt                                       >  2:363695/1‑251 (MQ=255)                                                                                                                                                                                                  gaCTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGTCGCGCTCg                                     >  2:13183/1‑251 (MQ=255)                                                                                                                                                                                                      tCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGTCGCGCTCGTc                                   <  1:21998/249‑1 (MQ=255)                                                                                                                                                                                                      tCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACGGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGGGTCGCGGTCGTCa                                  <  2:400692/250‑1 (MQ=255)                                                                                                                                                                                                        gATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGTCGCGCTCGTc                                   >  2:26900/1‑247 (MQ=255)                                                                                                                                                                                                                       gATTTGCAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGTCGCGCTCGTCAGATGCAGGTAGGCGGaa                 >  2:314980/1‑250 (MQ=255)                                                                                                                                                                                                                       gATTTGAAGACCATTATTCGCGCCGAACGGCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAAGCAGCAAAGTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTc                                                                                                                                                    <  1:390207/119‑1 (MQ=255)                                                                                                                                                                                                                       gATTTGAAGACCAGTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGTCGCGCTCGTCAGATGCAGGTAGGCGGaaa                <  2:246433/251‑1 (MQ=255)                                                                                                                                                                                                                           tGAAGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGTCGCGCTCGTCAGATGCAGGTAGGCGGAAAGGAt            >  1:282411/1‑251 (MQ=255)                                                                                                                                                                                                                            gAAGACCATTATTCGCGCCGAACGTCTGTGCCTGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGTCGCGCTCGTCAGATGCAGGTAGGCGGAAAGGATc           >  1:319904/1‑251 (MQ=255)                                                                                                                                                                                                                             aaGACCATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGTCGCGCTCGTCAGATGCAGGTAGGCGGAAAGGATcc          <  2:111961/251‑1 (MQ=255)                                                                                                                                                                                                                                  cATTATTCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGTCGCGCTCGTCAGATGCCGGTAGGCGGAAAGGCTCCGCTgg     >  1:17652/1‑251 (MQ=255)                                                                                                                                                                                                                                    ttattCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTTACCTGGTACTGGTCGCCGCGCGTCGCGCTCGTCAGATGCAGGTAGGCGGAAAGGATCCGCTGGt    >  2:10895/1‑250 (MQ=255)                                                                                                                                                                                                                                     tattCGCGCCGAACGTCTGTGCATGAGCCGCCAAAAGCAGCGTCACGACGCTTTAATCAGCAAATTGTTGACAGCCTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCCCCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGCCGCGCTCGTCCGATGCAGGTAGGCGGAAAGGATCCGCTGGTAc  >  1:59960/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                        GGGCGCGCAGCAAATTCGCCAGAAGATGCCGCACGCGCGGAGTATCTTTATTTTACCGCCGTCCAAAATTGAACTGGACCGCCGTCTACGCGGTCGCGGTCAGGACAGCGAAGAGGTCATTGCAAAGCGTATGGCGCAAGCTGTTGCAGAAATGAGCCATTACGCCGAATATGATTATCTGATTGTGAATGATGACTTCGATACCGCGTTGACCGATTTGAAGACCATTATTCGCGCCGAACGTCTGCGCATGAGCCGCCAAAAGCAGCGTCATGACGCTTTAATCAGCAAATTGTTGGCAGACTGAACCTGATTTCAGTATCATGCCCAGTCATTTCTTCACCTGTGGAGCTTTTTAAGTATGGCACGCGTAACTGTTCAGGACGCTGTAGAGAAAATTGGTAACCGTTTTGACCTGGTACTGGTCGCCGCGCGTCGCGCTCGTCAGATGCAGGTAGGCGGAAAGGATCCGCTGGTAC  >  NC_000913/3821745‑3822223

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