Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,635,159	C→T	G51D (GGC→GAC)	rob ←	right oriC‑binding transcriptional activator, AraC family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,635,159	0	C	T	100.0%	233.7 / NA	61	G51D (GGC→GAC)	rob	right oriC‑binding transcriptional activator, AraC family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/25);  total (36/25)

AGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTCCAGCCAGATTAAAAGGTCGCGAATAATGCCGGCCTGAT  >  NC_000913/4635014‑4635306                                                                                                                                                  |                                                                                                                                                    aGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAAt                                                                                                                                                 >  2:420070/1‑150 (MQ=255)      gggCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCAt                                                                                                                                            <  1:666149/150‑1 (MQ=255)            aCTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAg                                                                                                                                      >  1:110350/1‑150 (MQ=255)             cTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGt                                                                                                                                     >  2:501591/1‑150 (MQ=255)              tGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTg                                                                                                                                    >  2:422157/1‑150 (MQ=255)                 ttcttGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACAt                                                                                                                                >  1:626992/1‑151 (MQ=255)                   cttGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCt                                                                                                                              >  1:20817/1‑151 (MQ=255)                     tGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATc                                                                                                                               <  1:599660/148‑1 (MQ=255)                        aTGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTaa                                                                                                                          >  1:495211/1‑150 (MQ=255)                         tGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAAc                                                                                                                        >  1:661454/1‑151 (MQ=255)                              gggTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATtc                                                                                                                    <  1:294541/150‑1 (MQ=255)                                  aaaTGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGt                                                                                                               >  2:303692/1‑151 (MQ=255)                                    aTGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTa                                                                                                              >  1:387142/1‑150 (MQ=255)                                          gTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGcc                                                                                                        >  1:438262/1‑150 (MQ=255)                                           ttGAGAGTCGAAGCGGTATTGCAGCGCGTTGCCAAGAATCGGCCGCGCAGTCAGGCGTAGTGGGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAg                                                                                                            <  2:91036/145‑1 (MQ=255)                                                agTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTgg                                                                                                  >  2:420265/1‑150 (MQ=255)                                                 gTCGAGGGGGTATGGCAGCGCGATGCCCAAAACGGGCCGCGCGGTCGGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGaa                                                                                                <  2:604284/151‑1 (MQ=255)                                                 gTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGa                                                                                                 >  2:488377/1‑150 (MQ=255)                                                      aGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAAcc                                                                                           >  2:619974/1‑151 (MQ=255)                                                       gCGGTATTGCAGCGCGATGCCCAGAATCGGACGCGCAGTCAGGCTTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCt                                                                                          <  2:748721/151‑1 (MQ=255)                                                         ggTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGc                                                                                        >  2:133927/1‑151 (MQ=255)                                                         ggTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGc                                                                                        >  2:368146/1‑151 (MQ=255)                                                            aTTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCttt                                                                                     >  2:469274/1‑151 (MQ=255)                                                            aTTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCttt                                                                                     >  2:485478/1‑151 (MQ=255)                                                             ttGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCtt                                                                                      <  1:362738/149‑1 (MQ=255)                                                              tGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTc                                                                                    <  1:663710/150‑1 (MQ=255)                                                                   cgcgATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCt                                                                              <  2:275125/151‑1 (MQ=255)                                                                     cgATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTAc                                                                            >  2:760498/1‑151 (MQ=255)                                                                     cgATGTCCAGAACCGGCCGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTa                                                                             <  1:396438/150‑1 (MQ=255)                                                                       aTGTCCAAAATCGGACGCGCAGTCGGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACAt                                                                          <  1:124802/151‑1 (MQ=255)                                                                         gTCCAGAATCGGACGCGCAGTCGGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACAtt                                                                         <  2:86479/150‑1 (MQ=255)                                                                             aGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCg                                                                     >  2:679650/1‑150 (MQ=255)                                                                                   gggggcgcAGTCGGGCGTAGGGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACa                                                              <  1:712087/147‑1 (MQ=255)                                                                                    gACGCGCGTTCAGGCGGAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGa                                                                <  2:661376/148‑1 (MQ=255)                                                                                       gcgcAGTCGGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAg                                                             <  1:420116/148‑1 (MQ=255)                                                                                       gcgcAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAggg                                                           <  1:384550/150‑1 (MQ=255)                                                                                          cAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTg                                                       >  1:23727/1‑151 (MQ=255)                                                                                                ggCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATcc                                                   >  2:727582/1‑149 (MQ=255)                                                                                                ggCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATcc                                                   >  1:358327/1‑149 (MQ=255)                                                                                                    tAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATg                                              >  1:103045/1‑150 (MQ=255)                                                                                                      gTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGAcc                                           >  1:492733/1‑151 (MQ=255)                                                                                                          gACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTcc                                       >  2:62421/1‑151 (MQ=255)                                                                                                           aCCGCCGATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTcc                                       <  2:188001/150‑1 (MQ=255)                                                                                                                cgATTTCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTccagcc                                   >  2:364277/1‑149 (MQ=255)                                                                                                                     tCGACAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTCCAGCCAGATTa                             >  2:216723/1‑150 (MQ=255)                                                                                                                        aCAAACGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTCCAGCCAGATTaaa                           >  2:619952/1‑149 (MQ=255)                                                                                                                          aaaCGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTCCAGCCAGATTAAAAGGt                       <  2:704487/151‑1 (MQ=255)                                                                                                                          aaaCGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTCCAGCCAGATTAAAAGGt                       >  1:706574/1‑151 (MQ=255)                                                                                                                          aaaCGACGAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTCCAGCCAGATTAAAAGGt                       <  2:702792/151‑1 (MQ=255)                                                                                                                            acgacgAGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTCCAGCCAGATTAAAAGGTcg                     <  1:189651/151‑1 (MQ=255)                                                                                                                             cgacgaGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTCCAGCCAGATTAAAAGGTcgc                    <  2:181633/151‑1 (MQ=255)                                                                                                                              gacgaGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTCCAGCCAGATTAAAAGGTcgcg                   >  1:586694/1‑151 (MQ=255)                                                                                                                              gacgaGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGCAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTCCAGCCAGATTAAAAGGTcgcg                   <  1:30205/151‑1 (MQ=255)                                                                                                                                cgaGCACGAATATACGCGTCAATAGCATGGCCAGAGAAATAGTTAAACATTCACTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTCCAGCCAGATTAAAAGGTCGCGa                  >  2:784538/1‑150 (MQ=255)                                                                                                                                 gaGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTCCAGCCAGATTAAAAGGTCGCGaat                >  2:607976/1‑151 (MQ=255)                                                                                                                                 gaGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTCCAGCCAGATTAAAAGGTCGCGaat                >  1:626597/1‑151 (MQ=255)                                                                                                                                  aGCACGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTCCAGCCAGATTAAAAGGTCGCGaata               <  2:490176/151‑1 (MQ=255)                                                                                                                                     aCGAATATACGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTCCAGCCAGATTAAAAGGTCGCGAATAATg            >  2:784858/1‑151 (MQ=255)                                                                                                                                           ataCGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTCCAGCCAGATTAAAAGGTCGCGAATAATGCCGGcc      <  2:564296/151‑1 (MQ=255)                                                                                                                                           ataCGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTCCAGCCAGATTAAAAGGTCGCGAATAATGCCGGcc      >  2:729461/1‑151 (MQ=255)                                                                                                                                             aCGCGTCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTCCAGCCAGATTAAAAGGTCGCGAATAATGCCGGCCt     <  1:783731/150‑1 (MQ=255)                                                                                                                                               gcgtCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTCCAGCCAGATTAAAAGGTCGCGAATAATGCCGGCCTGAt  >  1:551380/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                    AGTCTGGGCAAACTGCTTCTTGAATGCGCGGGTAAATGTCTGTTGAGAGTCGAAGCGGTATTGCAGCGCGATGTCCAGAATCGGACGCGCAGTCAGGCGTAGTGCGACCGCCGATTTCGACAAACGACGAGCACGAATATACGCGCCAATAGCATGGCCAGTGACATCTTTAAACATTCTCTGTAAGTGCCACTTGGAATAACCTGCTTTCGCCGCTACATTGTCGAGCGACAGGGGCTGATCCAGATGACCTTCCAGCCAGATTAAAAGGTCGCGAATAATGCCGGCCTGAT  >  NC_000913/4635014‑4635306

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