Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,003,952	T→C	5.5%	*172Q (TAA→CAA)  T2T (ACT→ACC)	ycbV → ycbF →	putative fimbrial‑like adhesin protein putative periplasmic pilin chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,003,952	0	T	C	5.5%	223.8 / 2.1	73	*172Q (TAA→CAA)  T2T (ACT→ACC)	ycbV ycbF	putative fimbrial‑like adhesin protein putative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base T (36/33);  new base C (2/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTCCGCTTAACCAGGCGAGTAAAGACTACAGCCTACATAGCGGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATAC  >  NC_000913/1003808‑1004098                                                                                                                                                 |                                                                                                                                                   aTTCCGCTTAACCAGGCGAGTAAAGACTACAGCCTACATAGCGGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATa                                                                                                                                                >  1:144945/1‑149 (MQ=255)  ttCCGCTTAACCAGGCGAGTAAAGACTACAGCCTACATAGCGGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTTCTTCATTATGACTAATAc                                                                                                                                               >  4:195078/1‑149 (MQ=255)     cGCTTAACCAGGCGAGTAAAGACTACAGCCTACATAGCGGTAATGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTg                                                                                                                                            <  5:293259/149‑1 (MQ=255)      gCTTAACCAGGCGAGTAAAGACTACAGCCTACATAGCGGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTgg                                                                                                                                           <  6:214772/149‑1 (MQ=255)      gCTTAACCAGGCGAGTAAAGACTACAGCCTACATAGCGGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTg                                                                                                                                            >  2:244046/1‑148 (MQ=255)            ccAGGCGAGTATAGACTACATCCTACATAGCTGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTTTCAATATCGATGTGCAGTCGTGTAAAGTGAATGCCAGCGCGACATTTGTACTTAATTATGACTATTACGTGGAACAGa                                                                                                                                     <  2:56498/149‑1 (MQ=255)            ccAGGCGAGTAAAGACTACAGCCTACATAGCGGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGa                                                                                                                                     <  1:244046/149‑1 (MQ=255)             cAGGCGAGTAAAGACTACAGCCTACATAGCGGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGAt                                                                                                                                    >  5:42849/1‑149 (MQ=255)             cAGGCGAGTAAAGACTACAGCCTACATAGCGGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGAt                                                                                                                                    >  3:108667/1‑149 (MQ=255)              aGGCGAGTAAAGACTACAGCCTACATAGCGGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCAATATGACTAATACGTGGAACAGAtt                                                                                                                                   >  8:289667/1‑149 (MQ=255)                gCGAGTAAAGACTACAGCCTACATAGCGGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTgg                                                                                                                                 <  2:144945/149‑1 (MQ=255)                 cGAGTAAAGACTACAGCCTACATAGCGGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGgc                                                                                                                                <  5:72435/149‑1 (MQ=255)                    gTAAAGACTACAGCCTACATAGCGGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTgg                                                                                                                                 >  7:213288/1‑145 (MQ=255)                    gTAAAGACTACAGCCTACATAGCGGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTgg                                                                                                                                 <  8:213288/145‑1 (MQ=255)                     tAAAGACTACAGCCTACATAGCGGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGctc                                                                                                                            >  7:126908/1‑149 (MQ=255)                     tAAAGACTACAGCCTACATAGCGGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGctc                                                                                                                            >  3:43948/1‑149 (MQ=255)                     tAAAGACTACAGCCTACATAGCGGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGctc                                                                                                                            >  7:80994/1‑149 (MQ=255)                                gCCTACATAGCGGTAAGGTGCCATTATCTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGc                                                                                                                 >  3:253766/1‑149 (MQ=255)                                 ccTACATAGCGGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGcc                                                                                                                >  2:171704/1‑149 (MQ=255)                                   tACATAGCGGTAAGGTGCCATTATCTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATCATGACTAATACGTGGAACAGATTGGCGCTGTTGATTTTCGCCGt                                                                                                              >  3:164212/1‑149 (MQ=255)                                   tACATAGCGGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGt                                                                                                              >  6:194450/1‑149 (MQ=255)                                   tACATAGCGGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGt                                                                                                              >  8:23992/1‑149 (MQ=255)                                         gCGGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTAt                                                                                                         <  7:23992/148‑1 (MQ=255)                                         gCGGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATc                                                                                                        >  1:306627/1‑149 (MQ=255)                                             tAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCgctg                                                                                                    <  6:327916/149‑1 (MQ=255)                                                ggTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGctgctg                                                                                                 <  3:195078/149‑1 (MQ=255)                                                ggTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGctgctg                                                                                                 >  5:225070/1‑149 (MQ=255)                                                 gTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCtgg                                                                                                <  3:120762/149‑1 (MQ=255)                                                 gTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCtgg                                                                                                <  2:214288/149‑1 (MQ=255)                                                      aTTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCg                                                                                           >  8:16785/1‑149 (MQ=255)                                                       ttAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAATTGAATGCCAGCGCGAAATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTAGCCGTTTTATCGCTGCTGGTGGctg                                                                                          >  1:217082/1‑147 (MQ=255)                                                       ttAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCgg                                                                                          >  3:54276/1‑149 (MQ=255)                                                       ttAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCgg                                                                                          <  4:108667/149‑1 (MQ=255)                                                          aCTTTTTATGCGCAATTGCGACCTGTCAATATCGATGTGCATTCGGGTAAAGTGAATTAAAGCTCGACATTTGTACTTAATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGt                                                                                        <  2:75703/148‑1 (MQ=255)                                                          aCTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGt                                                                                        <  4:43948/148‑1 (MQ=255)                                                            tttttATGCGCAATTGCGCCCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGaa                                                                                     >  3:78921/1‑149 (MQ=255)                                                              tttATGCGCAATTGCTAACTGTCAATAGCGATTTGCAGTCGGGTAAAGTGAATTACAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTCGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAAt                                                                                    <  2:217082/148‑1 (MQ=255)                                                                  tGCGCAATTGCGACCTGTCAATAGCGATGTTCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTgca                                                                                <  6:250140/148‑1 (MQ=255)                                                                  tGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTgcag                                                                               >  4:251983/1‑149 (MQ=255)                                                                    cgcAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTgcagg                                                                              <  5:300016/148‑1 (MQ=255)                                                                    cgcAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTCGTGGCGGGTGAATTgcagg                                                                              >  5:263138/1‑148 (MQ=255)                                                                      cAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTgcaggcag                                                                           <  4:253766/149‑1 (MQ=255)                                                                              aCCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAgtggtg                                                                   >  5:114082/1‑149 (MQ=255)                                                                                 tGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTc                                                                >  3:283078/1‑149 (MQ=255)                                                                                 tGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACCAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTc                                                                >  1:136795/1‑149 (MQ=255)                                                                                 tGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACCAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTc                                                                >  1:141493/1‑149 (MQ=255)                                                                                   tCAAAAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTAATTCAATATGACTAAAACGTGGAACAGATAGGCGCTCTTGATTTTCGCCGTTTTATAGCTGATTGTGGCGGATGAATTACAGGCAGGAGTGGTGGTctg                                                              >  3:54796/1‑147 (MQ=255)                                                                                          cGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAAcg                                                       >  2:171122/1‑149 (MQ=255)                                                                                            aTGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTgg                                                                                                                                 <  5:257018/73‑1 (MQ=255)                                                                                            aTGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTgg                                                                                                                                 >  6:257018/1‑73 (MQ=255)                                                                                               tgCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGAt                                                   <  6:114082/148‑1 (MQ=255)                                                                                                        gTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTcc                                         <  7:16785/149‑1 (MQ=255)                                                                                                        gTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATAAGTGGAACAGATTGGCGCTCTTGATTTTCGCAGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGATGTCGGCGGAACGCGATTAATCTTTcc                                         >  8:128314/1‑149 (MQ=255)                                                                                                         tAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCg                                        <  6:42849/149‑1 (MQ=255)                                                                                                           aaGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTTGTGGTCGGCGGAACGCGATTTATCTTTCCGGc                                      <  2:306627/149‑1 (MQ=255)                                                                                                            aGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGca                                     <  4:78921/149‑1 (MQ=255)                                                                                                                     cAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGaa                            <  4:283078/149‑1 (MQ=255)                                                                                                                       gcgcGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGttt                                                                                                            >  5:47230/1‑67 (MQ=255)                                                                                                                       gcgcGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGttt                                                                                                            <  6:47230/67‑1 (MQ=255)                                                                                                                          cgACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGat                       <  6:263138/149‑1 (MQ=255)                                                                                                                          cgACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGGCGGCGGAACGCGATTTATCTTTCCAGCAGACAGAGAATCGat                       <  8:126908/149‑1 (MQ=255)                                                                                                                           gACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGata                      >  3:251290/1‑149 (MQ=255)                                                                                                                            aCATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGatat                     >  5:310980/1‑149 (MQ=255)                                                                                                                               tttGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTa                  >  2:152008/1‑149 (MQ=255)                                                                                                                               tttGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTAGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTa                  >  5:141008/1‑149 (MQ=255)                                                                                                                                     cTTCATTCTGACCAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTa             <  2:141493/148‑1 (MQ=255)                                                                                                                                     cTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTAc            >  7:289874/1‑149 (MQ=255)                                                                                                                                     cTTCATTATGACCAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTa             <  2:136795/148‑1 (MQ=255)                                                                                                                                        cATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACt           <  8:80994/147‑1 (MQ=255)                                                                                                                                         atcaTGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGa         <  4:164212/145‑1 (MQ=255)                                                                                                                                         aTTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGa         <  4:54276/148‑1 (MQ=255)                                                                                                                                              gACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATa   >  7:128205/1‑149 (MQ=255)                                                                                                                                               aCTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCGGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATa   <  5:194450/148‑1 (MQ=255)                                                                                                                                                 tAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATAc  >  7:239737/1‑147 (MQ=255)                                                                                                                                                 |                                                                                                                                                   ATTCCGCTTAACCAGGCGAGTAAAGACTACAGCCTACATAGCGGTAAGGTGCCATTAACTTTTTATGCGCAATTGCGACCTGTCAATAGCGATGTGCAGTCGGGTAAAGTGAATGCCAGCGCGACATTTGTACTTCATTATGACTAATACGTGGAACAGATTGGCGCTCTTGATTTTCGCCGTTTTATCGCTGCTGGTGGCGGGTGAATTGCAGGCAGGAGTGGTGGTCGGCGGAACGCGATTTATCTTTCCGGCAGACAGAGAATCGATATCTATTTTACTGACTAATAC  >  NC_000913/1003808‑1004098

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