Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	972,118	T→A	7.9%	L136Q (CTG→CAG)	secD →	protein translocase subunit SecD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	972,118	0	T	A	7.9%	196.1 / 7.1	63	L136Q (CTG→CAG)	secD	protein translocase subunit SecD
Reads supporting (aligned to +/- strand):  ref base T (27/31);  new base A (2/3);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

CTGCCAGCCAAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGC  >  NC_002947/971970‑972265                                                                                                                                                     |                                                                                                                                                    gggccagcggAGGAGGCGGTGCGCAAGGCATTGGGGGATGATTACGTCGTGGCCCTGAACCTGGGCCAGACTACCACGCAATGGCTGCGCAACCGGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCctg                                                                                                                                                    <  2:327711/148‑1 (MQ=255)  tGCCAGCCAAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCctgc                                                                                                                                                   <  1:166269/150‑1 (MQ=255)   gccagccaAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCctgct                                                                                                                                                  <  2:507945/150‑1 (MQ=255)    ccagccaAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCctgctg                                                                                                                                                 >  1:406488/1‑150 (MQ=255)          aaGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtg                                                                                                                                           >  2:61571/1‑150 (MQ=255)           aGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtgg                                                                                                                                          <  2:814319/150‑1 (MQ=255)             gATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggac                                                                                                                                        <  2:604268/150‑1 (MQ=255)               tGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggacat                                                                                                                                      <  2:406488/150‑1 (MQ=255)                gTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggacatg                                                                                                                                     <  2:1057412/150‑1 (MQ=255)                 tAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggacatgg                                                                                                                                    <  1:331691/150‑1 (MQ=255)                    tGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTCCACTTCCAGCAGGAAGtgcacaccaccg                                                                                                                                 <  2:1422626/150‑7 (MQ=255)                    tGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTCCACTTCCAGCAGGAAGtgcacaccaccg                                                                                                                                 <  2:1422611/150‑7 (MQ=255)                             cACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGa                                                                                                                        <  2:932440/150‑1 (MQ=255)                                        gatTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggacat                                                                                                                                      <  1:1405890/125‑1 (MQ=255)                                        gatTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggacat                                                                                                                                      >  2:1405890/1‑125 (MQ=255)                                               tcgtGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtg                                                                                                                                           >  2:414114/1‑113 (MQ=255)                                               tcgtGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtg                                                                                                                                           <  1:414114/113‑1 (MQ=255)                                               tcgtGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGaaa                                                                                                          >  2:1323743/1‑146 (MQ=255)                                               tcgtGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGaaa                                                                                                          <  1:1323743/146‑1 (MQ=255)                                                 gtGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTAc                                                                                                    <  2:714430/150‑1 (MQ=255)                                                      ccTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAgg                                                                                               >  1:1295770/1‑150 (MQ=255)                                                        tGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCg                                                                                             >  2:71223/1‑150 (MQ=255)                                                          aaCCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAg                                                                                           <  2:1295770/150‑1 (MQ=255)                                                                gCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAg                                                                                     >  2:694737/1‑150 (MQ=255)                                                                 cccAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGTCCTGGACCTGTCCGGTGGTATGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGa                                                                                    <  1:804590/150‑1 (MQ=255)                                                                           cccGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAg                                                                                                                               >  2:844253/1‑97 (MQ=255)                                                                           cccGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAg                                                                                                                               <  1:844253/97‑1 (MQ=255)                                                                           cccGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTgcg                                                                          >  2:698497/1‑150 (MQ=255)                                                                            ccGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAg                                                                                                                               >  2:582357/1‑96 (MQ=255)                                                                            ccGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAg                                                                                                                               <  1:582357/96‑1 (MQ=255)                                                                                           aaCCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGaa                                                                                                                                              >  2:1063775/1‑66 (MQ=255)                                                                                           aaCCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGaa                                                                                                                                              <  1:1063775/66‑1 (MQ=255)                                                                                           aaCCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTtg                                                                               >  2:53329/1‑129 (MQ=255)                                                                                           aaCCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTtg                                                                               <  1:53329/129‑1 (MQ=255)                                                                                              cTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtgga                                                                                                                                         >  1:1078577/1‑68 (MQ=255)                                                                                              cTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtgga                                                                                                                                         <  2:1078577/68‑1 (MQ=255)                                                                                              cTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGa                                                                                                                        >  1:452537/1‑85 (MQ=255)                                                                                              cTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGa                                                                                                                        <  2:452537/85‑1 (MQ=255)                                                                                               tGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAgc                                                                  >  2:951759/1‑138 (MQ=255)                                                                                               tGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAgc                                                                  <  1:951759/138‑1 (MQ=255)                                                                                                    gcAAGCCCGATGAAACGGGGCCTCGACCTGCCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAgc                                                 <  1:559590/150‑1 (MQ=255)                                                                                                     cAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGGCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAgcc                                                >  1:5362/1‑150 (MQ=255)                                                                                                            gATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTca                                         >  1:997486/1‑150 (MQ=255)                                                                                                              tGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTcagc                                       >  2:1120837/1‑150 (MQ=255)                                                                                                                  aCTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGa                                   <  1:694737/150‑1 (MQ=255)                                                                                                                     gggCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCCCAGCatgatgc                                >  1:460727/1‑149 (MQ=255)                                                                                                                      ggCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATggc                               <  1:71223/150‑1 (MQ=255)                                                                                                                             aCCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGcatca                        >  1:1091825/1‑150 (MQ=255)                                                                                                                             aCCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGcatca                        >  1:899666/1‑150 (MQ=255)                                                                                                                               cTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggacatggac                                                                                                                                  >  1:1417876/1‑42 (MQ=255)                                                                                                                               cTGTCCGGTGGTGTGCACTTCCTGCTGGAAGtggacatggac                                                                                                                                  <  2:1417876/42‑1 (MQ=255)                                                                                                                               cTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATggcgg                             <  2:423903/143‑1 (MQ=255)                                                                                                                               cTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATggcgg                             >  1:423903/1‑143 (MQ=255)                                                                                                                                tGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCATGATGCCGCCATCCTGc                     <  1:555308/150‑1 (MQ=255)                                                                                                                                 gTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCt                    <  2:1268414/150‑1 (MQ=255)                                                                                                                                  tCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTAc                                                       <  2:207108/114‑1 (MQ=255)                                                                                                                                  tCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTAc                                                       >  1:207108/1‑114 (MQ=255)                                                                                                                                    cGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCCCAGCATGATGCCGCCATCctgctgag                 <  2:460727/150‑3 (MQ=255)                                                                                                                                         gtgtGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCg            <  2:1192076/150‑1 (MQ=255)                                                                                                                                          tgtgCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGc           <  1:698497/150‑1 (MQ=255)                                                                                                                                           gtgCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTcc                                                            >  2:836563/1‑100 (MQ=255)                                                                                                                                           gtgCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTcc                                                            <  1:836563/100‑1 (MQ=255)                                                                                                                                           gtccACTTCCAGCAGGAAGTGCACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTcagca                                      >  1:360915/4‑122 (MQ=255)                                                                                                                                           gtccACTTCCAGCAGGAAGTGCACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTcagca                                      <  2:360915/119‑1 (MQ=255)                                                                                                                                            tccACTTCCAGCAGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTg         >  2:463060/3‑150 (MQ=255)                                                                                                                                              cACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGGCTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGAc       >  1:1343263/1‑150 (MQ=255)                                                                                                                                                   ccTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGc  >  1:854569/1‑150 (MQ=255)                                                                                                                                                     tgctgGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTAc                                                                                                    <  1:676015/50‑1 (MQ=255)                                                                                                                                                     tgctgGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTAc                                                                                                    >  2:676015/1‑50 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CTGCCAGCCAAGGATGTAGTGCGCAAGGCACTGGGCGATGATTACGTCGTGGCCCTGAACCTGGCCCAGACTACCCCGCAATGGCTGCGCAACCTGGGCGCAAGCCCGATGAAACTGGGCCTGGACCTGTCCGGTGGTGTGCACTTCCTGCTGGAAGTGGACATGGACAAGGCCATGACTGCCCGCATGAAAGTCTACGAAGGCGAGGTCAAGACCTTGCTGCGCAAAGAGCGCGTCCGCTACCGCAGCCTGCCTCAGCAGGATGGCGGCATCATGCTGGGCTTCGCTGACGATGC  >  NC_002947/971970‑972265

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