Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,726,046	A→C	5.5%	intergenic (‑300/+23)	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,046	0	A	C	5.5%	315.0 / 9.5	91	intergenic (‑300/+23)	kgtP/rrfG	alpha‑ketoglutarate transporter/5S ribosomal RNA of rrnG operon
Reads supporting (aligned to +/- strand):  ref base A (42/44);  new base C (2/3);  total (44/47)
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.21e-01

GATAAATTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAA  >  NC_000913/2725898‑2726190                                                                                                                                                     |                                                                                                                                                   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                                                                                                                                             <  1:393698/148‑1 (MQ=255)        ttCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATgg                                                                                                                                             >  2:550223/1‑148 (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                                                                                                        >  1:784094/1‑146 (MQ=255)                                              gggAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTAct                                                                                                        <  2:784094/146‑1 (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                                                                                       >  3:746065/1‑149 (MQ=255)                                                            gTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGAccc                                                                                       <  4:398639/149‑1 (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:233611/1‑148 (MQ=255)                                                                      cTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTAcc                                                                              >  1:466627/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                                                                                                                      >  4:446783/1‑105 (MQ=255)                                                                         gCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCt                                                                                                                      <  3:446783/105‑1 (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)                                                                                                                                                     |                                                                                                                                                   GATAAATTTCTGGTAAGGAGGACACGTATGGAAGTGGGCAAGTTGGGGAAGCCGTATCCGTTGCTGAATCTGGCATATGTGGGAGTATAAGACGCGCAGCGTCGCATCAGGCATTTTTTTCTGCGCCAATGCAAAAAGGCCATCCGTCA‑GGATGGCCTTTC‑GCATAATTTGATGCCTGGCAGTTCCCTACTCTCGCATGGGGAGACCCCACACTACCATCGGCGCTACGGCGTTTCACTTCTGAGTTCGGCATGGGGTCAGGTGGGACCACCGCGCTACGGCCGCCAGGCAAA  >  NC_000913/2725898‑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: ‑