Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,620,278	C→A	42.6%	R226S (CGT→AGT)	deoB →	phosphopentomutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,620,278	0	C	A	42.6%	47.0 / 85.1	75	R226S (CGT→AGT)	deoB	phosphopentomutase
Reads supporting (aligned to +/- strand):  ref base C (23/20);  new base A (18/14);  total (41/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

GATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAA  >  NC_000913/4620133‑4620424                                                                                                                                                  |                                                                                                                                                   gATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTa                                                                                                                                                 >  6:165444/1‑149 (MQ=255)   ttGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTAcc                                                                                                                                               <  8:294758/149‑1 (MQ=255)          ccATGAAGAAACTTTCTGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTAATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTaa                                                                                                                                          <  7:260896/147‑1 (MQ=255)           cATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGCCCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTACCCg                                                                                                                                       >  4:7336/1‑149 (MQ=255)                  aaaCTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGAc                                                                                                                                <  7:34312/149‑1 (MQ=255)                    aCTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCt                                                                                                                              <  6:2872/149‑1 (MQ=255)                      tttCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTgg                                                                                                                            >  8:28858/1‑149 (MQ=255)                            tCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTg                                                                                                                      <  8:122814/149‑1 (MQ=255)                                gATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAgcc                                                                                                                  >  2:17859/1‑149 (MQ=255)                                gATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAgcc                                                                                                                  <  5:1908/149‑1 (MQ=255)                                gATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAgcc                                                                                                                  >  1:33716/1‑149 (MQ=255)                                gATAAACTCTACGAACGGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATCTCGGGCGGGTTATCGCTGTTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCGGTTGAgcc                                                                                                                  <  8:273321/149‑1 (MQ=255)                                    aaCTCTACGAACTGTGCGAAATCGCACGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTACAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCa                                                                                                              >  8:17307/1‑149 (MQ=255)                                     aCTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAg                                                                                                             >  1:153740/1‑149 (MQ=255)                                         tACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGtt                                                                                                                       >  7:78211/1‑135 (MQ=255)                                         tACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGtt                                                                                                                       <  8:78211/135‑1 (MQ=255)                                         tACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCacc                                                                                                         >  7:49915/1‑149 (MQ=255)                                           cGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCac                                                                                                          <  2:267389/146‑1 (MQ=255)                                             aaCTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGGTCGTCCGTTTATCGGCGACAAAGCCGGTATCTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCaccgacc                                                                                                     <  5:124348/149‑1 (MQ=255)                                                  tgCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGtgct                                                                                                <  2:153740/149‑1 (MQ=255)                                                    cGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGtgctgc                                                                                              <  3:18930/149‑1 (MQ=255)                                                    cGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGtgctgc                                                                                              <  7:28858/149‑1 (MQ=255)                                                    cGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGtgctg                                                                                               >  3:83234/1‑148 (MQ=255)                                                     gAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGtgctg                                                                                               >  3:225608/1‑147 (MQ=255)                                                     gAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGtg                                                                                                  <  6:386040/144‑1 (MQ=255)                                                     gAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGtg                                                                                                  >  5:386040/1‑144 (MQ=255)                                                       aaTCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGa                                                                                           <  4:259157/149‑1 (MQ=255)                                                         tCGCCCGTGAAGAGCTGACCACCGGCGGCTACAATATCGGGCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGaaa                                                                                         <  5:127802/149‑1 (MQ=255)                                                         tCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGaaa                                                                                         >  3:219207/1‑149 (MQ=255)                                                           gCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATAGCTCGTCCGTTTATCTGCTACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACt                                                                                       >  6:122879/1‑149 (MQ=255)                                                             ccGTGAAGAGCTGACCAACGGCGTCTACAATATCGGTCTTGTTATCGCTCGTCCGTTTATCGGAGACAAAGCCGGAAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTgg                                                                                     >  6:70540/1‑149 (MQ=255)                                                               gTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGtt                                                                                   >  1:346250/1‑149 (MQ=255)                                                                    gagCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTtgatg                                                                               >  6:237340/1‑148 (MQ=255)                                                                          aCCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAaca                                                                        >  4:29991/1‑149 (MQ=255)                                                                          aCCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAac                                                                         <  4:83234/148‑1 (MQ=255)                                                                            cAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAaca                                                                        >  5:54164/1‑147 (MQ=255)                                                                            cAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACg                                                                      >  5:194548/1‑149 (MQ=255)                                                                             aaCGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTg                                                                                      >  3:102893/1‑132 (MQ=255)                                                                             aaCGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTg                                                                                      <  4:102893/132‑1 (MQ=255)                                                                                      tACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAggtggt                                                            >  1:99907/1‑149 (MQ=255)                                                                                         aaTATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTc                                                         >  5:194153/1‑149 (MQ=255)                                                                                             tCGGGCGGGGTAGCGCGCGTCCGTTTAGCGGCGACAAAGCCGGTACCTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAACCACGGCCAGGTGGTTTCTGTc                                                     <  5:294075/149‑1 (MQ=255)                                                                                                gTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGCCCAGGTGGTTTCTGTCgg                                                   >  6:273866/1‑148 (MQ=255)                                                                                                    tgtTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTaaaa                                              <  3:29991/149‑1 (MQ=255)                                                                                                       tATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTg                                           >  7:234064/1‑149 (MQ=255)                                                                                                        aTCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGc                                          >  7:67407/1‑149 (MQ=255)                                                                                                         tCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGTTAAAATTGCg                                         <  7:181817/149‑1 (MQ=255)                                                                                                         tCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCg                                         <  4:219207/149‑1 (MQ=255)                                                                                                         tCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGc                                          <  6:194548/148‑1 (MQ=255)                                                                                                         tCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCg                                         >  7:111060/1‑149 (MQ=255)                                                                                                             tCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACa                                     >  2:231700/1‑149 (MQ=255)                                                                                                                tCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACGGTGCTGCAGAAACTGGTTGATTAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCt                                  <  8:49915/149‑1 (MQ=255)                                                                                                                tCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGAc                                      >  4:109443/1‑145 (MQ=255)                                                                                                                tCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGAc                                      <  3:109443/145‑1 (MQ=255)                                                                                                                 ccGTTTATCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGACGCCAGCACCGACCGTGCAGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTa                                 >  7:284875/1‑149 (MQ=255)                                                                                                                  cGTTTATCGTCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGCCCTTGCTGTTGAGCCGCCAGCACCGACCGTGCTGAAGAACCTGGTGGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTa                                 <  7:164927/148‑1 (MQ=255)                                                                                                                        tCGGCGACAAAGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAAc                          >  8:249388/1‑149 (MQ=255)                                                                                                                              aCAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAACCACGGCCAggtgg                                                             >  5:58633/1‑108 (MQ=255)                                                                                                                                aaaGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGt                                                      >  7:276505/1‑113 (MQ=255)                                                                                                                                aaaGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGt                                                      <  8:276505/113‑1 (MQ=255)                                                                                                                                aaaGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGAAGAAACTGGATGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTAt                  >  7:228472/1‑149 (MQ=255)                                                                                                                                aaaGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAAtt                                            >  5:109284/1‑123 (MQ=255)                                                                                                                                aaaGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAAtt                                            <  6:109284/123‑1 (MQ=255)                                                                                                                                 aaGCCGGTAACTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATc                 <  7:72270/149‑1 (MQ=255)                                                                                                                                    ccGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCAc               >  1:184417/1‑148 (MQ=255)                                                                                                                                       gTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCaa            >  8:205071/1‑148 (MQ=255)                                                                                                                                          aCTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTa                                 >  6:212633/1‑124 (MQ=255)                                                                                                                                          aCTTCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTa                                 <  5:212633/124‑1 (MQ=255)                                                                                                                                             tCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGt      >  1:267945/1‑148 (MQ=255)                                                                                                                                             tCCAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACTGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGt      >  6:83073/1‑148 (MQ=255)                                                                                                                                              ccAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACAGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGt      <  6:54164/147‑1 (MQ=255)                                                                                                                                              ccAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGt      <  2:81904/147‑1 (MQ=255)                                                                                                                                              ccAGAGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGt      <  2:33716/147‑1 (MQ=255)                                                                                                                                                aGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGaaa  <  5:316060/149‑1 (MQ=255)                                                                                                                                                aGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGa    >  5:363307/1‑147 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GATTGCCTGCCATGAAGAAACTTTCGGTCTGGATAAACTCTACGAACTGTGCGAAATCGCCCGTGAAGAGCTGACCAACGGCGGCTACAATATCGGTCGTGTTATCGCTCGTCCGTTTATCGGCGACAAAGCCGGTAACTTCCAGCGTACCGGTAACCGTCACGACCTGGCTGTTGAGCCGCCAGCACCGACCGTGCTGCAGAAACTGGTTGATGAAAAACACGGCCAGGTGGTTTCTGTCGGTAAAATTGCGGACATCTACGCCAACTGCGGTATCACCAAAAAAGTGAAA  >  NC_000913/4620133‑4620424

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: ‑