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	98.5%	252.6 / ‑5.1	68	intergenic (‑104/‑166)	ybjL/ybjM	putative transport protein YbjL/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (25/42);  minor base C (1/0);  total (26/42)
Fisher's exact test for biased strand distribution p-value = 3.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCCATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAACGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGAAA  >  NC_000913/889786‑890053                                                                                                                                          |                                                                                                                                   tCCCATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCg                                                                                                                                  >  1:709799/1‑140 (MQ=255)   ccATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGAt                                                                                                                                >  1:504224/1‑140 (MQ=255)        tAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATagtaa                                                                                                                           <  1:9771/140‑1 (MQ=255)           cAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATagtaaagt                                                                                                                        <  2:589726/140‑1 (MQ=255)           cAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATagtaaagt                                                                                                                        >  2:1025580/1‑140 (MQ=255)             aTTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATagtaaagtaa                                                                                                                      <  1:474784/140‑1 (MQ=255)               tCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACa                                                                                                                    <  1:1185915/140‑1 (MQ=255)               tCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACa                                                                                                                    <  1:1419727/140‑1 (MQ=255)               tCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCACAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACa                                                                                                                    <  1:999281/140‑1 (MQ=255)                cGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACaa                                                                                                                   <  1:1063867/140‑1 (MQ=255)                cGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACaa                                                                                                                   <  2:97753/140‑1 (MQ=255)                cGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACaa                                                                                                                   <  1:816138/140‑1 (MQ=255)                 ggCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAat                                                                                                                  >  1:702182/1‑140 (MQ=255)                 ggCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAat                                                                                                                  <  1:1398384/140‑1 (MQ=255)                  gCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAata                                                                                                                 <  2:279270/140‑1 (MQ=255)                        tttATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATAtttt                                                                                                           >  2:1326547/1‑140 (MQ=255)                         ttATATTCACGTAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTa                                                                                                          <  1:607469/140‑1 (MQ=255)                               tCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAgtg                                                                                                    <  1:336921/140‑1 (MQ=255)                                 aCGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTa                                                                                                  >  2:298442/1‑140 (MQ=255)                                  cGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTaa                                                                                                 <  1:729013/140‑1 (MQ=255)                                  cGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTaa                                                                                                 <  2:1092730/140‑1 (MQ=255)                                   ggAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTAcactagt                                                                                                   <  1:1104785/137‑7 (MQ=255)                                   ggAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTAcactagt                                                                                                   >  2:1104785/1‑131 (MQ=255)                                   ggAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTaa                                                                                                 <  2:1211949/139‑1 (MQ=255)                                   ggAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTaa                                                                                                 >  1:1211949/1‑139 (MQ=255)                                   ggAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAAt                                                                                                >  1:501588/1‑140 (MQ=255)                                    gAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATc                                                                                               <  1:1197090/140‑1 (MQ=255)                                       gCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATcaca                                                                                            >  1:762784/1‑140 (MQ=255)                                         tAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACAtt                                                                                          <  1:1025580/140‑1 (MQ=255)                                           aCTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAg                                                                                        >  1:1132544/1‑140 (MQ=255)                                             ttcttGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGt                                                                                      <  1:937692/140‑1 (MQ=255)                                               cttGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTAt                                                                                    >  2:503467/1‑140 (MQ=255)                                                 tGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCa                                                                                  <  2:504224/140‑1 (MQ=255)                                                  gTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCaa                                                                                 <  1:1326547/140‑1 (MQ=255)                                                  gTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCaa                                                                                 <  2:337418/140‑1 (MQ=255)                                                  gTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCaa                                                                                 <  2:1132544/140‑1 (MQ=255)                                                    ttACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACg                                                                               <  2:709799/140‑1 (MQ=255)                                                      aCTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGc                                                                             <  2:1295782/140‑1 (MQ=255)                                                       cTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACg                                                                               <  1:375948/137‑1 (MQ=255)                                                       cTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACg                                                                               >  2:375948/1‑137 (MQ=255)                                                       cTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAACCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCt                                                                            >  2:691970/1‑140 (MQ=255)                                                             aGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGa                                                                      >  1:112418/1‑140 (MQ=255)                                                               cTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACg                                                                               >  1:93855/1‑129 (MQ=255)                                                               cTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACg                                                                               <  2:93855/129‑1 (MQ=255)                                                                  ttGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATc                                                                                   <  2:608712/122‑1 (MQ=255)                                                                  ttGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATc                                                                                   >  1:608712/1‑122 (MQ=255)                                                                   tGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCg                                                                <  2:1317110/140‑1 (MQ=255)                                                                       agaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTgg                                                            <  1:503467/140‑1 (MQ=255)                                                                       agaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTgg                                                            <  1:869752/140‑1 (MQ=255)                                                                        gaaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGc                                                           <  2:112418/140‑1 (MQ=255)                                                                         aaaTGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGcc                                                          >  1:624384/1‑140 (MQ=255)                                                                            tGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGt                                                               <  2:785265/132‑1 (MQ=255)                                                                            tGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGt                                                               >  1:785265/1‑132 (MQ=255)                                                                            tGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCtat                                                       >  2:899570/1‑140 (MQ=255)                                                                            tGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCtat                                                       >  1:1076652/1‑140 (MQ=255)                                                                             ggTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGGGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCtata                                                      <  1:298442/140‑1 (MQ=255)                                                                             ggTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCtata                                                      <  1:193367/140‑1 (MQ=255)                                                                                    tACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTAATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATAATTTACTAGTGTAATCACATTAGGTATCAACGGTTATATGAATTGCGTTGGCCTATATTAGCAt                                               <  1:726563/140‑1 (MQ=255)                                                                                                    aCCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCt                               <  1:986892/140‑1 (MQ=255)                                                                                                    aCCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCt                               <  2:624384/140‑1 (MQ=255)                                                                                                    aCCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATCAGCATGGAATGCGAAGCGGCt                               <  1:392396/140‑1 (MQ=255)                                                                                                                ggCGCAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGa                   >  1:305468/1‑140 (MQ=255)                                                                                                                    cAGATTCATTCAGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGc               >  2:113629/1‑140 (MQ=255)                                                                                                                        ttcattcaGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCcaca           <  1:899570/140‑1 (MQ=255)                                                                                                                            ttcaGCGCACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTgg       >  2:745309/1‑140 (MQ=255)                                                                                                                                 cgcACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGaaa  >  2:388904/1‑140 (MQ=255)                                                                                                                                   cACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCg     >  2:296650/1‑135 (MQ=255)                                                                                                                                   cACCTAAGCGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCg     <  1:296650/135‑1 (MQ=255)                                                                                                                                          |                                                                                                                                   TCCCATTTAACAATTCGGCGACGTTTATATTCACGGAGGCTAACTTCTTGTTTACTAGTAAGCTGTTGAAAGAAATGGTAATTTACGATAATGTTTTTTACCAGAATTCAGGGCGCAGATTCATTCAGCGCACCTAAACGATAGTAAAGTAACAATATATTTTACTAGTGTAATCACATTAGGTATCAACGGCTATATGAATTGCGTTGGCCTATATTAGCATGGAATGCGAAGCGGCTTTATCTTACTGAACGCCACACTGGCGAAA  >  NC_000913/889786‑890053

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: ‑