Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,835,648	C→A	7.1%	intergenic (‑90/+42)	PP_1639 ← / ← yohC	SprT family protein/inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,835,648	0	C	A	7.1%	212.5 / 4.3	68	intergenic (‑90/+42)	PP_1639/yohC	SprT family protein/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (31/31);  new base A (2/3);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

AAAAGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAAC  >  NC_002947/1835501‑1835795                                                                                                                                                    |                                                                                                                                                    aaaaGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTgcg                                                                                                                                                    >  4:293599/1‑149 (MQ=255) aaaaGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTgcg                                                                                                                                                    >  2:57804/1‑149 (MQ=255)             cTTGCTGGTAACAGGTTTCGACGCGTTGTGTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGt                                                                                                                                        <  3:224578/149‑1 (MQ=255)              ttGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGtt                                                                                                                                       >  1:53311/1‑149 (MQ=255)                gCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGGTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTg                                                                                                                                      <  7:185594/148‑1 (MQ=255)                 cTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGGGGCCGCCGCAAGGCGGCCTCTGGCTgt                                                                                                                                     <  7:200463/148‑2 (MQ=255)                  tGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATg                                                                                                                                   >  1:243169/1‑149 (MQ=255)                     tAACCGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGccc                                                                                                                                <  6:216495/149‑1 (MQ=255)                       aCAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCa                                                                                                                              <  5:60919/149‑1 (MQ=255)                       aCAGGGTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCa                                                                                                                              <  7:297104/149‑1 (MQ=255)                              tCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTg                                                                                                                        <  2:184530/148‑1 (MQ=255)                              tCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTg                                                                                                                        <  5:330425/148‑1 (MQ=255)                                        gttTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTg                                                                                                              <  7:346099/148‑1 (MQ=255)                                           tGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTGCGGTCTGTTc                                                                                                           >  6:177682/1‑148 (MQ=255)                                             aGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAg                                                                                                        <  5:210168/149‑1 (MQ=255)                                              gctgctCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGt                                                                                                       >  5:244246/1‑149 (MQ=255)                                              gctgctCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGt                                                                                                       >  8:107331/1‑149 (MQ=255)                                              gctgctCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGt                                                                                                       <  6:244246/149‑1 (MQ=255)                                              gctgctCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGt                                                                                                       >  5:107225/1‑149 (MQ=255)                                                  ctCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGgt                                                                                                   >  5:44695/1‑149 (MQ=255)                                                       gCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGGGGCCGCCGCAAGGCGGCCTCTgg                                                                                                                                         >  2:359171/1‑106 (MQ=38)                                                       gCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGGGGCCGCCGCAAGGCGGCCTCTgg                                                                                                                                         <  1:359171/106‑1 (MQ=38)                                                       gCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGa                                                                                              <  5:24616/149‑1 (MQ=255)                                                         aTGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGgt                                                                                                   >  3:237529/1‑142 (MQ=255)                                                         aTGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGgt                                                                                                   <  4:237529/142‑1 (MQ=255)                                                           ggcggggcggTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACggg                                                                                          >  5:306999/1‑149 (MQ=255)                                                            gcggggcggTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGc                                                                                         <  3:293599/149‑1 (MQ=255)                                                              ggggcggTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGccc                                                                                       >  6:57154/1‑149 (MQ=255)                                                                       ctcCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGcc                                                                                 >  8:260840/1‑146 (MQ=255)                                                                       ctcCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCa                                                                               <  6:107225/148‑1 (MQ=255)                                                                       ctcCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCa                                                                               <  2:53311/148‑1 (MQ=255)                                                                         ccTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGACGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTc                                                                            >  8:160464/1‑149 (MQ=255)                                                                          cTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTcc                                                                           >  8:321281/1‑149 (MQ=255)                                                                          cTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTcc                                                                           >  3:133830/1‑149 (MQ=255)                                                                                aGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGa                                                                      >  4:296020/1‑148 (MQ=255)                                                                                aGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGa                                                                      >  1:239922/1‑148 (MQ=255)                                                                                    aTTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTc                                                                            >  1:154426/1‑138 (MQ=255)                                                                                    aTTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTc                                                                            <  2:154426/138‑1 (MQ=255)                                                                                     ttATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCa                                                                >  4:199003/1‑149 (MQ=255)                                                                                       aTGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCAc                                                               <  1:57804/148‑1 (MQ=255)                                                                                        tGCCGCAAGCGCAGGTGGCTTGCCAAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACgg                                                             >  7:231107/1‑149 (MQ=255)                                                                                          ccGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGTGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACgg                                                             >  5:280954/1‑147 (MQ=255)                                                                                           cGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAg                                                          >  4:331733/1‑149 (MQ=255)                                                                                           cGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAg                                                          <  5:57154/149‑1 (MQ=255)                                                                                           cGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAg                                                          >  2:344873/1‑149 (MQ=255)                                                                                                      gTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGAtt                                               <  7:107331/149‑1 (MQ=255)                                                                                                        ggCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGc                                             <  8:188914/149‑1 (MQ=255)                                                                                                        ggCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGc                                             <  7:260840/149‑1 (MQ=255)                                                                                                                   gcACCGGTCAGCGCACAGCCAGATGCCGCCTTGGGGCGGCCCCTGGGTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGtacta                                  <  5:177682/148‑1 (MQ=255)                                                                                                                      ccGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTAcc                                >  6:137856/1‑148 (MQ=255)                                                                                                                       cGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCaa                              <  7:160464/149‑1 (MQ=255)                                                                                                                         gTCAGCGCACAGCCAGAGGCGGCATTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGc                            <  5:305794/149‑1 (MQ=255)                                                                                                                          tCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGcc                           <  2:239922/149‑1 (MQ=255)                                                                                                                          tCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGcc                           <  3:331733/149‑1 (MQ=255)                                                                                                                              cgcACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACt                       >  6:217905/1‑149 (MQ=255)                                                                                                                                cacaGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGc                     >  1:231986/1‑149 (MQ=255)                                                                                                                                 tcaaccAGGGGCCGCCGCAAGGCGGCCTCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGc                     >  8:109242/5‑148 (MQ=255)                                                                                                                                 tcaaccAGGGGCCGCCGCAAGGCGGCCTCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGc                     <  7:109242/144‑1 (MQ=255)                                                                                                                                        gagGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCAcc              >  6:295723/1‑148 (MQ=255)                                                                                                                                        gagGCCGCCTTGCGGCGGCCCCTGGTAGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCCGCAcc              >  3:310430/1‑148 (MQ=255)                                                                                                                                         agGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCAcc              >  6:64357/1‑147 (MQ=255)                                                                                                                                          ggCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGc                 <  2:21524/143‑1 (MQ=255)                                                                                                                                          ggCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGc                 >  1:21524/1‑143 (MQ=255)                                                                                                                                          ggCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAg           <  5:64357/149‑1 (MQ=255)                                                                                                                                           gccggcTTGCGGGGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGTTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAg           <  4:310430/148‑1 (MQ=255)                                                                                                                                           gccgccTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGcc                                                     >  3:99085/1‑106 (MQ=255)                                                                                                                                           gccgccTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGcc                                                     <  4:99085/106‑1 (MQ=255)                                                                                                                                           gccgccTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGc          >  1:306271/1‑149 (MQ=255)                                                                                                                                              gccTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTgg       <  3:199003/149‑1 (MQ=255)                                                                                                                                                cTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGgag     <  6:44695/149‑1 (MQ=255)                                                                                                                                                    cggcggCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAAc  <  1:344873/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    AAAAGGCTTCGGCTTGCTGGTAACAGGTTTCGACGCGTTGTTTGAGCTGCTCGGGCATGGCGGGGCGGTTCTCCTGACCAGGCATTATGCCGCAAGCGCAGGTGGCTTGCCGAGTCACCGGTCAGCGCACAGCCAGAGGCCGCCTTGCGGCGGCCCCTGGTTGATGCCCCAGTGTTGGTTCGGTCTGTTCTAGTTGGTGTAGACGGGCCCCACGCCCAGTCCCCAGATAATCACGGTAGCGGCCATGATTGCCACCAGTACTACCAAGCCCACTGCCAGCACCGAGCTGGAGAAC  >  NC_002947/1835501‑1835795

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