Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,436,195	T→A	7.5%	W44R (TGG→AGG)	gpI →	tail formation protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,436,195	0	T	A	7.5%	204.4 / 5.1	67	W44R (TGG→AGG)	gpI	tail formation protein
Reads supporting (aligned to +/- strand):  ref base T (31/31);  new base A (2/3);  total (33/34)
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.00e-01

GAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGG  >  NC_002947/3436054‑3436338                                                                                                                                                |                                                                                                                                                gAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCtggctggc                                                                                                                                          >  3:30684/1‑149 (MQ=255) gAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCtggctggc                                                                                                                                          >  8:321022/1‑149 (MQ=255)         tgAAATGAGCCTGTTACCGCACAACGCCACCCATTTGGAGCGCGCCCTGGAGGCGGCCGCCGTCCTTGGCCTTTACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTGTCCCGCCGAACTTCC‑TGCCCAGGCCAGCCAGGGCa                                                                                                                                   <  1:229480/149‑1 (MQ=255)           aCATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATgg                                                                                                                                <  2:80552/149‑1 (MQ=255)           aCATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATgg                                                                                                                                <  7:321022/149‑1 (MQ=255)               gAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAg                                                                                                                             <  1:334050/148‑1 (MQ=255)                     ggtACCGCACAACGCCACCCAGTGGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGCCCCGGACATTAGCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAgg                                                                                                                      <  1:182392/147‑1 (MQ=255)                      ttACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAggg                                                                                                                     >  6:189796/1‑149 (MQ=255)                        aCCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCt                                                                                                                   >  2:269261/1‑149 (MQ=255)                        aCCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCt                                                                                                                   >  7:327920/1‑149 (MQ=255)                         ccGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTg                                                                                                                  >  6:181116/1‑149 (MQ=255)                         ccGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTg                                                                                                                  <  2:184178/149‑1 (MQ=255)                            cacaACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTggg                                                                                                                                     >  5:211678/1‑127 (MQ=255)                            cacaACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTggg                                                                                                                                     <  6:211678/127‑1 (MQ=255)                                 cGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGTCATGGAAGGTCGAGGGCTGGGAggcgg                                                                                                          >  6:93603/1‑149 (MQ=255)                                 cGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggcgg                                                                                                          <  3:247645/149‑1 (MQ=255)                                  gCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggcggc                                                                                                         >  5:346036/1‑149 (MQ=255)                                  gCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAggcggc                                                                                                         <  7:93351/149‑1 (MQ=255)                                     aCCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTa                                                                                                      >  7:240956/1‑149 (MQ=255)                                       ccAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGTGCTGGGAGGCGGCTTaca                                                                                                    >  8:235029/1‑149 (MQ=255)                                        cAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAggg                                                                                                                     >  2:213787/1‑131 (MQ=255)                                        cAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAggg                                                                                                                     <  1:213787/131‑1 (MQ=255)                                                  gcgcCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAg                                                                                         >  5:185657/1‑149 (MQ=255)                                                  gcgcCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAg                                                                                         >  1:339101/1‑149 (MQ=255)                                                               ggcCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTgg                                                                                                                                      >  3:105558/1‑91 (MQ=255)                                                               ggcCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTgg                                                                                                                                      <  4:105558/91‑1 (MQ=255)                                                                  cgccgACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATcc                                                                         >  1:192239/1‑149 (MQ=255)                                                                  cgccgACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTAACACCGACGACCAGCGCCGCGCGCTGATcc                                                                         >  2:283236/1‑149 (MQ=255)                                                                   gccgACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCcg                                                                        >  4:347440/1‑149 (MQ=255)                                                                      gACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCcgc                                                                       >  5:160145/1‑147 (MQ=255)                                                                              ccTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGACCTTCCATGCCCAGGCCAGCCAGGGCA‑GGAAGTTCggcg                                                                                                                      <  2:106397/92‑4 (MQ=17)                                                                              ccTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGACCTTCCATGCCCAGGCCAGCCAGGGCA‑GGAAGTTCggcg                                                                                                                      >  1:106397/1‑89 (MQ=17)                                                                              ccTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTc                                                             >  1:126526/1‑149 (MQ=255)                                                                               cTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTcc                                                            <  6:346036/149‑1 (MQ=255)                                                                                    cccGGATAATGTCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGtt                                                        >  1:83437/1‑148 (MQ=255)                                                                                     ccGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTca                                                      <  5:181116/149‑1 (MQ=255)                                                                                     ccGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCCACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTca                                                      >  3:240550/1‑149 (MQ=255)                                                                                      cGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTcac                                                     <  3:165315/149‑1 (MQ=255)                                                                                       ggACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTCGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTcac                                                     >  7:3783/1‑148 (MQ=255)                                                                                        gACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACaa                                                   <  5:189796/149‑1 (MQ=255)                                                                                        gACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCACGCGCTGATCCGCGAGGCCATTCCGGTTcac                                                     >  6:212971/1‑147 (MQ=255)                                                                                                 cGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAgg                                           <  3:347440/148‑1 (MQ=255)                                                                                                 cGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAgg                                           <  8:240956/148‑1 (MQ=255)                                                                                                  ggggCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGACTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAggg                                          >  4:287458/1‑148 (MQ=255)                                                                                                   gggCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCa                                        <  6:211347/149‑1 (MQ=255)                                                                                                          gCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCg                                  <  2:339101/148‑1 (MQ=255)                                                                                                          gCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCGTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCgg                                 <  4:30684/149‑1 (MQ=255)                                                                                                            cGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGTCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgcg                               <  5:93603/149‑1 (MQ=255)                                                                                                            cGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgcg                               >  4:167404/1‑149 (MQ=255)                                                                                                            cGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgcg                               >  5:44564/1‑149 (MQ=255)                                                                                                             gACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTca                                                      >  1:170596/1‑125 (MQ=255)                                                                                                             gACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTca                                                      <  2:170596/125‑1 (MQ=255)                                                                                                             gACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgcgc                              <  2:126526/149‑1 (MQ=255)                                                                                                                   tCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACcga                                                                                                <  4:151277/77‑1 (MQ=255)                                                                                                                   tCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACcga                                                                                                >  3:151277/1‑77 (MQ=255)                                                                                                                     aCGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggc                      <  2:192239/149‑1 (MQ=255)                                                                                                                       gCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGgcg                               >  2:128265/1‑138 (MQ=255)                                                                                                                       gCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg                    <  6:160145/149‑1 (MQ=255)                                                                                                                       gCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg                    >  4:5491/1‑149 (MQ=255)                                                                                                                       gCTG‑CCCGCCGAACTTCC‑TGCCCTGACTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTcggcgg                    <  3:287458/149‑1 (MQ=255)                                                                                                                                  gAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCTGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAgg          <  1:269261/149‑1 (MQ=255)                                                                                                                                  gAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAgg          <  8:327920/149‑1 (MQ=255)                                                                                                                                      ttCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGGGCGATTCGGCGGGTg                 >  3:238205/1‑138 (MQ=255)                                                                                                                                      ttCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGGGCGATTCGGCGGGTg                 <  4:238205/138‑1 (MQ=255)                                                                                                                                        ccatGCCCAGGCCAGCCAGGGCA‑GGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCgcgcgc                                                                               <  6:185954/71‑1 (MQ=18)                                                                                                                                        ccatGCCCAGGCCAGCCAGGGCA‑GGAAGTTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCgcgcgc                                                                               >  5:185954/4‑74 (MQ=18)                                                                                                                                           tGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgc  >  5:117469/1‑148 (MQ=255)                                                                                                                                             cccTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg  <  2:83437/147‑1 (MQ=255)                                                                                                                                             cccTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCgg  <  7:235029/147‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                GAAGGTGGTGACATGAGCCTGTTACCGCACAACGCCACCCAGTTGGAGCGCGCCCTGGAGGCGGCCGCCGACCTTGGCCTTGACCCGGACATTATCCGGGGCGTGGCCGACTCGTCACGCTG‑CCCGCCGAACTTCC‑TGCCCTGGCTGGCCTGGGCATGGAAGGTCGAGGGCTGGGAGGCGGCTTACACCGACGACCAGCGCCGCGCGCTGATCCGCGAGGCCATTCCGGTTCACAAGACCAAGGGCACGGTCGGCGCGATTCGGCGGGTGCTCAAGGCGGTTCGG  >  NC_002947/3436054‑3436338

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