Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,110,741	A→C	F93L (TTT→TTG)	rfbA ←	glucose‑1‑phosphate thymidylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,110,741	0	A	C	96.2%	178.5 / ‑2.0	53	F93L (TTT→TTG)	rfbA	glucose‑1‑phosphate thymidylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (21/30);  total (22/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.36e-01

CTTTGGTTCTAACGGTTTTTCTTCCAGACTGATTGCCGTACCGTTTTTATCAAACTCAACGACACCATAGCGTTCTGGATCATTAACGTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATAAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACTGATAATCAAAATATCGCGAATACCCGCCAACATCAGTGTAGAGAGCGGGTAATAGATCATCGGTTTATCATAAATAGGTAATAGCTGTTTACTGACAGCCATAGTCACA  >  NC_000913/2110516‑2110963                                                                                                                                                                                                                                  |                                                                                                                                                                                                                               ctTTGGTTCTAACGGTTTTTCTTCCAGACTGATTGCCGTACCGTTTTTATCAAACTCAACGACACCATAGCGTTCTGGATCATTAACGTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGACTgcg                                                                                                                                                                                                                      >  2:465704/1‑236 (MQ=255) ctTTGGTTCTAACGGTTTTTCTTCCAGACTGATTGCCGTACCGTTTTTATCAAACTCAACGACACCATAGCGTTCTGGATCATTAACGTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATAAATGCCTgcg                                                                                                                                                                                                                      >  2:36490/1‑236 (MQ=255)    tGGTTCTAACGGTTTTTCTTCCAGACTGATTGCCGTACCGTTTTTATCAAACTCAACGACACCATAGCGTTCTGGATCATTAACGTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGcgcg                                                                                                                                                                                                                    >  1:694020/1‑235 (MQ=255)           aaCGGTTTTTCTTCCAGACTGATTGCCGTACCGTTTTTATCAAACTCAACGACACCATAGCGTTCTGGATCATTAACGTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCAt                                                                                                                                                                                                             <  1:423416/235‑1 (MQ=255)               gTTTTTCTTCCAGACTGATTGCCTTACCTTTTTTAGCAAACTCAACGACACCATAGCTTTCTGGATCATTAACGTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTggg                                                                                                                                                                                                        <  2:383672/236‑1 (MQ=255)                      ttcCAGACTGATTGCCTTACCGTTTTTATCAAACTCAACGACACCATAGCGTTCTGGATCATTAACGTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGtt                                                                                                                                                                                                 <  1:401609/236‑1 (MQ=255)                         catgcTGATTGCCGTACCGGTTTTAGCAATCTCAACGACCCCATAGCGTTCTGGATCATTAACGTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCa                                                                                                                                                                                              <  2:270927/232‑1 (MQ=255)                          aGACTGATTGCCGTACCGTTTTGATCAAACTCAACGACTCCATAGCGTTCTGGATCATTAACGTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCCGCGCAAGGtctggc                                                                                                                                                                                             <  2:636627/236‑6 (MQ=255)                                       tACCGTTTTTATCAAACTCAACGACACCATAGCGTTCTGGATCATTAACGTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGa                                                                                                                                                                                >  1:613459/1‑236 (MQ=255)                                         ccGTTTTTATCAAACTCAACGACACCATAGCGTTCTGGATCATTAACGTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGAtt                                                                                                                                                                              <  2:486391/236‑1 (MQ=255)                                           gTTTTTATCAAACTCAACGACACCATAGCGTTCTGGATCATTAACGTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATAAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCa                                                                                                                                                                            <  1:36490/236‑1 (MQ=255)                                             ttttATCAAACTCAACGACACCATAGCGTTCTGGATCATTAACGTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAgg                                                                                                                                                                          <  1:715047/236‑1 (MQ=255)                                                  tCAAACTCAACGACACCATAGCGTTCTGGATCATTAACGTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCa                                                                                                                                                                     <  2:733835/236‑1 (MQ=255)                                                  tCAAACTCAACGACACCATAGCGTTCTGGATCATTAACGTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATACCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCa                                                                                                                                                                     <  2:694020/236‑1 (MQ=255)                                                          aaCGACACCATAGCGTTCTGGATCATTAACGTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTg                                                                                                                                                                  >  2:425263/1‑231 (MQ=255)                                                               cacCATAGCGTTCTGGATCATTAACGTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCa                                                                                                                                                        <  1:345495/236‑1 (MQ=255)                                                                    tAGCGTTCTGGATCATTAACGTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTg                                                                                                                                                                  <  1:582327/221‑1 (MQ=255)                                                                    tAGCGTTCTGGATCATTAACGTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTg                                                                                                                                                                  >  2:582327/1‑221 (MQ=255)                                                                                   ttAACGTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACg                                                                                                                                    >  2:673295/1‑236 (MQ=255)                                                                                   ttAACGTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCAACTGGCTACCGTCACCCAGCAATTGTTTAAAAct                                                                                                                                    >  2:647743/1‑235 (MQ=255)                                                                                        gTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACgaggag                                                                                                                               <  1:4856/236‑1 (MQ=255)                                                                                           aTAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGAt                                                                                                                                                                               <  1:53077/185‑1 (MQ=255)                                                                                           aTAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGAt                                                                                                                                                                               >  2:53077/1‑185 (MQ=255)                                                                                           aTAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACgaggag                                                                                                                               >  1:567822/1‑233 (MQ=255)                                                                                                   aTACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTg                                                                                                                        <  1:425263/232‑1 (MQ=255)                                                                                                         ttGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACTg                                                                                                              >  1:676759/1‑236 (MQ=255)                                                                                                             tccacTTTCATGGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACTg                                                                                                              <  2:676759/231‑1 (MQ=255)                                                                                                                 ctttctttGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACTGATAATCaa                                                                                                      >  1:592172/1‑236 (MQ=255)                                                                                                                              cagcagACTCCATTAACTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACATTCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACTGATAATCAAAATATCGCGAAt                                                                                          <  2:613459/235‑1 (MQ=255)                                                                                                                               ggCGGCCTACATTAGCTGGGGCAGATCGTGACCGTAAATGATATAATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAAGGCCTGCGCGAGGCCATCTGGGCTAGGTTGCCCTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCGGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACGGATAATCAAAATATCGCGAATAc                                                                                        <  2:163456/235‑1 (MQ=255)                                                                                                                                gCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGGGTACTGATAATCAAAATATCTCGAATaa                                                                                        >  2:265377/1‑234 (MQ=255)                                                                                                                                    ccTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACTGATAATCAAAATATCGCGAATACCCGcc                                                                                   >  2:419409/1‑236 (MQ=255)                                                                                                                                    ccTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTAGCCTTAGGTGTACTGATAATCAAAATATCGCGAATAccacac                                                                                   >  2:30176/1‑232 (MQ=255)                                                                                                                                      tCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACTGATAATCAAAATATCGCGAATACCCGCCaa                                                                                 <  2:592172/236‑1 (MQ=255)                                                                                                                                       ccATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACTGATAATCAAAATATCGCGAATACCCGCCaa                                                                                 <  1:465704/235‑1 (MQ=255)                                                                                                                                                    aaGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACTGATAATCAAAATATCGCGAATACCCGCCAACATCAGTGTagaga                                                                   <  1:647743/235‑1 (MQ=255)                                                                                                                                                     aGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACTGATAATCAAAATATCGCGAATACCCGCCAACATCAGTGTagagag                                                                  <  1:419409/236‑1 (MQ=255)                                                                                                                                                     aGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACTGATAATCAAAATATCGCGAATACCCGCCAACATCAGTGTagagag                                                                  <  1:673295/236‑1 (MQ=255)                                                                                                                                                     aGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGCCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACTGATAATCAAAATATCGCGAATACCCGCCAACATCAGTGTagagag                                                                  <  1:685180/236‑1 (MQ=255)                                                                                                                                                         cGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGcccc                                                                                                                                                                      >  2:230581/1‑132 (MQ=255)                                                                                                                                                         cGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGcccc                                                                                                                                                                      <  1:230581/132‑1 (MQ=255)                                                                                                                                                                        tattatCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACTGATAATCAAAATATCGCGAATACCCGCCAACATCAGTGTAGAGAGCGGGTAATAGATCATCGGt                                               <  2:486885/236‑1 (MQ=255)                                                                                                                                                                         attatCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACTGATAATCAAAATATCGCGAATACCCGCCAACATCAGTGTAGAGAGCGGGTAATAGATCATCGGt                                               >  1:173352/1‑235 (MQ=255)                                                                                                                                                                                  aaGAACCAAAGCAAAATCATCACCACCAATAATCCCTTCACCGATGATCAATGCATGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACTGATAATCAAAATATCGCGAATACCCGCCAACATCAGTGTAGAGAGCGGGTAATAGATCATCGGTTTATCATaaa                                     <  2:444096/236‑1 (MQ=255)                                                                                                                                                                                               caATCATCACCACCAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACTGATAATCAAAATATCGCGAATACCCGCCAACATCAGTGTAGAGAGCGGGTAATAGATCATCGGTTTATCATAAATAGGTa                               >  1:160829/1‑229 (MQ=255)                                                                                                                                                                                                       accaccAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTg                                                                                                                                                                                    >  1:217998/1‑72 (MQ=255)                                                                                                                                                                                                       accaccAATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTg                                                                                                                                                                                    <  2:217998/72‑1 (MQ=255)                                                                                                                                                                                                            caATAAACTCTTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACTGATAATCAAAATATCGCGAATACCCGCCAACATCAGTGTAGAGAGCGGGTAATAGATCATCGGTTTATCATAAATAGGTAATAGCTGTTTACTGACAGcc           >  1:419886/1‑236 (MQ=255)                                                                                                                                                                                                                   ctctTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGATAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACTGATAATCAAAATATCGCGAATACCCGCCAACATCAGTGTAGAGAGCGGGTAATAGATCATCGGTTTATCATAAATAGGTAATAGCTGTTTACTGACAGCCATAGTca    <  1:265377/236‑1 (MQ=255)                                                                                                                                                                                                                   ctctGCACGGATGATCAATCCCTGCACGAGTCCATCGGGTCTAGGTTGCTATTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTAAAAACGAGGAGTATCCTGAGGTGTACTGATAATCAAAATATCGCGAATCCCCGCCAACATCAGTGTAGAGAGCGGGTAATAGATCATCGGTTTATCATAAATAGGTAATAGCTGTTTACTGACAGCCATAGTca    <  2:173352/236‑1 (MQ=255)                                                                                                                                                                                                                     ctTCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACTGATAATCAAAATATCGCGAATACCCGCCAACATCAGTGTAGAGAGCGGGTAATAGATCATCGGTTTATCATAAATAGGTAATAGCTGTTTACTGACAGCCATAGTcaca  >  1:599019/1‑236 (MQ=255)                                                                                                                                                                                                                       tCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACTGATAATCAAAATATCGCGAATACCCGCCAACa                                                                               >  1:663959/1‑157 (MQ=255)                                                                                                                                                                                                                       tCACCGATGATCAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACTGATAATCAAAATATCGCGAATACCCGCCAACa                                                                               <  2:663959/157‑1 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                               CTTTGGTTCTAACGGTTTTTCTTCCAGACTGATTGCCGTACCGTTTTTATCAAACTCAACGACACCATAGCGTTCTGGATCATTAACGTGATAGGCAAATACCGTTGCACCACTTTCTTTGTTAACAGCGGCCTCCATTAGCTTCGGCAGATCGTGACCGTAAAAGATATTATCACCAAGAACCAAAGCACAATCATCACCACCAATAAACTCTTCACCGATGATAAATGCCTGCGCGAGGCCATCTGGGCTAGGTTGCACTTTGTACTGAAGATTCAGGCCCCACTGGCTACCGTCACCCAGCAATTGTTGAAAACGAGGAGTATCCTGAGGTGTACTGATAATCAAAATATCGCGAATACCCGCCAACATCAGTGTAGAGAGCGGGTAATAGATCATCGGTTTATCATAAATAGGTAATAGCTGTTTACTGACAGCCATAGTCACA  >  NC_000913/2110516‑2110963

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