Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	626,569	C→A	6.4%	L208L (CTG→CTT)	dedA ←	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	626,569	0	C	A	6.4%	248.2 / 6.0	78	L208L (CTG→CTT)	dedA	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (34/39);  new base A (3/2);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.05e-01

CCAGCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACCCAAAGCAG  >  NC_002947/626425‑626712                                                                                                                                                  |                                                                                                                                                ccAGCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGgc                                                                                                                                             <  3:268714/149‑1 (MQ=255)  cAGCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGgcg                                                                                                                                            >  4:55581/1‑149 (MQ=255)   aGCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGAGCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGgcgc                                                                                                                                           <  7:221424/149‑1 (MQ=255)   aGCAATGGCAGGCCGTACGCGAGCGG‑CTGCCAATGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGCGCCTTGGCCGTACGGCACAGgcgc                                                                                                                                           >  4:178080/1‑149 (MQ=37)    gCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTCGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGcc                                                                                                                                          >  3:151400/1‑149 (MQ=255)     ctttGGCA‑GCCGCTCGCGTAGGGGCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCGAGAGTACGCTTGCTGCGCGTGGCTCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTTCGGCCCAGGCGcc                                                                                                                                          <  4:151400/146‑1 (MQ=255)     cATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGccc                                                                                                                                         <  2:223673/149‑1 (MQ=255)     cATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGcc                                                                                                                                          >  8:188871/1‑148 (MQ=255)     cATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGCCGTACGGCCAAGgctcac                                                                                                                                         <  4:22423/149‑5 (MQ=255)          gCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACgc                                                                                                                                     <  4:143090/148‑1 (MQ=255)              gCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACgcagc                                                                                                                                  >  5:160726/1‑148 (MQ=255)                 gCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGc                                                                                                                              <  3:293430/149‑1 (MQ=255)                      cgTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCgaga                                                                                                                         <  5:71012/149‑1 (MQ=255)                      cgTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCgaga                                                                                                                         <  6:12423/149‑1 (MQ=255)                       gTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCgaga                                                                                                                         >  2:270838/1‑148 (MQ=255)                       gTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGaa                                                                                                                        >  5:28328/1‑149 (MQ=255)                       gTACGGCCAGCCATAGCTGAGGCGTGAACGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAAAGACCACATGCCGTTCCTTCTTTACCAGTTCACTGACCCTAGGCCGTACGGCCCACGCGCCCACGCAGCAAGCCcaga                                                                                                                         >  4:281016/1‑148 (MQ=255)                        tACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCga                                                                                                                           >  6:228241/1‑145 (MQ=255)                        tACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCga                                                                                                                           <  5:228241/145‑1 (MQ=255)                          cGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAatc                                                                                                                      <  1:140822/148‑1 (MQ=255)                             ccTGCCATTGCTGTGGCGTGATCGGGTAACNCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCg                                                                                                                  >  1:23709/1‑149 (MQ=255)                             ccTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCg                                                                                                                  >  5:167045/1‑149 (MQ=255)                             ccTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCg                                                                                                                  <  7:83408/149‑1 (MQ=255)                             ccTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCg                                                                                                                  >  5:175237/1‑149 (MQ=255)                              cTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCTTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCgg                                                                                                                 <  1:68196/149‑1 (MQ=255)                              cTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCgg                                                                                                                 <  1:270838/149‑1 (MQ=255)                              cTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCgg                                                                                                                 <  6:160726/149‑1 (MQ=255)                               tGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCgg                                                                                                                 >  2:227820/1‑148 (MQ=255)                                 ccATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAatcatcaga                                                                                                                >  5:271240/1‑144 (MQ=255)                                 ccATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAatcat                                                                                                                    <  6:271240/143‑1 (MQ=255)                                          tgGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAGTCATCGGCACCAGCGAc                                                                                                       <  3:286599/147‑1 (MQ=255)                                            gCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCTGCACCAGCGACAGga                                                                                                   >  6:214448/1‑149 (MQ=255)                                            gCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGga                                                                                                   >  6:138388/1‑149 (MQ=255)                                               tGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCAc                                                                                                              >  1:223163/1‑135 (MQ=255)                                               tGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCAc                                                                                                              <  2:223163/135‑1 (MQ=255)                                                  tCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGgatgatga                                                                                             >  2:272624/1‑149 (MQ=255)                                                       tAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGccc                                                                                        >  4:139219/1‑149 (MQ=255)                                                         aCGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCccac                                                                                      >  7:120474/1‑149 (MQ=255)                                                           gcgcCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCccacca                                                                                    <  8:120474/149‑1 (MQ=255)                                                            cgcCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAg                                                                                                                               >  4:286901/1‑105 (MQ=255)                                                            cgcCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAg                                                                                                                               <  3:286901/105‑1 (MQ=255)                                                            cgcCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCcccccac                                                                                   >  7:119446/1‑149 (MQ=255)                                                                      cGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGAc                                                                         <  6:28328/149‑1 (MQ=255)                                                                       gCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGa                                                                          >  4:206196/1‑147 (MQ=255)                                                                        cTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAg                                                                       <  3:55581/149‑1 (MQ=255)                                                                            cgccgcCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGtgc                                                                   <  1:227820/149‑1 (MQ=255)                                                                             gccgccAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGtgct                                                                  >  6:124639/1‑149 (MQ=255)                                                                               cgccAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGtgctgc                                                                <  1:272624/149‑1 (MQ=255)                                                                                 ccAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCtt                                                              >  6:16773/1‑149 (MQ=255)                                                                                 ccAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCtt                                                              >  8:278350/1‑149 (MQ=255)                                                                                  cAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTtg                                                             <  5:138388/149‑1 (MQ=255)                                                                                  cAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTtg                                                             <  5:214448/149‑1 (MQ=255)                                                                                    ggCGCTGAATGACCTCATGGCGTTCGTTCTTTACCTGTTCAGTGAGCGTTGGCCGTACGGCCCAGGCGCCCACGCAGGAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGAGAGGTGCTGCTtgat                                                           <  2:23709/149‑1 (MQ=255)                                                                                       gcTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGaa                                                        >  7:173711/1‑149 (MQ=255)                                                                                         tGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATg                                                      <  3:206196/149‑1 (MQ=255)                                                                                                ccACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGGGCCGACGCAGCAAGCCGAGAATCATCGGCACCTGCGACAGGATGATGATGCCCACCAGCATCAGCGACAGGGGCTGCTTGATGAATGGCACGtt                                               <  3:281016/149‑1 (MQ=255)                                                                                                  acaTGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGc                                             <  8:119446/149‑1 (MQ=255)                                                                                                   caTGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGcc                                            <  3:139219/149‑1 (MQ=255)                                                                                                                   tACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGccc                               <  6:175237/146‑1 (MQ=255)                                                                                                                   tACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAg                            >  5:278726/1‑149 (MQ=255)                                                                                                                       tGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTg                                              <  6:239861/127‑1 (MQ=255)                                                                                                                       tGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTg                                              >  5:239861/1‑127 (MQ=255)                                                                                                                         ttCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGAc                       <  6:167045/148‑1 (MQ=255)                                                                                                                           cAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCa                     <  7:188871/148‑1 (MQ=255)                                                                                                                            tgggcgCCTGGGCCGTACGGCCAAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAgg                   >  7:35597/6‑149 (MQ=255)                                                                                                                             gTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAg                                  >  1:283731/1‑133 (MQ=255)                                                                                                                             gTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAg                                  <  2:283731/133‑1 (MQ=255)                                                                                                                               gAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGgccg                >  8:33785/1‑149 (MQ=255)                                                                                                                                 gCCTGGGCCGTACGGCCAAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGgccgcc              >  3:259190/1‑149 (MQ=255)                                                                                                                                 gCCTGGGCCGTACGGCCAAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGgccgcc              >  2:27260/1‑149 (MQ=255)                                                                                                                                   cTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGaa                                                        >  5:82191/1‑105 (MQ=255)                                                                                                                                   cTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGaa                                                        <  6:82191/105‑1 (MQ=255)                                                                                                                                   cTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGa            <  5:124639/149‑1 (MQ=255)                                                                                                                                    ttGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGAc           >  1:308674/1‑149 (MQ=255)                                                                                                                                      ggCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGAcc          <  2:235488/148‑1 (MQ=255)                                                                                                                                            aCGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACCCAAagca   >  2:50658/1‑149 (MQ=255)                                                                                                                                             cGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACCCAAagca   <  7:278350/148‑1 (MQ=255)                                                                                                                                             cGGCCAAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACCCAAagcag  <  8:35597/149‑1 (MQ=255)                                                                                                                                                 ccAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAgg                   <  6:56991/128‑1 (MQ=255)                                                                                                                                                 ccAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAgg                   >  5:56991/1‑128 (MQ=255)                                                                                                                                                  |                                                                                                                                                CCAGCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACCCAAAGCAG  >  NC_002947/626425‑626712

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