Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,030,818	A→T	8.2%	K28* (AAG→TAG) ‡	mgtA →	ATP‑dependent magnesium transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,030,818	0	A	T	8.2%	176.3 / 7.3	61	K28* (AAG→TAG) ‡	mgtA	ATP‑dependent magnesium transporter
Reads supporting (aligned to +/- strand):  ref base A (30/26);  new base T (2/3);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01

GACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCAG‑CGAACAGGGCTTGACCGAGCACGAAGCCGCCAAGCGCCTGGTGCGCGATGGGGCCAACC  >  NC_002947/3030671‑3030965                                                                                                                                                    |                                                                                                                                                       gacagacaACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAg                                                                                                                                                      >  2:188517/1‑150 (MQ=255) gacagacaACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGATTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAg                                                                                                                                                      >  1:1093723/1‑150 (MQ=255)  acagacaACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAgg                                                                                                                                                     >  2:1073048/1‑150 (MQ=255)        aaCTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAAc                                                                                                                                               >  2:915128/1‑150 (MQ=255)        aaCTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAAc                                                                                                                                               >  2:465849/1‑150 (MQ=255)          cTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCTTGGCCGatgag                                                                                                                                             <  2:1125848/150‑5 (MQ=255)             tCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCgag                                                                                                                                          >  2:806579/1‑150 (MQ=255)                 aaaGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGAccc                                                                                                                                      >  1:540227/1‑150 (MQ=255)                   aGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGc                                                                                                                                    <  1:206514/150‑1 (MQ=255)                       ccAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGt                                                                                                                                >  1:386958/1‑150 (MQ=255)                          gCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCTTGGCCGatga                                                                                                                                              >  2:165259/1‑130 (MQ=255)                          gCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCTTGGCCGatga                                                                                                                                              <  1:165259/133‑4 (MQ=255)                             ttattCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCAccc                                                                                                                          >  2:408054/1‑150 (MQ=255)                                ttCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCcgcg                                                                                                                       >  2:1126905/1‑150 (MQ=255)                                   aGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGcgccg                                                                                                                    >  1:106566/1‑150 (MQ=255)                                             cATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAgc                                                                                                          <  1:617087/150‑1 (MQ=255)                                              aTGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAgcg                                                                                                         >  1:1065780/1‑150 (MQ=255)                                               tGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAgcgc                                                                                                        <  2:1065780/150‑1 (MQ=255)                                                   ttCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCaa                                                                                                    >  2:935861/1‑150 (MQ=255)                                                   ttCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCaa                                                                                                    >  1:324803/1‑150 (MQ=255)                                                       ttCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGc                                                                                                <  1:1003665/150‑1 (MQ=255)                                                       ttCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGc                                                                                                >  1:1069810/1‑150 (MQ=255)                                                        tCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCTTGGCCGatga                                                                                                                                              >  1:118325/1‑100 (MQ=255)                                                        tCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCTTGGCCGatga                                                                                                                                              <  2:118325/103‑4 (MQ=255)                                                                    tGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCAccc                                                                                   >  1:692154/1‑150 (MQ=255)                                                                      aaaCGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTc                                                                                 <  1:216763/150‑1 (MQ=255)                                                                          gTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCa                                                                             <  2:540227/150‑1 (MQ=255)                                                                             tgtgtgTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCgg                                                                                         <  1:657365/135‑1 (MQ=255)                                                                             tgtgtgTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCgg                                                                                         >  2:657365/1‑135 (MQ=255)                                                                                 tgtCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGa                                                                      >  1:174130/1‑150 (MQ=255)                                                                                  gtCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑t                                                                    <  1:408054/150‑1 (MQ=255)                                                                                      tGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCa                                                               <  1:188517/150‑1 (MQ=255)                                                                                      tGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCa                                                               <  1:915128/150‑1 (MQ=255)                                                                                      tGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCa                                                               <  2:106566/150‑1 (MQ=255)                                                                                           ccgcAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCAG‑CGaa                                                         <  1:935861/150‑1 (MQ=255)                                                                                                gCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCAG‑CGAACAggg                                                    <  1:320243/150‑1 (MQ=255)                                                                                                gCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGAATGTGCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCAG‑CGAACAggg                                                    <  1:465849/150‑1 (MQ=255)                                                                                                     aGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCAG‑CGAACAGGGCTTGa                                               <  2:386958/150‑1 (MQ=255)                                                                                                         gCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCAG‑CGAACAGGGCTTGACCGa                                           <  1:1073048/150‑1 (MQ=255)                                                                                                            caGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCAG‑CGAACAGGGCTTGACCGAGCa                                        <  1:174570/150‑1 (MQ=255)                                                                                                             aGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGTCGGGCACCCTCGCCAACCGCTA‑T‑GTTTG‑CGAACAGGTCTTGACCGACCAc                                       >  2:790170/1‑150 (MQ=255)                                                                                                               cAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCAG‑CGAACAGGGCTTGACCGAGCACGa                                     <  1:189363/150‑1 (MQ=255)                                                                                                                    tGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCAG‑CGAACAGGGCTTGACCGAGCACGAAgcct                                <  2:1093723/150‑2 (MQ=255)                                                                                                                        tGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCTT‑T‑GCCAG‑CGAACAGGGCTTGACCGAGCACGAAGCCGCCaa                            >  1:851422/1‑150 (MQ=255)                                                                                                                        tGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCAG‑CGAACAGGGCTTGACCGAGCACGAAGCCGCCaa                            >  1:388152/1‑150 (MQ=255)                                                                                                                        tGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCAG‑CGAACAGGGCTTGACCGAGCACGAAGCCGCCaa                            >  2:792345/1‑150 (MQ=255)                                                                                                                        tGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCAG‑CGAACAGGGCTTGACCGAGCACGAAGCCGCCaa                            <  2:324803/150‑1 (MQ=255)                                                                                                                        tGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCAG‑CGAACAGGGCTTGACCGAGCACGAAGCCGCCaa                            >  2:1013312/1‑150 (MQ=255)                                                                                                                              aTACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCg                                                                                          <  2:681281/85‑1 (MQ=255)                                                                                                                              aTACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCg                                                                                          >  1:681281/1‑85 (MQ=255)                                                                                                                               tcccGTAACACCTCATCGGCCACGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGAGGTCACCCTCGCCAACCTCGA‑T‑GCCAG‑CGAACAGGGCTTGACCGAGCACGAAGCCGCCAAGCGCCTg                     <  1:790170/148‑1 (MQ=255)                                                                                                                                 ccttAACGCCTCACCGGACAGGGCCGAACGCGTGCTTCTCCTAGCCACCTTCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCAG‑CGAACAGGGCTTGACCGAGCACGAAGCCGCCAAGCGCCTGGt                   <  2:851422/147‑1 (MQ=255)                                                                                                                                   gTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCAG‑CGAACAGGGCTTGACCGAGCACGAAGCCGCCAAGCGCCTGGTgc                 <  2:692154/150‑1 (MQ=255)                                                                                                                                       cacCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCGCGA‑T‑GCCAG‑CGAACAGGGCTTGACCGAGCACAAAGCCGCCAAGCGCCTGGTGCGCGa             <  1:150781/150‑1 (MQ=255)                                                                                                                                        acCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGCT‑ACCCCCG‑‑GACTTAAACTGGCCGGACAAACAAGAGGTGACCGAGCACGAAGCCGCCAAGCGCCTGGTGCGCGAt            >  1:130625/1‑150 (MQ=255)                                                                                                                                         ccTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCGG‑CGAACAGGGCTTGACCGAGCACGAAGCCGCCAAGCGCCTGGTGCGCGATg           >  1:910521/1‑150 (MQ=255)                                                                                                                                          cTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCAG‑CGAACAGGGCTTGACCGAGCACGAAGCCGCCAAGCGCCTGGTGCGCGATgg          >  2:839778/1‑150 (MQ=255)                                                                                                                                           tCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCAG‑CGAACAGGGCTTGACCGAGCACGAAGCCGCCAAGCGCCTGGTGCGCGATggg         >  2:272835/1‑150 (MQ=255)                                                                                                                                                ggCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCAG‑CGAACAGGGCTTGACCGAGCACGAAGCCGCCAAGCGCCTGGTGCGCGATGGGGCCaa    >  1:91633/1‑150 (MQ=255)                                                                                                                                                  ccAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCAG‑CGAACAGGGCTTGACCGAGCACGAAGCCGCCAAGCGCCTGGTGCGCGATGGGGCCAAcc  >  1:97720/1‑150 (MQ=255)                                                                                                                                                  ccAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCAG‑CGAACAGGGCTTGACCGAGCACGAAGCCGCCAAGCGCCTGGTGCGCGATGGGGCCAAcc  <  1:1126905/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                       GACAGACAACTGTCGGAAAGAACCAGCATTATTCAGCGACCCGACATGCCTTCGTTCACCGCCTGAGTGAAACGTCTGTGTGTCGTGCCGCCGCAGCGGCAGGAGCACAGCAACATGACTGTCAAATACCGTAACACCTCATCGGCCAAGGCCGAACGCGAGACCCGCCTGTCCACCCGCGCCGCCCGTGAAGCGCGCAATGGCCTGGCGGTCACCCTCGCCAACCTCGA‑T‑GCCAG‑CGAACAGGGCTTGACCGAGCACGAAGCCGCCAAGCGCCTGGTGCGCGATGGGGCCAACC  >  NC_002947/3030671‑3030965

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