Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	889,923	A→G	intergenic (‑104/‑166)	ybjL ← / → ybjM	putative transport protein YbjL/putative inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	889,923	0	A	G	100.0%	265.1 / NA	68	intergenic (‑104/‑166)	ybjL/ybjM	putative transport protein YbjL/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/42);  total (26/42)

CCATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAACGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTG  >  NC_000913/889788‑890047                                                                                                                                        |                                                                                                                             ccATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGAt                                                                                                                          <  2:1185415/140‑1 (MQ=255) ccATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGAt                                                                                                                          <  2:149674/140‑1 (MQ=255) ccATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGAt                                                                                                                          >  2:644114/1‑140 (MQ=255) ccATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGAt                                                                                                                          <  2:306037/140‑1 (MQ=255) ccATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGAt                                                                                                                          <  1:728908/140‑1 (MQ=255)   aTTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATag                                                                                                                        >  2:964630/1‑140 (MQ=255)   aTTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATag                                                                                                                        >  2:805834/1‑140 (MQ=255)   aTTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATag                                                                                                                        >  2:765443/1‑140 (MQ=255)     ttAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATagta                                                                                                                      <  2:395404/140‑1 (MQ=255)     ttAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATagta                                                                                                                      >  2:1153596/1‑140 (MQ=255)         cAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATagtaaagt                                                                                                                  <  2:711928/140‑1 (MQ=255)         cAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATagtaaagt                                                                                                                  <  1:46113/140‑1 (MQ=255)              cGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACaa                                                                                                             >  1:732577/1‑140 (MQ=255)              cGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACaa                                                                                                             <  1:765443/140‑1 (MQ=255)              cGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACaa                                                                                                             <  1:1018736/140‑1 (MQ=255)                gCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAata                                                                                                           <  1:1046647/140‑1 (MQ=255)                    cGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATAtt                                                                                                       >  2:1015741/1‑140 (MQ=255)                       ttATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTa                                                                                                    <  1:25272/140‑1 (MQ=255)                       ttATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTa                                                                                                    <  2:1254807/140‑1 (MQ=255)                            ttCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAgt                                                                                               <  2:1162038/140‑1 (MQ=255)                                 ggAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAAt                                                                                          <  1:350699/140‑1 (MQ=255)                                 ggAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAAt                                                                                          <  2:884464/140‑1 (MQ=255)                                 ggAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAAt                                                                                          <  1:1336108/140‑1 (MQ=255)                                    ggCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAAt                                                                                          <  1:428419/137‑1 (MQ=255)                                    ggCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAAt                                                                                          >  2:428419/1‑137 (MQ=255)                                     gCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATcaca                                                                                      <  1:644114/140‑1 (MQ=255)                                        aaCTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATcac                                                                                       <  1:978657/136‑1 (MQ=255)                                        aaCTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATcac                                                                                       >  2:978657/1‑136 (MQ=255)                                         aCTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAg                                                                                  >  1:568940/1‑140 (MQ=255)                                              ttGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATc                                                                             <  1:692190/140‑1 (MQ=255)                                                     cTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATcaca                                                                                      >  2:87275/1‑124 (MQ=255)                                                     cTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATcaca                                                                                      <  1:87275/124‑1 (MQ=255)                                                        gTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCtata                                                                   >  1:457115/1‑140 (MQ=255)                                                           aGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGa                                                                >  1:843951/1‑140 (MQ=255)                                                           aGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGa                                                                >  2:369235/1‑140 (MQ=255)                                                            gCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGaa                                                               <  1:1153596/140‑1 (MQ=255)                                                                ttGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGc                                                           >  2:1047759/1‑140 (MQ=255)                                                                ttGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGc                                                           >  1:838399/1‑140 (MQ=255)                                                                ttGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGc                                                           >  2:353073/1‑140 (MQ=255)                                                                 tGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCg                                                          <  1:1195152/140‑1 (MQ=255)                                                                 tGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCg                                                          <  1:805834/140‑1 (MQ=255)                                                                 tGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCg                                                          <  1:353073/140‑1 (MQ=255)                                                                  gaaagaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGt                                                         <  2:658212/140‑1 (MQ=255)                                                                  gaaagaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGt                                                         <  2:304498/140‑1 (MQ=255)                                                                     agaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTgg                                                      <  2:96948/140‑1 (MQ=255)                                                                      gaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGc                                                     <  2:732577/140‑1 (MQ=255)                                                                      gaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGc                                                     <  1:1047759/140‑1 (MQ=255)                                                                      gaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGc                                                     <  1:369235/140‑1 (MQ=255)                                                                      gaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGc                                                     <  2:1097478/140‑1 (MQ=255)                                                                      gaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGc                                                     <  1:964630/140‑1 (MQ=255)                                                                       aaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGcc                                                    <  2:46215/140‑1 (MQ=255)                                                                          tGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGaa                                                               >  2:1341260/1‑126 (MQ=255)                                                                          tGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGaa                                                               <  1:1341260/126‑1 (MQ=255)                                                                          tGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCta                                                  >  2:204657/1‑139 (MQ=255)                                                                          tGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCta                                                  <  1:204657/139‑1 (MQ=255)                                                                           ggTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCtata                                                <  2:568940/140‑1 (MQ=255)                                                                                             tttttACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAg                              <  1:133918/140‑1 (MQ=255)                                                                                                 tACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGaa                               <  1:311603/135‑1 (MQ=255)                                                                                                 tACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGaa                               >  2:311603/1‑135 (MQ=255)                                                                                                 tACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGc                          >  1:942218/1‑140 (MQ=255)                                                                                                 tACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGc                          <  2:942218/140‑1 (MQ=255)                                                                                                      gAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTAt                     >  2:1302261/1‑140 (MQ=255)                                                                                                                  cAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACt               <  2:1025864/134‑1 (MQ=255)                                                                                                                  cAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACt               >  1:1025864/1‑134 (MQ=255)                                                                                                                  cAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGc         >  2:469596/1‑140 (MQ=255)                                                                                                                   aGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGaa            >  2:544729/1‑136 (MQ=255)                                                                                                                   aGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGaa            <  1:544729/136‑1 (MQ=255)                                                                                                                         attcaGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTg  >  1:1246064/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                             CCATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAACGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTG  >  NC_000913/889788‑890047

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