Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	155,546	A→C	10.0%	intergenic (‑251/+32)	yjbB ← / ← ygdQ	Na+/Pi‑cotransporter/transport protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	155,546	0	A	C	10.0%	205.7 / 13.4	70	intergenic (‑251/+32)	yjbB/ygdQ	Na+/Pi‑cotransporter/transport protein
Reads supporting (aligned to +/- strand):  ref base A (32/31);  new base C (5/2);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

TCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGC  >  NC_002947/155401‑155683                                                                                                                                                   |                                                                                                                                          tCGCCAGCTTCGCATGGCTGTCATTTATTTGAAACAGATCTTGAAAGGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCATGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAg                                                                                                                                         <  8:69331/148‑1 (MQ=255) tCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGaa                                                                                                                                          <  3:316914/147‑1 (MQ=255) tCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGaa                                                                                                                                          <  6:249782/147‑1 (MQ=255) tCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAgg                                                                                                                                        <  6:307437/149‑1 (MQ=255)   gCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTTCGCAGGGc                                                                                                                                      >  8:126165/1‑149 (MQ=255)   gCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTTCGCAGGGc                                                                                                                                      >  5:82994/1‑149 (MQ=255)   gCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGc                                                                                                                                      >  8:137126/1‑149 (MQ=255)     cAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGccc                                                                                                                                    <  5:13147/149‑1 (MQ=255)         ttCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCttgt                                                                                                                                >  1:90600/1‑149 (MQ=255)           cGCATGGCAGTCATTTATTTGAATCAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGATTCGCGGGCAAGCAACCACAGTGACAACGCGCAGGAAAAAACGCACTGAAAAAAGGCATTAAAAAGGGCCCTTCGAAGGGCCCttgtt                                                                                                                               <  2:84734/148‑1 (MQ=255)            gCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTTCGCAGGGCCCttttta                                                                                                                              <  7:126165/148‑2 (MQ=255)            gCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTtgttgt                                                                                                                             >  6:225314/1‑149 (MQ=255)             cATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTtgtg                                                                                                                            <  8:306947/149‑1 (MQ=255)               tGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTTCGCAGGGCCCtttttaatgc                                                                                                                          <  6:82994/149‑6 (MQ=255)                  cAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCa                                                                                                                       <  2:90600/149‑1 (MQ=255)                       atttatttGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCa                                                                                                                       <  3:73746/144‑1 (MQ=255)                       atttatttGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCa                                                                                                                       >  4:73746/1‑144 (MQ=255)                       atttatttGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTcc                                                                                                                  >  2:33323/1‑149 (MQ=255)                          tatttGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGtgc                                                                                                                           <  3:239150/137‑1 (MQ=255)                          tatttGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGtgc                                                                                                                           >  4:239150/1‑137 (MQ=255)                           atttGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAg                                                                                                                                         <  3:19566/122‑1 (MQ=255)                           atttGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAg                                                                                                                                         >  4:19566/1‑122 (MQ=255)                           atttGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTg                                                                                                               <  5:225314/148‑1 (MQ=255)                            tttGAAACAGATCTTGAAACGAATGACTGCCCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCATgg                                                                                                              >  1:209113/1‑149 (MQ=255)                                    aGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGa                                                                                                     >  2:250235/1‑149 (MQ=255)                                    aGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGa                                                                                                     <  2:209113/149‑1 (MQ=255)                                     gATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGGATTAAAAAGGGCCCTGCGAAg                                                                                                                                         >  7:95210/1‑112 (MQ=255)                                     gATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAg                                                                                                                                         <  8:95210/112‑1 (MQ=255)                                                gAATGACTG‑CCAGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGCGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGACTTATTGACCCGGAATGtac                                                                                         >  1:553/1‑147 (MQ=255)                                                 aaTGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCt                                                                                        <  7:137126/149‑1 (MQ=255)                                                  aTGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCtt                                                                                       <  2:270104/149‑1 (MQ=255)                                                     aCTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTgc                                                                                     <  1:33323/148‑1 (MQ=255)                                                       tG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTg                                                                                                               >  8:238818/1‑120 (MQ=255)                                                       tG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTg                                                                                                               >  8:234335/1‑120 (MQ=255)                                                       tG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTg                                                                                                               <  7:238818/120‑1 (MQ=255)                                                       tG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTg                                                                                                               <  7:234335/120‑1 (MQ=255)                                                          cccGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGt                                                                                >  8:82467/1‑149 (MQ=255)                                                                      gctTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCg                                                                                                 <  4:254956/120‑1 (MQ=255)                                                                      gctTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCg                                                                                                 >  3:254956/1‑120 (MQ=255)                                                                      actTCTCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCa                                                                                                                       <  6:196058/97‑1 (MQ=255)                                                                        ttCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCAt                                                                   >  7:314307/1‑148 (MQ=255)                                                                             gggCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCt                                                             >  4:25016/1‑149 (MQ=255)                                                                              ggCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCtt                                                            >  6:51420/1‑149 (MQ=255)                                                                              ggCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGCGCCCTTCTTCTGCTCCACTACCTTGGCCTTATTGACACGGAATGACCTTGCGCAGTATCACCGGTTCATGATCCtt                                                            >  4:78804/1‑149 (MQ=255)                                                                                   gCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCtt                                                                                                                >  8:197040/1‑92 (MQ=255)                                                                                   gCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCtt                                                                                                                <  7:197040/92‑1 (MQ=255)                                                                                                       cAGGAAAAAACGCACACAAAAAAGGCATTAAAAAGGGCCCTGCCAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTg                                   <  3:201908/149‑1 (MQ=255)                                                                                                        aGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTgc                                  >  4:224899/1‑149 (MQ=255)                                                                                                        aGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTgc                                  >  1:104210/1‑149 (MQ=255)                                                                                                        aGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTgc                                  >  2:163050/1‑149 (MQ=255)                                                                                                        aGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTgc                                  <  2:267217/149‑1 (MQ=255)                                                                                                         ggAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTgcg                                 <  2:104210/149‑1 (MQ=255)                                                                                                         ggAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTgcg                                 >  5:196463/1‑149 (MQ=255)                                                                                                            aaaaaCGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCAt                              >  2:232873/1‑149 (MQ=255)                                                                                                              aaaCGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGc                            <  7:205918/149‑1 (MQ=255)                                                                                                                    acaGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGt                      >  7:203137/1‑149 (MQ=255)                                                                                                                    acaGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGt                      <  7:82467/149‑1 (MQ=255)                                                                                                                     caGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTgcgc                                >  3:39505/1‑138 (MQ=255)                                                                                                                     caGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTgcgc                                <  4:39505/138‑1 (MQ=255)                                                                                                                      aGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGtt                     >  6:145578/1‑148 (MQ=255)                                                                                                                        aaaaaaGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGAt                  <  6:196463/149‑1 (MQ=255)                                                                                                                         aaaaaGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCTGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATg                 <  4:441/149‑1 (MQ=255)                                                                                                                                cATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCa          <  3:25016/149‑1 (MQ=255)                                                                                                                                 aTTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCAc         >  8:10577/1‑149 (MQ=255)                                                                                                                                   caacaaGGGCCCTTCGCAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCg       >  8:1983/5‑149 (MQ=255)                                                                                                                                    aacaaGGGCCCTTCGCAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGc      >  4:3338/4‑149 (MQ=255)                                                                                                                                     aaaaGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACACGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGc      >  5:344139/1‑148 (MQ=255)                                                                                                                                      aaaGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGc      >  1:247625/1‑147 (MQ=255)                                                                                                                                        aGGGCCCTTCGCAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGc  >  6:335322/1‑149 (MQ=255)                                                                                                                                        aGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGc  <  5:51420/149‑1 (MQ=255)                                                                                                                                               tGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCa    >  6:94616/1‑140 (MQ=255)                                                                                                                                               tGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCa    <  5:94616/140‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                          TCGCCAGCTTCGCATGGCAGTCATTTATTTGAAACAGATCTTGAAACGAATGACTG‑CCCGTGGGGGCTGCTTCGCGGGCAAGCCACCACAGTGACAACGCGCAGGAAAAAACGCACAGAAAAAAGGCATTAAAAAGGGCCCTGCGAAGGGCCCTTGTTGTGCTGCACTTCCTTGGCCTTATTGACCCGGAATGTCCTTGCGCAGTTTCACCGGTTCATGATCCTTGCCATCCTTGCGGGCCATGGCGGTGCGCATGCGGATGTTGATGGCTTCCACCGCCAGC  >  NC_002947/155401‑155683

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