Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,120,913	C→A	6.3%	G15W (GGG→TGG)	bssS ←	biofilm regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,120,913	0	C	A	6.3%	227.4 / 2.7	63	G15W (GGG→TGG)	bssS	biofilm regulator
Reads supporting (aligned to +/- strand):  ref base C (36/23);  new base A (2/2);  total (38/25)
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.48e-01

TTGTCTGGCAACATCAGTGGTTAACCAGAGCGTCTGACCAACTTCAGTCCCTTCCTGCTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTTTATCTGAAGCAACTTCAGTATACCAATAAGTAAGGGTATTCA  >  NC_000913/1120779‑1121052                                                                                                                                       |                                                                                                                                            ttGTCTGGCAACATCAGTGGTTAACCAGAGCGTCTGACCAACTTCAGTCCCTTCCTGCTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCTCCCTACGa                                                                                                                                        >  2:1684631/1‑140 (MQ=255)      tGGCAACATCAGTGGTTAACCAGAGCGTCTGACCAACTTCAGTCCCTTCCTGCTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGa                                                                                                                                   >  2:526908/1‑140 (MQ=255)           aCATCAGTGGTTAACCAGAGCGTCTGACCAACTTCAGTCCCTTCCTGCTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGt                                                                                                                              <  1:1537419/140‑1 (MQ=255)                 gTGGTTAACCAGAGCGTCTGACCAACTTCAGTCCCTTCCTGCTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAAt                                                                                                                        <  1:1553792/140‑1 (MQ=255)                   ggTTAACCAGAGCGTCTGACCAACTTCAGTCCCTTCCTGCTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGa                                                                                                                      <  1:1771352/140‑1 (MQ=255)                       aaCCAGAGCGTCTGACCAACTTCAGTCCCTTCCTGCTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTc                                                                                                                  >  1:1384884/1‑140 (MQ=255)                       aaCCAGAGCGTCTGACCAACTTCAGTCCCTTCCTGCTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTc                                                                                                                  >  1:1875044/1‑140 (MQ=255)                             agCGTCTGACCAACTTCAGTCCCTTCCTGCTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCattatt                                                                                                            >  1:905947/1‑140 (MQ=255)                             agCGTCTGACCAACTTCAGTCCCTTCCTGCTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCattatt                                                                                                            >  2:552673/1‑140 (MQ=255)                             agCGTCTGACCAACTTCAGTCCCTTCCTGCTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCATACGAGCGGATGAGTCTGAATGACTTCattatt                                                                                                            >  1:1590689/1‑140 (MQ=255)                             agCGTCTGACCAACTTCAGTCCCTTCCTGCTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCATACGAGCGGATGAGTCTGAATGACTTCattatt                                                                                                            >  1:1590351/1‑140 (MQ=255)                                  cTGACCAACTTCAGTCCCTTCCTGCTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTAttttttt                                                                                                       >  2:518933/1‑140 (MQ=255)                                  cTGACCAACTTCAGTCCCTTCCTGCTCGGACTTATGTGGGGTCTGGTAGTGCAAACGCAACATCGGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTAttttttt                                                                                                       >  2:519140/1‑140 (MQ=255)                                         aCTTCAGTCCCTTCCTGCTCGTACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTACCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTccgacat                                                                                                >  2:216119/1‑135 (MQ=255)                                          cTTCAGTCCCTTCCTGCTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATg                                                                                               <  1:479080/140‑1 (MQ=255)                                          cTTCAGTCCCTTCCTGCTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATg                                                                                               <  1:1683794/140‑1 (MQ=255)                                               gTCCCTTCCTGCTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGcccc                                                                                          <  2:1673320/140‑1 (MQ=255)                                               gTCCCTTCCTGCTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGcccc                                                                                          <  2:1565771/140‑1 (MQ=255)                                               gTCCCTTCCTGCTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGcccc                                                                                          >  2:1562718/1‑140 (MQ=255)                                                     tCCTGCTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAAtt                                                                                    >  1:665746/1‑140 (MQ=255)                                                       cTGCTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCATACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCg                                                                                  <  2:1590351/140‑1 (MQ=255)                                                       cTGCTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCATACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCg                                                                                  <  2:1590689/140‑1 (MQ=255)                                                          cTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCgtgt                                                                               <  1:28349/140‑1 (MQ=255)                                                                 ttATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGa                                                                        <  1:1396352/140‑1 (MQ=255)                                                                 ttATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGa                                                                        <  1:660759/140‑1 (MQ=255)                                                                 ttATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGa                                                                        <  1:754043/140‑1 (MQ=255)                                                                 ttATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGa                                                                        <  1:1653470/140‑1 (MQ=255)                                                                 ttATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGa                                                                        <  1:1572975/140‑1 (MQ=255)                                                                     ttGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGTGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAAc                                                                    <  2:722393/140‑1 (MQ=255)                                                                     ttGGGGTCTGGTAGTGAAACCGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTATTTTCCATCATGGCCCCCTAAGTCGTGTTACTAGACAAc                                                                    <  2:723246/140‑1 (MQ=255)                                                                         ggTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTa                                                                >  2:275533/1‑140 (MQ=255)                                                                                     cAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATTAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCaa                                                    >  2:1206653/1‑140 (MQ=255)                                                                                     cAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATTAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCaa                                                    >  2:1207229/1‑140 (MQ=255)                                                                                         cGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCg                                                >  2:1177168/1‑140 (MQ=255)                                                                                             aCATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGtttt                                            >  1:1073083/1‑140 (MQ=255)                                                                                             aCATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGtttt                                            >  1:1578669/1‑140 (MQ=255)                                                                                             aCATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGtttt                                            >  1:1781504/1‑140 (MQ=255)                                                                                                tCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTTTAt                                         <  2:84084/140‑1 (MQ=255)                                                                                                    cgcATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTTTATCTGa                                     >  2:1755610/1‑140 (MQ=255)                                                                                                    cgcATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTTTATCTGa                                     >  1:369273/1‑140 (MQ=255)                                                                                                    cgcATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAAGGCGTTTTTATCTGa                                     >  1:1420299/1‑140 (MQ=255)                                                                                                     gcATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTTTATCTGaa                                    <  2:281389/140‑1 (MQ=255)                                                                                                                aTCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACCTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGTGGTTCAATGCGTTTTTATCTGAAGCAACTTCAGt                         <  1:1072770/140‑1 (MQ=255)                                                                                                                 tCAACGGTGCTGATGTCCCACGCTACGAGCGGATGAGTCTGAATGTCTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCATTGCGTTTTTATCTGTTGCAACTTCAGta                        >  1:1372304/1‑140 (MQ=255)                                                                                                                 tCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTTTATCTGAAGCAACTTCAGta                        >  1:654331/1‑140 (MQ=255)                                                                                                                 tCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTTTATCTGAAGCAACTTCAGta                        >  1:37689/1‑140 (MQ=255)                                                                                                                 tCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTTTATCTGAAGCAACTTCAGta                        >  1:37653/1‑140 (MQ=255)                                                                                                                 tCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTTTATCTGAAGCAACTTCAGta                        >  1:36936/1‑140 (MQ=255)                                                                                                                 tCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTTTATCTGAAGCAACTTCAGta                        >  1:1499255/1‑140 (MQ=255)                                                                                                                 tCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAAGGCGTTTTTATCTGAAGCAACTTCAGta                        >  1:278193/1‑140 (MQ=255)                                                                                                                 tCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCAAAATCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTATATCTGAAGCAACTTCAGtc                        >  1:1902329/1‑139 (MQ=255)                                                                                                                 tCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCAATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTTTATCTGAAGCAACTTCAGta                        >  1:1195023/1‑140 (MQ=255)                                                                                                                   aaCGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTTTATCTGAAGCAACTTCAGtata                      >  1:102082/1‑140 (MQ=255)                                                                                                                   aaCGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTTTATCTGAAGCAACTTCAGtata                      >  1:101977/1‑140 (MQ=255)                                                                                                                      ggTGCAGATGTCCCACCCCACGAGCGGATGATTGTGAATGACTTATTTATTTTTTTCCATCATGGCCTCCTAATTCGTTGTACTAGAAAACGTAATTCGAGGTTCAATGCGTTTTTATCTGAAGGAACTTAAGTATAcct                   <  2:1016465/140‑2 (MQ=255)                                                                                                                        tGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTAATTCGAGGTTCAATGCGTTTTTATCTGAATCAACTTCAGTATACCAAt                 >  2:1706402/1‑140 (MQ=255)                                                                                                                        tGCTGATGTCACACCCTACGTGCGGCTGAGTCTGAATGACTTCAGTATTTTTTTCCATCATGGCCCGCTAATTCGTGTTACTAGTCAACGTTAGTCGAGGGTCAAGGCGTGTTTATCGGTAGCAACTTCAGTATCCCatg                 >  2:521213/1‑138 (MQ=255)                                                                                                                           tGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTTTATCTGAAGCAACTTCAGTATACCAAtaag              <  2:1384884/140‑1 (MQ=255)                                                                                                                             aTGTCCCTCCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTTTATCTGAAGCAACTTCAGTATACCAAtaagta            <  1:1684631/140‑1 (MQ=255)                                                                                                                             aTGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTTTATCTGAAGCAACTTCAGTATACCAAtaagta            <  2:905947/140‑1 (MQ=255)                                                                                                                             aTGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTGCCATCATGGCCCCCTAATTCGTGTTACTAGTCAACGTTATTCGAGGTTCAATGCGTTTTTATCTGAAGCAACTTCAGTATACCAAtaagta            <  2:961803/140‑1 (MQ=255)                                                                                                                                  ccACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTTTATCTGAAGCAACTTCAGTATACCAATAAGTAAGGGt       >  1:526449/1‑140 (MQ=255)                                                                                                                                  ccACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTTTATCTGAAGCAACTTCAGTATACCAATAAGTAAGGGt       >  1:1822811/1‑140 (MQ=255)                                                                                                                                       cTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTTTATCTGAATCAACTTCAGTATACCAATAAGTAAGGGTATTCa  >  1:922657/1‑140 (MQ=255)                                                                                                                                       cTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTTTATCTGAAGCAACTTCAGTATACCAATAAGTAAGGGTATTCa  >  2:1345674/1‑140 (MQ=255)                                                                                                                                       cTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTTTATCTGAAGCAACTTCAGTATACCAATAAGTAAGGGTATTCa  >  1:922109/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                            TTGTCTGGCAACATCAGTGGTTAACCAGAGCGTCTGACCAACTTCAGTCCCTTCCTGCTCGGACTTATTTGGGGTCTGGTAGTGCAAACGCAACATCAGCGCATCATAGCTATCAACGGTGCTGATGTCCCACCCTACGAGCGGATGAGTCTGAATGACTTCATTATTTTTTTCCATCATGGCCCCCTAATTCGTGTTACTAGACAACGTTATTCGAGGTTCAATGCGTTTTTATCTGAAGCAACTTCAGTATACCAATAAGTAAGGGTATTCA  >  NC_000913/1120779‑1121052

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