Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	557,635	A→T	L36Q (CTG→CAG)	folD ←	bifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	557,635	0	A	T	97.8%	150.0 / NA	46	L36Q (CTG→CAG)	folD	bifunctional methylenetetrahydrofolate dehydrogenase/methenyltetrahydrofolate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (17/28);  minor base C (1/0);  total (18/28)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACC  >  NC_000913/557391‑557870                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            ggcggaataaaaacttaccacgaCTTTGACTTTACCACTACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTg                                                                                                                                                                                                                                       <  2:93772/229‑1 (MQ=255)    cGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTg                                                                                                                                                                                                                                    >  2:217012/1‑251 (MQ=255)       aTGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGACGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCACCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGCGCGGGAGACGAACACGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTCCTAACCACCAG‑ACCACAGGCCTGTCCTGGTGccc                                                                                                                                                                                                                                 >  1:188566/1‑251 (MQ=255)         gAATACGTTCAAGAACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGCCTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGATGCCCGc                                                                                                                                                                                                                               >  1:72633/1‑251 (MQ=255)            tACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGAAGTTGAAACAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTTTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCCGTCCTGGTGCCCGCatg                                                                                                                                                                                                                            >  2:274056/1‑249 (MQ=255)             aCGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGCTAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTc                                                                                                                                                                                                                           >  1:408409/1‑251 (MQ=255)             aCGTTCAAGCACTTTGCAGTTATCAATACCCGCCGGTAACGGCAGCTGACACAGAATGCCCTCGCTGGCGATGCCGGCATTCAGCGTATAGATAAGATCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGCCGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGCTGCCGGGTTACTACCCCCCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTc                                                                                                                                                                                                                           >  2:404764/1‑251 (MQ=255)                  ccAGCTCTTTGGCTTTATCACTCCCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGc                                                                                                                                                                                                                       <  1:84082/250‑1 (MQ=255)                       aCTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAAt                                                                                                                                                                                                                 >  1:55715/1‑251 (MQ=255)                            tacgtcttccccccccGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGGCCGGCTGCAATAcgc                                                                                                                                                                                                             <  1:85947/237‑1 (MQ=255)                               gTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCt                                                                                                                                                                                                         >  2:39369/1‑251 (MQ=255)                                ttactaacccccGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTGGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGcc                                                                                                                                                                                                          <  2:187074/240‑1 (MQ=255)                                    aaaTACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACt                                                                                                                                                                                                    <  1:153351/250‑1 (MQ=255)                                      atccccGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACttt                                                                                                                                                                                                  <  2:95128/248‑1 (MQ=255)                                       cccccGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACtttt                                                                                                                                                                                                 <  1:39369/249‑1 (MQ=255)                                                gTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACt                                                                                                                                                                                        <  2:119313/251‑1 (MQ=255)                                                  aaCGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGGTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTTCCCGCAGTCCGGCTGCAATACGCGTCTGAACGTTTTGAGCAACtt                                                                                                                                                                                       <  2:249538/250‑1 (MQ=255)                                                   aCGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCa                                                                                                                                                                                     <  2:229135/251‑1 (MQ=255)                                                                cAGAATGCCATCGAGGGGTTTGGCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGGCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTGGCTTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTGCCCACCAGCACAAC‑GGGCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACctgct                                                                                                                                                                        <  1:208458/251‑1 (MQ=255)                                                                       ccATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGACCTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAAt                                                                                                                                                                 >  1:234539/1‑251 (MQ=255)                                                                            gATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCCTAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCCCCAGCACAAC‑GGCATGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCACCTTCAGAGCGCACCTGCTGCGCAATCGttt                                                                                                                                                            >  1:165410/1‑251 (MQ=255)                                                                             aTGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGTGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGtttt                                                                                                                                                           <  1:88431/251‑1 (MQ=255)                                                                                gtgtTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTTCTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTAcc                                                                                                                                                        <  1:216319/251‑1 (MQ=255)                                                                                     gtCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCa                                                                                                                                                    >  1:245774/1‑250 (MQ=255)                                                                                         gCATTCAGCGGATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCaataat                                                                                                                                               <  2:341520/251‑1 (MQ=255)                                                                                             tCAGCGTAACGAGAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTGCCGGGAGGTCATAAGAGCGGGAGACGAACTCTTCTTCTTCACAAGCCTTTCCTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACGCCCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTTCGGGTGCAATACGCTCCTTAAGGGTGTGAACAACTTCAGATCGCACCTGCTGCGCAATCGTATTACCGTCACTAATCtt                                                                                                                                            <  2:23823/250‑1 (MQ=255)                                                                                                       gatcaGCTACAGCAGAACCGCTTCGCTGGTGGATTCCGGGAGGTCATGAGCGCGGGAGAAGCACCCGACTTCTTCGCAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACACC‑GGCCTGTCCTGGTGACCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGACCGCACCTGCTGCGCACTCGTTTTACCGTCCATAATCTTTGCTGCCAta                                                                                                                                  >  2:138484/5‑249 (MQ=255)                                                                                                               ccAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGGCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAgg                                                                                                                         <  2:9929/251‑1 (MQ=255)                                                                                                                 agcagcTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGAt                                                                                                                       >  1:176061/1‑251 (MQ=255)                                                                                                                         cGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTg                                                                                                               >  1:83992/1‑251 (MQ=255)                                                                                                                              cGCTGGTGGGTGCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTAc                                                                                                           <  2:55715/250‑1 (MQ=255)                                                                                                                                    tggtTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGCTGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGt                                                                                                         >  2:174516/1‑246 (MQ=255)                                                                                                                                                        gagCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTAtttt                                                                                <  2:43782/251‑1 (MQ=255)                                                                                                                                                              gagaCGAACCCGACTTCTTCACAAGCCTTTCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGGCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGAGGCGCCTAGTTTGTCAGa                                                                          <  2:179630/251‑1 (MQ=255)                                                                                                                                                              gagaCGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGa                                                                          <  2:83992/251‑1 (MQ=255)                                                                                                                                                                  cGAACCCGACTTCTTCACAAGCCGTGCGTTTGCTTGCGACATAAATTGGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCg                                                                      <  2:408409/251‑1 (MQ=255)                                                                                                                                                                                  acaAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCCATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGAAGATCAAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgt                                                      >  2:299813/1‑251 (MQ=255)                                                                                                                                                                                       ccTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGATGATTCCATCTGTTACGTAGATCGAAGGTGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgtttcgt                                                 <  1:351901/251‑1 (MQ=255)                                                                                                                                                                                       ccTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgtttcgt                                                 >  2:343878/1‑251 (MQ=255)                                                                                                                                                                                       ccTTGCGTTTGCTTGCGACAGAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGgtttcgt                                                 <  2:157383/251‑1 (MQ=255)                                                                                                                                                                                                 gcttgcGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGAttt                                        <  1:22185/250‑1 (MQ=255)                                                                                                                                                                                                                 tGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTGACGTAGATCTAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGAcc                        <  2:66435/250‑1 (MQ=255)                                                                                                                                                                                                                      tGCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGc                  >  1:99243/1‑251 (MQ=255)                                                                                                                                                                                                                       gCAGGGTTACTACCCACCAGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCg                 <  2:35106/251‑1 (MQ=255)                                                                                                                                                                                                                                    ccaccaGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCa    <  1:217012/251‑1 (MQ=255)                                                                                                                                                                                                                                      cccaGCACAAC‑GGCCTGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCAcc  <  1:343878/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            GTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACCCGCCGGTAACGGCAGTTGAACCAGAATGCCATCGATGGTGTTGTCGGCATTCAGCGTATCGATAAGCTCCAGCAGCTCCGCTTCGCTGGTGGTTTCCGGGAGGTCATAAGAGCGGGAGACGAACCCGACTTCTTCACAAGCCTTGCGTTTGCTTGCGACATAAATTTGCGATGCAGGGTTACTACCCACCAGCACAAC‑GGCCAGTCCTGGTGCCCGCAGTCCGGCTGCAATACGCGCCTGAACTTTTTGAGCAACTTCAGAGCGCACCTGCTGCGCAATCGTTTTACCGTCAATAATCTTTGCTGCCATCAGAGAGAGGATTCCATCTGTTACGTAGATCGAAGGGGATGCGCCTATTTTGTCAGAAGCGGGGCGCGCTGTCAGGTTTCGTTTCAGATTTATCGCGTGAAGCGACCTCTTGCGAAGGTGAGGCGCACC  >  NC_000913/557391‑557870

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