Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,253,896	C→T	6.0%	R686Q (CGA→CAA)	ycgV ←	putative adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,253,896	0	C	T	6.0%	211.6 / 3.0	66	R686Q (CGA→CAA)	ycgV	putative adhesin
Reads supporting (aligned to +/- strand):  ref base C (28/34);  new base T (2/2);  total (30/36)
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.69e-01

AACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAGGCGTGGGAGCGGGTTCAGGCGTAGGATCGGGGTTGAC  >  NC_000913/1253749‑1254038                                                                                                                                                    |                                                                                                                                               aaCGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCg                                                                                                                                               >  7:154682/1‑149 (MQ=255)    ggCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAg                                                                                                                                            <  8:156775/149‑1 (MQ=255)    ggCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAg                                                                                                                                            <  4:177303/149‑1 (MQ=255)       aTTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTTGCAGGt                                                                                                                                          >  3:265699/1‑148 (MQ=255)       aTTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTTGCAGGTc                                                                                                                                         >  3:244132/1‑149 (MQ=255)         tACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGTTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACGCTGATTTCGCAGGTCAc                                                                                                                                       >  8:255162/1‑149 (MQ=255)          aCATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCAc                                                                                                                                       <  6:84354/148‑1 (MQ=255)             tCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCAt                                                                                                                                   >  1:82174/1‑149 (MQ=255)                  gagaCGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGTCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGtt                                                                                                                              <  5:11955/149‑1 (MQ=255)                           aTCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGc                                                                                                                      <  8:272473/148‑1 (MQ=255)                              cccACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTa                                                                                                                   <  3:291848/148‑1 (MQ=255)                              cccACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTAc                                                                                                                  <  8:154682/149‑1 (MQ=255)                               ccaccaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACg                                                                                                                 >  3:11300/1‑149 (MQ=255)                               ccaccaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACg                                                                                                                 >  8:316378/1‑149 (MQ=255)                                caccaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAATGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACgg                                                                                                                >  1:9707/1‑149 (MQ=255)                                caccaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACgg                                                                                                                >  7:181786/1‑149 (MQ=255)                                caccaAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACgg                                                                                                                <  7:74623/149‑1 (MQ=255)                                   caAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGttt                                                                                                             <  8:181786/149‑1 (MQ=255)                                      aCTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTc                                                                                                           >  5:274277/1‑148 (MQ=255)                                      aCTCGATACCGCTGTAAACCATGTCAAAGCAGCTCAGTCTGCCACTGGCAAAGGAGTCCACGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCAATCCGTTCCATCAGCGTACGGTTCTCa                                                                                                          <  2:70386/149‑1 (MQ=255)                                         ggATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTc                                                                                                           >  7:105618/1‑145 (MQ=255)                                         ggATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTc                                                                                                           <  8:105618/145‑1 (MQ=255)                                          gATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACa                                                                                                       >  3:307055/1‑148 (MQ=255)                                          gATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGACACTGGCAAAGGAGTCAAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACa                                                                                                       >  6:226721/1‑148 (MQ=255)                                                  tGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCaa                                                                                               >  7:40648/1‑148 (MQ=255)                                                  tGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCaa                                                                                               >  7:136308/1‑148 (MQ=255)                                                    tAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACg                                                                                                                 <  8:51501/128‑1 (MQ=255)                                                    tAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACg                                                                                                                 >  7:51501/1‑128 (MQ=255)                                                    tAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATaa                                                                                            <  8:38734/149‑1 (MQ=255)                                                      aCCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAg                                                                                           >  6:21035/1‑148 (MQ=255)                                                        ccATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGtt                                                                                                  >  2:340813/1‑139 (MQ=255)                                                        ccATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGtt                                                                                                  <  1:340813/139‑1 (MQ=255)                                                         cATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAAc                                                                                       <  3:223580/149‑1 (MQ=255)                                                         cATGTCAAAGCCCCTCAGTTTTCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCCACCAACCGTTGCATCAGCGTACGGTTTTCAACATTGTTAAATAATTAAc                                                                                       >  2:129749/1‑149 (MQ=255)                                                          aTGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAAcc                                                                                      <  8:75466/149‑1 (MQ=255)                                                                aaGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCTCAGGTCACCAATCCTTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACAtt                                                                                <  4:174748/149‑1 (MQ=255)                                                                aaGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACAtt                                                                                <  4:11300/149‑1 (MQ=255)                                                                         aGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCTGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAAtt                                                                       <  2:82174/149‑1 (MQ=255)                                                                                cACTGGCAATGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCa                                                                <  2:9707/149‑1 (MQ=255)                                                                                cACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTgccgc                                                                  >  3:155886/1‑147 (MQ=255)                                                                                 aCTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAg                                                               >  6:223882/1‑149 (MQ=255)                                                                                     gCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATc                                                           >  3:206893/1‑149 (MQ=255)                                                                                            aGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTc                                                    <  8:279022/149‑1 (MQ=255)                                                                                             gTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTACCCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCg                                                   >  1:150351/1‑149 (MQ=255)                                                                                                cAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCgtg                                                 <  6:274277/148‑1 (MQ=255)                                                                                                  ggCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCgtggtg                                              <  5:21035/149‑1 (MQ=255)                                                                                                   gCTTCCCCAATAACTGCGCAACCAGATATTACAGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCATCGTACGGTTTTCAACATAGTTCAATAAGTAACCGAGATTGAGATAATTGCCGCCAGCATCTGCGGTCGtggtgg                                             <  2:314686/149‑1 (MQ=255)                                                                                                     ttCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGt                                           >  1:206444/1‑149 (MQ=255)                                                                                                       cccccATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGtt                                          <  8:40648/148‑1 (MQ=255)                                                                                                         cccATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAg                                       <  5:191896/149‑1 (MQ=255)                                                                                                          ccATAACTGCGCCACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAgg                                      >  3:133836/1‑149 (MQ=255)                                                                                                                 tGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAGGCGTgg                                 >  8:285547/1‑147 (MQ=255)                                                                                                                     ctaCCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCCGCGTACGGTTTTCAACATAGTTCAATAAGTATCCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTTGTGGGTTTAGGCGTGGGAGCgg                           <  5:226721/147‑1 (MQ=255)                                                                                                                      aaCCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAGGCGTGGGAGCggg                          >  4:57728/1‑149 (MQ=255)                                                                                                                      aaCCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAGGCGTGGGAGCggg                          <  1:129749/149‑1 (MQ=255)                                                                                                                        ccAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGGGGGTTTAGGCGTGGGAGCGGGtt                        <  4:155886/149‑1 (MQ=255)                                                                                                                         cAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAGGCGTGGGAGCGGGTTc                       <  7:316378/149‑1 (MQ=255)                                                                                                                             tatTACCGTCTTTACTCTGATTTTGCAGGTCACCCAACCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAGGCGTGGGAGCGGGTTCAGGc                   <  4:244132/149‑1 (MQ=255)                                                                                                                             tatTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAGGCGTGGGAGCGGGTTCAGGc                   <  4:206893/149‑1 (MQ=255)                                                                                                                              atTACCGTCTTTACTCTGATTTTGCAGGTCACCCATCCGTTGCATCAGCGGACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAGGCGTGGGAGCGGGTTCAGGc                   <  4:265699/148‑1 (MQ=255)                                                                                                                                 aCCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAGGCGTGGGAGCGGGTTCAGGCGTa                <  7:285547/148‑1 (MQ=255)                                                                                                                                  ccGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAGGCGTGGGAGCGGGTTCAGGCGTAg               >  7:284641/1‑148 (MQ=255)                                                                                                                                         tACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAGGCGTGGGAGCGGGTTCAGGCGTAGGATCggg        <  8:136308/148‑1 (MQ=255)                                                                                                                                           ctctGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAGGCGTGGGAGCGGGTTCAGGCGTAGGATCGGGGtt     <  4:307055/149‑1 (MQ=255)                                                                                                                                           ctctGATTTCGCAGGACACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCTTGGTGGGTTTAGGCGTGGGAGCGGGTGCAGGCGGAGGAACGGGGtt     >  7:165463/1‑149 (MQ=255)                                                                                                                                              tGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAGGCGTGGGAGCGGGTTCAGGCGTAGGATCGGGGTTGAc  >  7:32539/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                               AACGGCATTACATCAGAGAGACGTTTATCCCCACCAAACTGGATACCGCTGTAACCCATGTCAAAGCCGCTCAGTTTGCCACTGGCAAAGGAGTCCAGGCTTCCCCCATAACTGCGCAACCAGATATTACCGTCTTTACTCTGATTTCGCAGGTCACCCATCCGTTGCATCAGCGTACGGTTTTCAACATAGTTCAATAAGTAACCGACATTGAGATAATTGCCGCCAGCATCTGCGGTCGTGGTGGGTTTAGGCGTGGGAGCGGGTTCAGGCGTAGGATCGGGGTTGAC  >  NC_000913/1253749‑1254038

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