Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,524,463	Δ1 bp	5.5%	intergenic (+95/+18)	yncE → / ← ansP	ATP‑binding protein, periplasmic, function unknown/L‑asparagine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,524,463	0	A	.	5.5%	228.7 / 6.4	73	intergenic (+95/+18)	yncE/ansP	ATP‑binding protein, periplasmic, function unknown/L‑asparagine transporter
Reads supporting (aligned to +/- strand):  ref base A (44/25);  new base . (2/2);  total (46/27)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

ACTAAACAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCAAAATA  >  NC_000913/1524315‑1524607                                                                                                                                                      |                                                                                                                                                 aCTAAACAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTa                                                                                                                                                  >  3:198486/1‑149 (MQ=255)  cTAAACAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAg                                                                                                                                                 <  1:312719/149‑1 (MQ=255)     aaCAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCg                                                                                                                                              >  6:72106/1‑149 (MQ=255)      aCAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACAGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCgg                                                                                                                                             >  7:270816/1‑149 (MQ=255)        agcagGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAACACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCggg                                                                                                                                            >  5:230134/1‑148 (MQ=255)        agcagGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCggg                                                                                                                                            >  1:161885/1‑148 (MQ=255)        agcagGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGGAATAGATATGGCGAGGGTGAACACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGc                                                                                                                                           <  5:23701/149‑1 (MQ=255)          cagGAAGATACCCATCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCa                                                                                                                                          <  4:101401/148‑1 (MQ=255)          aagGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGTGTGAAAACGGTGAGACTTGCATGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCa                                                                                                                                          >  3:101401/2‑148 (MQ=255)                gCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTg                                                                                                                                    >  3:216872/1‑148 (MQ=255)                 cTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTACCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGT‑GCCGGGCAAattcaga                                                                                                                                  <  5:253602/149‑5 (MQ=255)                   aCCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAAt                                                                                                                                >  8:172120/1‑149 (MQ=255)                    cccAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCTGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAAt                                                                                                                                >  2:300077/1‑148 (MQ=255)                    cccAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAAt                                                                                                                                >  1:183436/1‑148 (MQ=255)                         ccagACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTc                                                                                                                          <  1:245446/149‑1 (MQ=255)                            gACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGcc                                                                                                                       <  3:150669/149‑1 (MQ=255)                                   agATTCGTATTGCGCTGTAATCGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCCCGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAg                                                                                                                <  3:261141/148‑1 (MQ=255)                                               cgcTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACaa                                                                                                    >  1:248339/1‑149 (MQ=255)                                                    tAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTcc                                                                                               >  7:16965/1‑149 (MQ=255)                                                    tAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTcc                                                                                               >  5:294036/1‑149 (MQ=255)                                                    tAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTAATGCTTAAACACAATTTcc                                                                                               >  5:79285/1‑149 (MQ=255)                                                    tAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACAGCGTATCCGGCACTTTCGTCAGTAATCTGATTTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTc                                                                                                <  8:270816/148‑1 (MQ=255)                                                       tAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGt                                                                                            >  5:263105/1‑149 (MQ=255)                                                       tAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGACGGGCAATTTTGCAATTAACTCGCTGTTTCAGGCTTAAACACAATTTCCTGt                                                                                            <  4:216872/149‑1 (MQ=255)                                                             tGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttct                                                                                      <  2:161885/149‑1 (MQ=255)                                                             tGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttc                                                                                       <  5:72106/148‑1 (MQ=255)                                                             tGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttc                                                                                       <  7:172120/148‑1 (MQ=255)                                                                 gAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTAACACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑CCGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCGTACACCCAATTTCCTGTTTTTCTTCAt                                                                                  >  2:206826/1‑149 (MQ=255)                                                                    ggTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCtt                                                                               >  2:48497/1‑149 (MQ=255)                                                                     gTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCATTTAACTCGCCGTTTCAGGCTTAAACTCAATTTCCTGTTTTTCTTCATCTTc                                                                              >  2:296104/1‑149 (MQ=255)                                                                     gTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTc                                                                              >  2:300874/1‑149 (MQ=255)                                                                      tGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGTAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTcc                                                                             >  5:77440/1‑149 (MQ=255)                                                                        aaaaCGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCTCTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGGAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTc                                                                           <  4:198486/149‑1 (MQ=255)                                                                               tGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACt                                                                    >  1:226784/1‑149 (MQ=255)                                                                                       gcctgcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGt                                                            >  6:114166/1‑149 (MQ=255)                                                                                         ctgcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCGGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGc                                                          >  5:31683/1‑149 (MQ=255)                                                                                           gcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACAGAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCtg                                                        >  2:29436/1‑149 (MQ=255)                                                                                              ttttATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTg                                                                                                                                    <  6:38417/70‑1 (MQ=255)                                                                                                ttATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑TCGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACAGATTTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGaa                                                    <  2:72828/148‑1 (MQ=255)                                                                                                 tATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAAt                                                   <  1:175295/148‑1 (MQ=255)                                                                                                   tCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCAGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCATTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAAttt                                                 >  7:331794/1‑148 (MQ=255)                                                                                                     aCGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCa                                               <  6:31683/148‑1 (MQ=255)                                                                                                     aCGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCa                                               <  5:113188/148‑1 (MQ=255)                                                                                                     aCGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAg                                              >  2:253699/1‑149 (MQ=255)                                                                                                     aCGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAg                                              <  1:48497/149‑1 (MQ=255)                                                                                                             tCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGt                                      >  2:214171/1‑149 (MQ=255)                                                                                                             tCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGt                                      >  3:202141/1‑149 (MQ=255)                                                                                                                ggCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTg                                   >  6:283664/1‑149 (MQ=255)                                                                                                                 gCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTgc                                  <  6:263105/149‑1 (MQ=255)                                                                                                                  cACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTgcg                                 >  3:46553/1‑149 (MQ=255)                                                                                                                    cTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGca                               >  2:300121/1‑149 (MQ=255)                                                                                                                     tttCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGcac                              >  6:123792/1‑149 (MQ=255)                                                                                                                     tttCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGca                               >  2:25883/1‑148 (MQ=255)                                                                                                                     tttCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGca                               >  7:23444/1‑148 (MQ=255)                                                                                                                       tCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGcac                              <  2:248339/147‑1 (MQ=255)                                                                                                                        cGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCacac                            <  2:183436/148‑1 (MQ=255)                                                                                                                         gTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACAcc                           >  1:108190/1‑148 (MQ=255)                                                                                                                         gTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACAcc                           >  8:73493/1‑148 (MQ=255)                                                                                                                          tCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCa                          >  3:139427/1‑148 (MQ=255)                                                                                                                            aGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAAc                       <  1:214171/149‑1 (MQ=255)                                                                                                                            aGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCAGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACTACTGGCGCAGAGCTGTGAATTTCAGCAACACGTTTGCGCACACCAAAc                       >  6:25661/1‑149 (MQ=255)                                                                                                                             gTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAAcc                      >  1:129816/1‑149 (MQ=255)                                                                                                                                 tCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGcc                  >  4:213626/1‑149 (MQ=255)                                                                                                                                 tCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGcc                  >  1:197860/1‑149 (MQ=255)                                                                                                                                 tCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGc                   >  2:292282/1‑148 (MQ=255)                                                                                                                                   caaaattGCCCGGCTACGT‑GCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGGTGTGCGGCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAACCCCGCCta                >  5:36245/6‑149 (MQ=255)                                                                                                                                    gAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCta                >  1:169814/1‑148 (MQ=255)                                                                                                                                     aaTTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCtat               <  1:300077/148‑1 (MQ=255)                                                                                                                                     aaTTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCt                 <  1:253699/146‑1 (MQ=255)                                                                                                                                          gCCCGGCTACGT‑GCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAg         >  2:199246/1‑149 (MQ=255)                                                                                                                                            ccGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGc        >  8:286556/1‑148 (MQ=255)                                                                                                                                             cGGC‑ACGTAGCCGGGTAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCaa      <  6:77440/149‑1 (MQ=255)                                                                                                                                             cGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCaa      <  4:46553/149‑1 (MQ=255)                                                                                                                                                  aCGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCAAAATa  >  2:333771/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                 ACTAAACAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGC‑ACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCAAAATA  >  NC_000913/1524315‑1524607

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