Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,726,043	G→T	5.5%	intergenic (‑297/+26)	kgtP ← / ← rrfG	alpha‑ketoglutarate transporter/5S ribosomal RNA of rrnG operon

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,726,043	0	G	T	5.5%	326.0 / 6.4	94	intergenic (‑297/+26)	kgtP/rrfG	alpha‑ketoglutarate transporter/5S ribosomal RNA of rrnG operon
Reads supporting (aligned to +/- strand):  ref base G (43/44);  new base T (2/3);  total (45/49)
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.56e-01

AAGGATAAATTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAA  >  NC_000913/2725895‑2726190                                                                                                                                                     |                                                                                                                                                      aaGGATAAATTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCg                                                                                                                                                       >  3:477811/1‑149 (MQ=255)   ggATAAATTTCTGGTAAGGAGGACACGTATAGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGt                                                                                                                                                      <  2:231927/148‑1 (MQ=255)    gATAAATTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTc                                                                                                                                                     <  1:865682/148‑1 (MQ=255)    gATAAATTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCa                                                                                                                                                    >  4:23167/1‑149 (MQ=255)      tAAATTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑gg                                                                                                                                                 <  3:931540/149‑1 (MQ=255)        aaTTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGa                                                                                                                                                <  2:126144/148‑1 (MQ=255)         aTTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGAt                                                                                                                                               <  1:844896/148‑1 (MQ=255)          tttCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATgg                                                                                                                                             >  4:868928/1‑149 (MQ=255)           ttCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATgg                                                                                                                                             >  2:550223/1‑148 (MQ=255)           ttCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATgg                                                                                                                                             <  1:393698/148‑1 (MQ=255)            tCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGcc                                                                                                                                           <  3:516252/149‑1 (MQ=255)             cTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCTGAC‑GGATGGcc                                                                                                                                           <  2:14998/148‑1 (MQ=255)              tGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCTGAC‑GGATGGCCtt                                                                                                                                         <  3:71729/149‑1 (MQ=255)              tGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCtt                                                                                                                                         <  3:724478/149‑1 (MQ=255)               ggTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCC‑TGACGGATGGCCttt                                                                                                                                        <  4:348647/149‑1 (MQ=38)                gTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTc                                                                                                                                       >  4:399622/1‑149 (MQ=255)                 tAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCTGAC‑GGATGGCCtttt                                                                                                                                       <  1:809497/148‑2 (MQ=255)                 tAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCC‑TGACGGATGGCCtttt                                                                                                                                       <  2:67017/148‑2 (MQ=38)                   aggaggACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCa                                                                                                                                   >  3:905145/1‑149 (MQ=255)                       ggaCACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATaa                                                                                                                                >  1:315653/1‑148 (MQ=255)                           acGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTc                                                                                                                                       >  4:163906/1‑138 (MQ=255)                           acGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTc                                                                                                                                       <  3:163906/138‑1 (MQ=255)                            cGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGa                                                                                                                           >  1:448569/1‑148 (MQ=255)                             gTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATGTAGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGTATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATg                                                                                                                         >  4:793626/1‑149 (MQ=255)                               aTGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGc                                                                                                                        >  1:640766/1‑148 (MQ=255)                                   aaGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGc                                                                                                                   >  3:65346/1‑149 (MQ=255)                                    aGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCa                                                                                                                  <  4:883322/149‑1 (MQ=255)                                       gggCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGt                                                                                                                >  2:122787/1‑148 (MQ=255)                                          cAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTcc                                                                                                             <  2:448569/148‑1 (MQ=255)                                           aaGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCt                                                                                                           <  4:477811/149‑1 (MQ=255)                                            aGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTa                                                                                                          <  4:905145/149‑1 (MQ=255)                                               tGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTAct                                                                                                        >  1:143476/1‑148 (MQ=255)                                                 gggAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTc                                                                                                              <  1:293210/140‑1 (MQ=255)                                                 gggAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTc                                                                                                              >  2:293210/1‑140 (MQ=255)                                                 gggAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTAct                                                                                                        <  2:784094/146‑1 (MQ=255)                                                 gggAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTAct                                                                                                        >  1:784094/1‑146 (MQ=255)                                                 gggAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACtctc                                                                                                     <  4:758292/149‑1 (MQ=255)                                                     aGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCAt                                                                                                 <  3:18369/149‑1 (MQ=255)                                                      gCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATg                                                                                                <  3:865425/149‑1 (MQ=255)                                                       ccGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATg                                                                                                <  1:654985/148‑1 (MQ=255)                                                       ccGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATg                                                                                                >  2:807211/1‑148 (MQ=255)                                                        cGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATggg                                                                                              <  4:64995/149‑1 (MQ=255)                                                               gTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGAccc                                                                                       <  4:398639/149‑1 (MQ=255)                                                               gTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGAccc                                                                                       >  3:746065/1‑149 (MQ=255)                                                                  gCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCcac                                                                                    >  4:102692/1‑149 (MQ=255)                                                                   cTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCcaca                                                                                   >  4:868907/1‑149 (MQ=255)                                                                    tGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCcaca                                                                                   >  2:135260/1‑148 (MQ=255)                                                                         cTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTAcc                                                                              >  1:466627/1‑148 (MQ=255)                                                                         cTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTAcc                                                                              >  1:233611/1‑148 (MQ=255)                                                                          tGGCATATGTGGGAGTATAAGACGCGCATCGTCGCATCAGGCATTTTCTTCTGCGCCACTGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTTTCGCATGGGGAGACCCCACACTACCAt                                                                            >  4:719709/1‑149 (MQ=255)                                                                          tGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCAt                                                                            >  4:247031/1‑149 (MQ=255)                                                                            gCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCt                                                                                                                      <  3:446783/105‑1 (MQ=255)                                                                            gCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCt                                                                                                                      >  4:446783/1‑105 (MQ=255)                                                                                    gggAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTAc                                                                   >  1:200701/1‑148 (MQ=255)                                                                                    gggAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACg                                                                  >  4:833769/1‑149 (MQ=255)                                                                                     ggAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTc                                                                                                                                       >  3:936567/1‑80 (MQ=255)                                                                                     ggAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTc                                                                                                                                       <  4:936567/80‑1 (MQ=255)                                                                                           taAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGttt                                                            <  1:655599/148‑1 (MQ=255)                                                                                            aaGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCa                                                          >  3:100757/1‑149 (MQ=255)                                                                                              gACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCa                                                                                                                                   >  1:595912/1‑74 (MQ=255)                                                                                              gACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCa                                                                                                                                   <  2:595912/74‑1 (MQ=255)                                                                                               aCGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACt                                                        >  1:208397/1‑148 (MQ=255)                                                                                               aCGCGCAGCGTCGCATCAGGCATGTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACt                                                        >  1:823203/1‑148 (MQ=255)                                                                                                 gcgcAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                                                     <  4:772860/149‑1 (MQ=255)                                                                                                  cgcAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATACAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCt                                                     >  1:793553/1‑148 (MQ=255)                                                                                                  cgcAGCGTCGCATCAGGCATTTTTTTCTGCACCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTg                                                    <  4:901080/149‑1 (MQ=255)                                                                                                    cAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGa                                                   <  2:793553/148‑1 (MQ=255)                                                                                                       cGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTc                                               <  4:811325/149‑1 (MQ=255)                                                                                                          cGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCgg                                             >  2:168023/1‑148 (MQ=255)                                                                                                           gCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGc                                            <  2:200701/148‑1 (MQ=255)                                                                                                              tCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATgg                                        <  4:100757/149‑1 (MQ=255)                                                                                                                  gCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTc                                               <  2:631810/138‑1 (MQ=255)                                                                                                                  gCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTc                                               >  1:631810/1‑138 (MQ=255)                                                                                                                  gCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTc                                    >  4:178477/1‑149 (MQ=255)                                                                                                                   cATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTTTGAGTTCGGCATGGGGTc                                    >  2:371415/1‑148 (MQ=255)                                                                                                                    aTTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAg                                  <  3:532851/149‑1 (MQ=255)                                                                                                                      ttttttCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGt                                <  4:388581/149‑1 (MQ=255)                                                                                                                      ttttttCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGt                                <  3:571874/149‑1 (MQ=255)                                                                                                                       tttttCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGt                                <  2:937245/148‑1 (MQ=255)                                                                                                                        ttttCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTg                               <  2:143476/148‑1 (MQ=255)                                                                                                                        ttttCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTg                               <  1:135260/148‑1 (MQ=255)                                                                                                                         tttCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTggg                             <  3:178477/149‑1 (MQ=255)                                                                                                                                 ccAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACcg                      <  2:709463/148‑1 (MQ=255)                                                                                                                                    aTGCAAAAAGGCCATCCGTCA‑AGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCgcgc                   >  1:570816/1‑148 (MQ=255)                                                                                                                                        aaaaaGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACgg              >  4:104239/1‑149 (MQ=255)                                                                                                                                         gaaaGGCCATCCTGAC‑GGATGGCCTTTTTGCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACgg              >  3:511113/2‑149 (MQ=11)                                                                                                                                           aaGGCCATCCTGAC‑GGATGGCCTTTC‑GCATAATTTGATGCTTGGGAATTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACTACCGCGCTACGgccg           >  3:808017/1‑149 (MQ=14)                                                                                                                                             ggCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGgccgc          >  1:465812/1‑148 (MQ=255)                                                                                                                                              gCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCa        <  4:341814/149‑1 (MQ=255)                                                                                                                                               ccATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCa        >  2:268794/1‑148 (MQ=255)                                                                                                                                               ccATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAg       <  4:509982/149‑1 (MQ=255)                                                                                                                                                 aTCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAgg      >  2:814500/1‑148 (MQ=255)                                                                                                                                                  tCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGgc             >  1:950213/1‑140 (MQ=255)                                                                                                                                                  tCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGgc             <  2:950213/140‑1 (MQ=255)                                                                                                                                                    cGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCaaa  <  3:102692/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      AAGGATAAATTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAA  >  NC_000913/2725895‑2726190

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