Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	670,635	A→T	L324Q (CTG→CAG)	holA ←	DNA polymerase III subunit delta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	670,635	0	A	T	96.4%	202.0 / NA	56	L324Q (CTG→CAG)	holA	DNA polymerase III subunit delta
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (29/25);  minor base G (0/2);  total (29/27)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

TTCAAGCATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTT  >  NC_000913/670390‑670881                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       ttCAAGCATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGctct                                                                                                                                                                                                                                                   >  2:258924/1‑251 (MQ=255)                 cGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACtgactgac                                                                                                                                                                                                                                   >  2:327877/1‑251 (MQ=255)                  gCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGAcc                                                                                                                                                                                                                                  >  2:256439/1‑251 (MQ=255)                            gTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGt                                                                                                                                                                                                                        <  1:90391/251‑1 (MQ=255)                                gCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTg                                                                                                                                                                                                                     >  1:179004/1‑250 (MQ=255)                                              aCGATGCGGAGGAACGTTATGCGGGAGGATTGTGACCAGCGTCAGACCAATCGAGTTCGCCAGCGTTTCCACTGGTTTTAGATGACCGTAGTGCACCGGATCAAAGGTGCCGACAAACAGAGCCTGTAGAGATGTCATATCAACCGTCGATAAATACGTCCGCCAGGGGGTTATGGCAGAACATAAGAGATAACCCTTCCGGCTCTGCCCATACTGACTGAACGTAATCTTGGTTGAGTGTGAGttacgct                                                                                                                                                                                                      <  2:456218/251‑6 (MQ=255)                                                gATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCg                                                                                                                                                                                                    >  1:211850/1‑251 (MQ=255)                                                             attattAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGc                                                                                                                                                                                       <  1:327877/251‑1 (MQ=255)                                                                       atgatTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGTCCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCGGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGAc                                                                                                                                                                             <  2:282582/251‑1 (MQ=255)                                                                             gtgACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGAGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACt                                                                                                                                                                       >  2:71208/1‑251 (MQ=255)                                                                                        cAGACCAATCAAAGTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTTCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGACGATAAATGCGTCCGCCAGGGGTTTATGGCACAACAGAATAGATAACCCTTCCTGCTCTGCCCACACGGACTGACCGTAATCGTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTGGCACGGCCTGACGTAACTGCGTCTGACtt                                                                                                                                                            <  2:130733/251‑1 (MQ=255)                                                                                         aGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTAGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACtt                                                                                                                                                            >  1:196645/1‑250 (MQ=255)                                                                                              aaTCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCg                                                                                                                                                       >  2:79283/1‑250 (MQ=255)                                                                                               aTCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGAt                                                                                                                                                     >  2:211053/1‑251 (MQ=255)                                                                                               aTCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGAt                                                                                                                                                     <  1:256439/251‑1 (MQ=255)                                                                                                 cAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGAttt                                                                                                                                                   <  1:65157/251‑1 (MQ=255)                                                                                                   aaTTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTaa                                                                                                                                                 <  1:13925/251‑1 (MQ=255)                                                                                                    aTTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAAc                                                                                                                                                <  2:510371/251‑1 (MQ=255)                                                                                                        gCCAGCGTTTCCACGGGTTTTAGATGCCCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGCggg                                                                                                                                                                                                                 >  2:223462/1‑182 (MQ=255)                                                                                                        gCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAggg                                                                                                                                                                                                                 <  1:223462/182‑1 (MQ=255)                                                                                                        gCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCt                                                                                                                                            >  1:338610/1‑251 (MQ=255)                                                                                                         ccAGGGTTTCCTCGTGTTGTAGATGACCAAAGTGCACCGGATAAAAGGGGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCCGGGTTTTATGGCATAACAGAAGAGATAACCCTGCCGGCTCTGCCCACACTGACGGACCGTAAGCTTGTTTGAGGGTTAGTTCCGTTCGTGTCAGGAGTTGTACTGCCTGTCGTGACTGCGTTTTTCTTAAGCGTTTCAACGCCTc                                                                                                                                           <  2:436267/251‑1 (MQ=255)                                                                                                              gTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGcccc                                                                                                                                      >  1:60526/1‑250 (MQ=255)                                                                                                              gTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCca                                                                                                                                      <  2:36860/251‑1 (MQ=255)                                                                                                                 tCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCcatc                                                                                                                                    <  2:77595/250‑1 (MQ=255)                                                                                                                  ccACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCcatcat                                                                                                                                  >  1:118889/1‑251 (MQ=255)                                                                                                                    aCGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGc                                                                                                                                >  1:88294/1‑251 (MQ=255)                                                                                                                        gTTTTAGATGACCATAGAGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGTGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACAGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGCATTAAGCGATTTAACGCCTCGCCCATCATGCCCcg                                                                                                                            >  2:517474/1‑251 (MQ=255)                                                                                                                          tttAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCcgg                                                                                                                           <  1:306355/250‑1 (MQ=255)                                                                                                                            tAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCcggcg                                                                                                                         <  1:103656/250‑1 (MQ=255)                                                                                                                                          gCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGTAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATAccc                                                                                                           <  1:262249/250‑1 (MQ=255)                                                                                                                                                 aTCAAGGGTGCAGCAAAACAGAGCCTGTAAAGAGTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGGGAGATCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTACGCGATTTAACGCCTCGCCCATCGTGCCCCTGCGGTTCTGCCATAcccagtgcttt                                                                                                   >  2:240200/1‑243 (MQ=255)                                                                                                                                                     aaGGTGCCGCCAAACAGAGCCTGTAAAGCTTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCCCGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATcc                                                                                                >  2:408237/1‑249 (MQ=255)                                                                                                                                                            cgccAAACGGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGAAAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTACCTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGGTCTGCCATACCCGATGCTTAtacaaccacgc                                                                                        >  2:340185/1‑241 (MQ=255)                                                                                                                                                                    agatccTGTAAAGATTTCATATCACCCGTCGATAAATACGTCCGCCAGCGGTTTATGGCTCAACAGAAGAGATAACCCTTCCTGCTCTGCCCACGCTGACTGACCGTAATCTTGTTTGAGGTTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCGGACTTAAGCGATTTAACGCCTCTCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACTCTCGCAGTg                                                                                <  2:524431/247‑1 (MQ=255)                                                                                                                                                                       ggcTGTAAAGATTTCATATCAACCGTCGATAACTACGTCCGCCAGGGGTTTATGGCACAACATAAGAGATAACCCGTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCg                                                                             <  2:202515/249‑1 (MQ=255)                                                                                                                                                                           gTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATg                                                                         <  1:248288/251‑1 (MQ=255)                                                                                                                                                                                    tCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACt                                                                 <  2:51618/250‑1 (MQ=255)                                                                                                                                                                                    tCATATCAACCGTCGATAAATACGTCCGCCAGGCGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGAAGCTTATCAAACAACGCAAGCAGGTGCGTATGGGCAGACTg                                                                >  2:178804/1‑251 (MQ=255)                                                                                                                                                                                           aaCCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTACCGCCTCGCCCATCATGCCCCGGCGGTTATGCCATACCCGATGCTTCTCAACCAACGCCCCCAGTGGCGGATGGGCAGACTGGCGTTTc                                                         >  2:428961/1‑251 (MQ=255)                                                                                                                                                                                            aCCGTCGATAAATTCGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACGTACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAGCAACGCACGCAGTGGGGTATGGGTAGACTGGCGTTTCa                                                        <  2:73553/251‑1 (MQ=255)                                                                                                                                                                                                 cGATAAATTCTTCCGCCAGGTGTTTATGGCACAACAGAAGAGTTAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCCCGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGTAGACTGGCGTTTCAGGTTa                                                   <  2:74840/251‑1 (MQ=255)                                                                                                                                                                                                 cGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGtt                                                    <  1:211053/250‑1 (MQ=255)                                                                                                                                                                                                              cGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTaacaac                                       >  1:113980/1‑250 (MQ=255)                                                                                                                                                                                                                 cAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGtt                                   <  2:118889/251‑1 (MQ=255)                                                                                                                                                                                                                    gggTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCAc                                >  1:452982/1‑251 (MQ=255)                                                                                                                                                                                                                    gggTTTATGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCAc                                >  2:94751/1‑251 (MQ=255)                                                                                                                                                                                                                          aTGGCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGt                           >  1:7135/1‑250 (MQ=255)                                                                                                                                                                                                                            ggCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAAt                        <  2:88294/251‑1 (MQ=255)                                                                                                                                                                                                                             gCACAACAGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAAtg                       <  1:258924/251‑1 (MQ=255)                                                                                                                                                                                                                                  acaGAAGAGATAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAAAAGTTCACGTTGTAATGTgcgc                  >  1:28072/1‑251 (MQ=255)                                                                                                                                                                                                                                          gaTAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAAt          >  1:500566/1‑251 (MQ=255)                                                                                                                                                                                                                                            tAACCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATa         >  1:82734/1‑250 (MQ=255)                                                                                                                                                                                                                                             aaCCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAAc       >  1:339249/1‑251 (MQ=255)                                                                                                                                                                                                                                             aaCCCTTCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAAc       <  1:329336/251‑1 (MQ=255)                                                                                                                                                                                                                                                  ttCCTGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGGGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGGTCACGTTGTAATGTGCGCAACCAAATAACCGGtt  >  2:157275/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       TTCAAGCATGTGTTTACGCTGCACGCTGTTCGCTTCCGGCTGGGGACGATGCGGAGGAACATTATTAGGGATGATTGTGACCCGCGTCAGACCAATCAAATTCGCCAGCGTTTCCACGGGTTTTAGATGACCATAGTGCACCGGATCAAAGGTGCCGCCAAACAGAGCCTGTAAAGATTTCATATCAACCGTCGATAAATACGTCCGCCAGGGGTTTATGGCACAACAGAAGAGATAACCCTTCCAGCTCTGCCCACACTGACTGACCGTAATCTTGTTTGAGGGTGAGTTCCGTTCGTGTCAGGAGTTGCACGGCCTGACGTAACTGCGTCTGACTTAAGCGATTTAACGCCTCGCCCATCATGCCCCGGCGGTTCTGCCATACCCGATGCTTATCAAACAACGCACGCAGTGGCGTATGGGCAGACTGGCGTTTCAGGTTAACCAGTAACAACAGTTCACGTTGTAATGTGCGCAACAAAATAACCGGTT  >  NC_000913/670390‑670881

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