Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,887,033	G→A	G69D (GGC→GAC)	mnmE →	tRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,887,033	0	G	A	100.0%	226.8 / NA	60	G69D (GGC→GAC)	mnmE	tRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/37);  total (23/37)

GGCATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGGCCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGTTTTCCGAACGCGCGTTTCTTA  >  NC_000913/3886885‑3887179                                                                                                                                                     |                                                                                                                                                   ggCATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTACCCTGAcc                                                                                                                                                  <  1:761156/151‑1 (MQ=255)  gCATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGAccc                                                                                                                                                 <  2:53184/151‑1 (MQ=255)         gcgcATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCg                                                                                                                                          >  1:736244/1‑151 (MQ=255)          cgcATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAACCCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGt                                                                                                                                         >  2:112523/1‑151 (MQ=255)                 ccGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCa                                                                                                                                      <  2:65970/147‑1 (MQ=255)                 ccGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCa                                                                                                                                      <  2:380420/147‑1 (MQ=255)                 ccGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCg                                                                                                                                   <  2:448228/150‑1 (MQ=255)                        cAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGCTAAACTGCCTGAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACAGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCTTTCACCGGCGAAGAt                                                                                                                           <  2:821460/151‑1 (MQ=255)                           aGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGAtgtg                                                                                                                        <  2:841294/151‑1 (MQ=255)                            gCCCGTGAGGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGc                                                                                                                       >  2:308956/1‑151 (MQ=255)                            gCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGc                                                                                                                       <  2:407991/151‑1 (MQ=255)                            gCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGc                                                                                                                       <  1:277350/151‑1 (MQ=255)                                   aaGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACt                                                                                                                <  1:771442/151‑1 (MQ=255)                                   aaGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACt                                                                                                                <  2:558634/151‑1 (MQ=255)                                    aGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTg                                                                                                               >  1:53046/1‑151 (MQ=255)                                        gCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAg                                                                                                           <  1:460892/151‑1 (MQ=255)                                           gAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTc                                                                                                        >  1:251542/1‑151 (MQ=255)                                           gAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTc                                                                                                        <  1:472753/151‑1 (MQ=255)                                            aaaCCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCa                                                                                                       <  2:145640/151‑1 (MQ=255)                                              aCCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATg                                                                                                     <  1:191603/151‑1 (MQ=255)                                               ccGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATgg                                                                                                    <  1:372856/151‑1 (MQ=255)                                                       ggTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCgg                                                                                            >  1:545114/1‑151 (MQ=255)                                                         tAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTg                                                                                          >  2:506140/1‑151 (MQ=255)                                                           aaCTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGAt                                                                                        <  1:25681/151‑1 (MQ=255)                                                              tGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATcc                                                                                      >  2:695/1‑150 (MQ=255)                                                               gCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCAt                                                                                                      <  2:786020/133‑1 (MQ=255)                                                               gCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCAt                                                                                                      >  1:786071/1‑133 (MQ=255)                                                               gCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTc                                                                                    <  1:213435/151‑1 (MQ=255)                                                                   aaGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGa                                                                                  <  2:72420/149‑1 (MQ=255)                                                                         cgcTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGt                                                                          <  1:248287/151‑1 (MQ=255)                                                                             aCGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTaa                                                                       <  1:83760/150‑1 (MQ=255)                                                                             aCGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTa                                                                        <  1:894168/149‑1 (MQ=255)                                                                              cGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTaaaa                                                                     <  2:897086/151‑1 (MQ=255)                                                                               gCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAAc                                                                    <  2:887186/151‑1 (MQ=255)                                                                               gCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGGGATCCTCGACCTGCTGTTaaaa                                                                     >  2:199621/1‑150 (MQ=255)                                                                                          tCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGGGATCCTCGGACTGCTGTTAAAACGCATTCTGAc                                                          >  2:818396/1‑150 (MQ=255)                                                                                           ccGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGAcc                                                         <  1:144607/150‑1 (MQ=255)                                                                                            cGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGAcc                                                         >  1:905637/1‑149 (MQ=255)                                                                                               ttAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGGGATCCTTGACCTGCTGTTAAAACGAATTCTGACCATTcc                                                    >  1:254002/1‑151 (MQ=255)                                                                                                 aaaGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCg                                                  <  1:688344/151‑1 (MQ=255)                                                                                                    gACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGcc                                               <  2:8370/151‑1 (MQ=255)                                                                                                        ccGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCg                                           <  1:268682/151‑1 (MQ=255)                                                                                                         cGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCgg                                          >  1:716988/1‑151 (MQ=255)                                                                                                          gACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGa                                         <  1:200527/151‑1 (MQ=255)                                                                                                                  cGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGc                                 <  1:332322/151‑1 (MQ=255)                                                                                                                   gTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCCCGACCTGCTGTTAAAACGCATTCTGGCCATTCCCGGCCTGCGGATTGCt                                    >  1:455484/1‑147 (MQ=255)                                                                                                                        cGATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGt                           <  1:579766/151‑1 (MQ=255)                                                                                                                         gATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTg                          <  2:793389/151‑1 (MQ=255)                                                                                                                         gATCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTg                          >  1:138631/1‑151 (MQ=255)                                                                                                                          aTCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTgg                            >  2:839471/1‑148 (MQ=255)                                                                                                                          aTCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTgg                            <  1:839538/148‑1 (MQ=255)                                                                                                                          aTCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTg                          >  1:84008/1‑150 (MQ=255)                                                                                                                          aTCAGGGGATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTg                          >  1:218428/1‑150 (MQ=255)                                                                                                                              ggggATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGttt                     <  1:112525/151‑1 (MQ=255)                                                                                                                               gggATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGtttt                    <  1:185967/151‑1 (MQ=255)                                                                                                                                 gATTGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCCGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGGTTGCTCGCCCTGGTGAGTTTTcc                  >  1:858468/1‑151 (MQ=255)                                                                                                                                    tGCGCTATGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGGGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGTTTTCCg                 >  1:156379/1‑149 (MQ=255)                                                                                                                                           tGGTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGTTTTCCGAAcg             >  1:595584/1‑146 (MQ=255)                                                                                                                                             gTTCCCTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGTTTTCCGAACGCGCGttt      <  1:655293/151‑1 (MQ=255)                                                                                                                                                  cTGACCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGTTTTCCGCACGCGCGTTTCTTa  >  1:566421/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GGCATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGACGGCAGCGTGCTCGATCAGGGGATTGCGCTATGGTTCCCTGGCCCGAACTCGTTCACCGGCGAAGATGTGCTGGAACTGCAAGGTCATGGCGGTCCGGTGATCCTCGACCTGCTGTTAAAACGCATTCTGACCATTCCCGGCCTGCGGATTGCTCGCCCTGGTGAGTTTTCCGAACGCGCGTTTCTTA  >  NC_000913/3886885‑3887179

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