Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,628,430	G→T	7.8%	G17C (GGC→TGC)	clpP →	ATP‑dependent Clp protease proteolytic subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,628,430	0	G	T	7.8%	198.8 / 5.6	65	G17C (GGC→TGC)	clpP	ATP‑dependent Clp protease proteolytic subunit
Reads supporting (aligned to +/- strand):  ref base G (33/26);  new base T (2/3);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

TGCTGAGTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCAACTGCTG  >  NC_002947/2628281‑2628576                                                                                                                                                      |                                                                                                                                                   tgctgAGTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAg                                                                                                                                                    >  1:1083003/1‑150 (MQ=255)   ctgAGTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAggc                                                                                                                                                  >  1:1184795/1‑150 (MQ=255)         aaCGGTCTCCTGCAGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTg                                                                                                                                            >  2:562089/1‑150 (MQ=255)            ggTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCCTGCGGCCTggat                                                                                                                                         <  1:1268297/150‑3 (MQ=255)              tcttctCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTccc                                                                                                                                       >  1:1307832/1‑150 (MQ=255)               cttctCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCg                                                                                                                                      >  1:261026/1‑150 (MQ=255)                     tcGTAGGAAACCTCCACAAGCCAGCCTTCGAGTTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGtt                                                                                                                                >  2:444837/1‑150 (MQ=255)                     tcGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGtt                                                                                                                                <  1:1351920/150‑1 (MQ=255)                       gTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTAt                                                                                                                              >  1:829079/1‑150 (MQ=255)                        tAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATc                                                                                                                             <  2:1307832/150‑1 (MQ=255)                         aGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCg                                                                                                                            <  1:418586/150‑1 (MQ=255)                           gAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAg                                                                                                                          <  2:829079/150‑1 (MQ=255)                            aaaCCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGCCAACCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAga                                                                                                                         >  1:14469/1‑149 (MQ=255)                               ccTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGt                                                                                                                      <  2:1104940/150‑1 (MQ=255)                                cTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTc                                                                                                                     >  1:1101359/1‑150 (MQ=255)                                   cacaAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGc                                                                                                                  >  1:538636/1‑150 (MQ=255)                                        gccagccTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTg                                                                                                             >  1:428447/1‑150 (MQ=255)                                                     atGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTa                                                                                                >  2:950754/2‑150 (MQ=255)                                                      tGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCGGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTAc                                                                                               >  1:219569/1‑150 (MQ=255)                                                          ttGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACa                                                                                           >  1:650588/1‑150 (MQ=255)                                                          ttGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCACCATGGGTGATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACa                                                                                           >  2:126481/1‑150 (MQ=255)                                                                gTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGCCGAACGCGCTTACGACATCTatt                                                                                     >  1:650403/1‑148 (MQ=255)                                                                  aTTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAg                                                                                                                                                    <  2:507234/85‑1 (MQ=255)                                                                  aTTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAg                                                                                                                                                    >  1:507234/1‑85 (MQ=255)                                                                   ttCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTcgc                                                                                  <  2:650588/150‑1 (MQ=255)                                                                   ttCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTcgc                                                                                  >  1:1078716/1‑150 (MQ=255)                                                                     cAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCgcgc                                                                                <  2:1403894/150‑1 (MQ=255)                                                                       aGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTAATCGCGCCt                                                                              <  2:1078716/150‑1 (MQ=255)                                                                         ccATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCCTGCGGCCTgga                                                                                                                                          <  2:919575/88‑2 (MQ=255)                                                                         ccATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCCTGCGGCCTgga                                                                                                                                          >  1:919575/1‑87 (MQ=255)                                                                         ccATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGt                                                                            >  1:475200/1‑150 (MQ=255)                                                                           aTGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCCTGCGGCCTg                                                                                                                                            <  2:187836/84‑1 (MQ=255)                                                                           aTGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCCTGCGGCCTg                                                                                                                                            >  1:187836/1‑84 (MQ=255)                                                                            tGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAggcgg                                                                                                                                                <  2:386834/79‑1 (MQ=255)                                                                            tGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAggcgc                                                                                                                                                >  1:386834/1‑78 (MQ=255)                                                                            tGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGa                                                                         >  2:1512/1‑150 (MQ=255)                                                                             gACTATTTGGGAGTGAGCGCAGGACATGTCACCAAATCCGTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGaa                                                                        <  2:824465/150‑1 (MQ=255)                                                                                tATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGaa                                                                        >  1:1058262/1‑147 (MQ=255)                                                                                tATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGaa                                                                        <  2:1058262/147‑1 (MQ=255)                                                                                tATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGa                                                                     >  2:572208/1‑150 (MQ=255)                                                                                          tGAGCGCGGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATc                                                           <  2:219569/150‑1 (MQ=255)                                                                                              cgcAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGagatcgg                                                       >  1:1379886/1‑144 (MQ=255)                                                                                               gcAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCt                                                      <  2:428447/150‑1 (MQ=255)                                                                                                cAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTg                                                     <  2:1184795/150‑1 (MQ=255)                                                                                                   gACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGtt                                                  <  2:538636/150‑1 (MQ=255)                                                                                                     cATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTgg                                                >  1:1069105/1‑150 (MQ=255)                                                                                                          cccGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCgg                                           <  1:1512/150‑1 (MQ=255)                                                                                                           ccGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATc                                                                                                                             >  2:1405190/1‑67 (MQ=255)                                                                                                           ccGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATc                                                                                                                             <  1:1405190/67‑1 (MQ=255)                                                                                                            cGCAATTCTTATATTCAGCAGAGCTCTGGCATTCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAACAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTa                                         >  1:1327925/1‑150 (MQ=255)                                                                                                             gCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTag                                        >  2:1289162/1‑150 (MQ=255)                                                                                                             ccAATTCTCATATTCAGCCGAGACCTGCCATCCAGGCCGCACCCGGCCTGGTGCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACAGCGACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTag                                        <  1:562089/149‑1 (MQ=255)                                                                                                              cAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTaga                                       <  1:572208/150‑1 (MQ=255)                                                                                                                       tatTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACAt                              <  2:291699/150‑1 (MQ=255)                                                                                                                         ttCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGAAGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATgg                            <  2:261026/150‑1 (MQ=255)                                                                                                                          tCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGc                           <  1:1289162/150‑1 (MQ=255)                                                                                                                              gatctctctGACATCCAGGCCGCAGGCGGCCTGGTTTCGATTGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAAc                       <  2:1182936/145‑1 (MQ=255)                                                                                                                               agagCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAAcc                      >  2:74655/1‑150 (MQ=255)                                                                                                                                 cgCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATAAAGCCGTTCTCCCGTGGCGAACGCGCTTACGACATCTAGTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCtc                    <  1:87572/149‑1 (MQ=255)                                                                                                                                   ctctGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCtcgt                  >  2:411878/1‑150 (MQ=255)                                                                                                                                    tctGACATCCAGGCCGCAGGCGGCCTGGTCCCGCCGTTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAAccagatc                 >  1:1182936/1‑145 (MQ=255)                                                                                                                                         cATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGAAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTgc            <  1:727066/150‑1 (MQ=255)                                                                                                                                          aTCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTgcg           >  1:1056253/1‑150 (MQ=255)                                                                                                                                           tCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTgcgc          >  1:927947/1‑150 (MQ=255)                                                                                                                                           tCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTgcgc          >  1:194146/1‑150 (MQ=255)                                                                                                                                              aGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTcgc                                                                                  <  2:768771/75‑1 (MQ=255)                                                                                                                                              aGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTcgc                                                                                  >  1:768771/1‑75 (MQ=255)                                                                                                                                              aGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCAAc       <  2:475200/150‑1 (MQ=255)                                                                                                                                              aGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCAAc       <  2:194146/150‑1 (MQ=255)                                                                                                                                               ggCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTACTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTgaggaaat      >  2:516443/1‑143 (MQ=255)                                                                                                                                                   gCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCAActgctg  >  2:289633/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                   TGCTGAGTAACGGTCTTCTCTCGTAGGAAACCTCCACAAGCCAGCCTTCGAGCTGGCTTGTGCGTATTCAAGCCATGACTATTTGGGAGTGAGCGCAGGACATGTCCCGCAATTCTTATATTCAGCAGAGCTCTGACATCCAGGCCGCAGGCGGCCTGGTCCCGATGGTTATCGAGCAGTCCGCCCGTGGCGAACGCGCTTACGACATCTACTCGCGCCTGTTGAAGGAGCGTGTGATCTTCCTGGTTGGCCCGGTAGAGGACTACATGGCCAACCTCGTCGTTGCGCAACTGCTG  >  NC_002947/2628281‑2628576

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