Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,195,757	A→T	6.3%	S115T (TCC→ACC)	thiC ←	phosphomethylpyrimidine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,195,757	0	A	T	6.3%	170.9 / 2.3	65	S115T (TCC→ACC)	thiC	phosphomethylpyrimidine synthase
Reads supporting (aligned to +/- strand):  ref base A (33/27);  new base T (2/2);  total (36/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

TTCCATTTCCGGCGTGATGATGCCCTGGCGGGCGTAGTGCAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAGGTGTCGTAGACCGGAATCGCTTCGTTTTCTTCG  >  NC_000913/4195611‑4195905                                                                                                                                                   |                                                                                                                                                     ttCCATTTCCGGCGTGATGATGCCCTGGCGGGCGTAGTGCAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAAc                                                                                                                                                    >  1:212621/1‑149 (MQ=255)     aTTTCCGGCGTGATGATGCCCTGGCGGGCGTAGTGCAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTTGTGTAATCGGAACTgcg                                                                                                                                                >  4:114816/1‑149 (MQ=255)          cGGCGTGATGATGCCCTGGCGGGCGTAGTGCAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACgg                                                                                                                                           <  5:232085/149‑1 (MQ=255)           ggCGTGATGATGCCCTGGCGGGCGTAGTGCAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGTACTGCGCACgg                                                                                                                                           >  8:100542/1‑148 (MQ=255)             cGTGATGATGCCCTGGCGGGCGTAGTGCAGTTGGGTGACGCGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTaa                                                                                                                                        <  5:283013/149‑1 (MQ=255)               tgatgaTGCCCTGGCGGGCGTAGTGCAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGt                                                                                                                                      >  6:13222/1‑149 (MQ=255)                 atgatgCCCTGGCGGGCGTAGTGCAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGtt                                                                                                                                     <  5:242702/148‑1 (MQ=255)                    atgCCCTGGCGGGCGTAGTGCAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGtgcttc                                                                                                                                 <  2:232823/149‑1 (MQ=255)                        ccTGGCGGGCGTAGTGCAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTa                                                                                                                             >  2:49794/1‑149 (MQ=255)                                gcgTAGTGCAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGaa                                                                                                                                                     <  1:375430/117‑1 (MQ=255)                                gcgTAGTGCAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGaa                                                                                                                                                     >  2:375430/1‑117 (MQ=255)                                     gTGCAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCa                                                                                                                 <  7:64779/148‑1 (MQ=255)                                      tGCAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATc                                                                                                               >  5:349141/1‑149 (MQ=255)                                        cAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGa                                                                                                             <  8:98753/149‑1 (MQ=255)                                               gTGACACGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACgg                                                                                                      >  1:335605/1‑149 (MQ=255)                                                            cTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGTACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGtttt                                                                                         <  7:100542/149‑1 (MQ=255)                                                            cTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGtttt                                                                                         <  2:216515/149‑1 (MQ=255)                                                            cTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGtttt                                                                                         <  2:299411/149‑1 (MQ=255)                                                            cTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATACACGGCTGGCGTAGtttt                                                                                         >  4:101006/1‑149 (MQ=255)                                                            cTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTACTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGtttt                                                                                         >  4:261392/1‑149 (MQ=255)                                                                     cgcgTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCt                                                                                >  4:27450/1‑149 (MQ=255)                                                                      gcgTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCtt                                                                               >  5:76486/1‑149 (MQ=255)                                                                       cgTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTtg                                                                              <  5:13222/149‑1 (MQ=255)                                                                       cgTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTtg                                                                              <  3:261392/149‑1 (MQ=255)                                                                        gTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGTCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTtgc                                                                             >  4:345419/1‑149 (MQ=255)                                                                        gTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTtgc                                                                             <  6:76486/149‑1 (MQ=255)                                                                          tttGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGGACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTtgctg                                                                           >  6:9435/1‑149 (MQ=255)                                                                            cggTGTTAGTACGCCGCTAAAACGCAGTTCTTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTa                                                                                                                             <  8:187997/96‑1 (MQ=255)                                                                                 ttAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGtt                                                                     >  5:364690/1‑148 (MQ=255)                                                                                 ttAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGtt                                                                     >  7:160398/1‑148 (MQ=255)                                                                                    gTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAAt                                                                  <  3:122151/148‑1 (MQ=255)                                                                                      aCGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTAGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATgg                                                                >  1:19571/1‑148 (MQ=255)                                                                                       cgccgcTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCa                                                              <  8:136135/149‑1 (MQ=255)                                                                                        gccgcTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCcgcgc                                                                                                                    >  1:286264/1‑94 (MQ=255)                                                                                        gccgcTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCcgcgc                                                                                                                    <  2:286264/94‑1 (MQ=255)                                                                                          cgcTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATc                                                           >  6:180844/1‑149 (MQ=255)                                                                                           gcTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCTCTTTAGTGTAATCGGAACTGCGCAAGGTAAGTTCTTCGGTATAGCCTCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTTCAAGTTAATGGCAATCt                                                          >  2:365745/1‑149 (MQ=255)                                                                                             tAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGc                                                        >  3:91348/1‑149 (MQ=255)                                                                                             tAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGc                                                        >  1:151071/1‑149 (MQ=255)                                                                                               aaaCGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCg                                                       <  2:335605/148‑1 (MQ=255)                                                                                               aaaCGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATAGGAACTGCGCACGGTAAGATCTTCGCTATCGCCGCGCGCATCGATCCACGTCTGGCTTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTCCg                                                       >  1:204870/1‑148 (MQ=255)                                                                                                 aCGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCCCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGAt                                                    >  8:245769/1‑149 (MQ=255)                                                                                                         cgtcgAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATAc                                            >  6:85179/1‑149 (MQ=255)                                                                                                           tcgAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGTACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACgg                                          >  1:288121/1‑149 (MQ=255)                                                                                                                 ccATCATCTGCCAGCCGCGCTTTAGTGTAATCGGTACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGa                                    <  2:288121/149‑1 (MQ=255)                                                                                                                 ccATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGa                                    <  1:49794/149‑1 (MQ=255)                                                                                                                  catcatCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAg                                   >  1:20528/1‑149 (MQ=255)                                                                                                                  catcatCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAg                                   >  4:46613/1‑149 (MQ=255)                                                                                                                   atcatcTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATa                                             <  2:37344/138‑1 (MQ=255)                                                                                                                   atcatcTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATa                                             >  1:37344/1‑138 (MQ=255)                                                                                                                    tcatcTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAGgt                                 >  1:360151/1‑149 (MQ=255)                                                                                                                     catcTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAGgtg                                >  3:213848/1‑149 (MQ=255)                                                                                                                      atcTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAGgtgt                               <  6:349141/149‑1 (MQ=255)                                                                                                                      atcTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAGgtg                                <  3:46613/148‑1 (MQ=255)                                                                                                                       tcagCCAGCCGCGCTTTAGTGTAATCTGAACTGCGGAAGGTATTTTGTTGGGTATCGCCGCGAGCATCGATGGAGGGGTTTCGTAGTTTTGCTAGCGCTTGCTGCTGGTTATTGGCATTCTGCGGATTACCATACGGGCGGGAGGTGTc                              <  2:204870/146‑1 (MQ=255)                                                                                                                          gCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGt                                                                      >  6:239566/1‑106 (MQ=255)                                                                                                                          gCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGt                                                                      <  5:239566/106‑1 (MQ=255)                                                                                                                            cAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAGGTGTCGTAGa                         <  4:91348/149‑1 (MQ=255)                                                                                                                               ccGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGc                                                                                       >  7:49372/1‑84 (MQ=255)                                                                                                                               ccGCGCTTTAGTGTAATAGGATCTGCGCACTGTTAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGGTGGGGTAGTTTTGc                                                                                       <  8:49372/84‑1 (MQ=255)                                                                                                                                         gtgtAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGGTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGTAGGTGTCGTAGACCGGAATCGCTTc            <  5:85179/149‑1 (MQ=255)                                                                                                                                          tgtAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAGGTGTCGTa                           >  6:174147/1‑133 (MQ=255)                                                                                                                                          tgtAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAGGTGTCGTa                           <  5:174147/133‑1 (MQ=255)                                                                                                                                          tgtAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAGGTGTCGTAGACCGGAATCGCTTCg           <  2:360151/149‑1 (MQ=255)                                                                                                                                            tAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAGGTGTCGTAGACCGGAATCGCTTCGtt         >  7:278340/1‑149 (MQ=255)                                                                                                                                            tAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAGGTGTCGTAGACCGGAATCGCTTCGtt         >  5:251700/1‑149 (MQ=255)                                                                                                                                               tCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAGGTGTCGTAGACCGGAATCGCTTCGTTttc      >  6:219071/1‑149 (MQ=255)                                                                                                                                                   aaCTGCGCACGGTAAGGTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAGGTGTCGTAGACCGGAATCGCTTCGTTTTCTTCg  <  6:364524/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     TTCCATTTCCGGCGTGATGATGCCCTGGCGGGCGTAGTGCAGTTGGGTGACACGGCGTCCTGCTTTGGCGCGTTTTGGTGTTAGTACGCCGCTAAAACGCAGTTCGTCGAGGCCATCATCTGCCAGCCGCGCTTTAGTGTAATCGGAACTGCGCACGGTAAGTTCTTCGGTATCGCCGCGCGCATCGATCCACGGCTGGCGTAGTTTTGCCAGCCCTTGCTGCACGTTAATGGCAATCTGCGGATCACCATACGGGCCGGAGGTGTCGTAGACCGGAATCGCTTCGTTTTCTTCG  >  NC_000913/4195611‑4195905

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