Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,701,700	G→A	6.6%	S29S (TCG→TCA)	PP_3267 →	Clp protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,701,700	0	G	A	6.6%	174.2 / 2.7	61	S29S (TCG→TCA)	PP_3267	Clp protease
Reads supporting (aligned to +/- strand):  ref base G (32/25);  new base A (2/2);  total (34/27)
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.37e-01

GGTGCCAGGATGACATTCGAGCGTTTCGCCCACCGCCAGAGGGCCGGGTTGGGCGACGATGGGCAAAAGATCAAGCAATGACCTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTACCCTGTACAACTTTCTGCGGGCATTACCCGTGGAGGTGCA  >  NC_002947/3701556‑3701849                                                                                                                                                 |                                                                                                                                                      ggTGCCAGGATGACATTCGAGCGTTTCGCCCACCGCCAGAGGGCCGGGTTGGGCGACGATGGGCAAAAGATCAAGCAATGACCTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAAc                                                                                                                                                  <  2:912557/150‑1 (MQ=255)               aTTCGAGCGTTTCGCCCACCGCCAGAGGGCCGGGTTGGGCGACGATGGGCAAAAGATCAAGCAATGACCTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGa                                                                                                                                    >  2:700219/1‑150 (MQ=255)                     gCGTTTCGCCCACCGCCAGAGGGCCGGGTTGGGCGACGATGGGCAAAAGATCAAGCAATGACCTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGa                                                                                                                              >  2:57337/1‑150 (MQ=255)                       gTTTCGCCCACCGCCAGAGGGCCGGGTTGGGCGACGATGGGCAAAAGATCAAGCAATGACCTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGAcc                                                                                                                            >  1:704497/1‑150 (MQ=255)                         ttCGCCCACCGCCAGAGGGCCGGGTTGGGCGACGATGGGCAAAAGATCAAGCAATGACCTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAAc                                                                                                                                        >  1:583657/1‑136 (MQ=255)                         ttCGCCCACCGCCAGAGGGCCGGGTTGGGCGACGATGGGCAAAAGATCAAGCAATGACCTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAAc                                                                                                                                        <  2:583657/136‑1 (MQ=255)                                 ccGCCAGAGGGCCGGGTTGGGCGACGATGGGCAAAAGATCAAGCAATGACCTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCggg                                                                                                                  <  1:1370515/150‑1 (MQ=255)                                   gCCAGAGGGCCGGGTTGGGCGACGATGGGCAAAAGATCAAGCAATGACCTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCa                                                                                                                >  2:692901/1‑150 (MQ=255)                                   gCCAGAGGGCCGGGTTGGGCGACGATGGGCAAAAGATCAAGCAATGACCTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCa                                                                                                                >  2:1214304/1‑150 (MQ=255)                                      agagGGCCGGGTTGGGCGACGATGGGCAAAAGATCAAGCAATGACCTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTg                                                                                                             <  2:888158/150‑1 (MQ=255)                                          ggCCGGGTTGGGCGACGATGGGCAAAAGATCAAGCAATGACCTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAAc                                                                                                         >  1:506885/1‑150 (MQ=255)                                             cGGGTTGGGCGACGATGGGCAACCGATGAAGCAACGTCCTGCCACCCGTGGAGGCAAGCATGGCCAGCCATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAgg                                                                                                      <  2:435855/150‑1 (MQ=255)                                             cGGGTTGGGCGACGATGGGCAAAAGATCCAGCAATGTCCTGACACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGAACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGCCCTGCTCGCGGGCATTGCAACAgg                                                                                                      <  2:184397/150‑1 (MQ=255)                                             cGGGTTGGGCGACGATGGGCAAAAGATCAAGCAATGACCTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAgg                                                                                                      >  1:610948/1‑150 (MQ=255)                                             cGGGTTGGGCGACGATGGGCAAAAGATCAAGCAATGACCTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAgg                                                                                                      <  1:57337/150‑1 (MQ=255)                                             cGGGTTGGGCGACGATGGGCAAAAGATCAAGCAATGACCTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAgg                                                                                                      >  1:866639/1‑150 (MQ=255)                                                         cgaTGGGCAAAAGATCAAGCAATGACCTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTc                                                                                          <  1:1214304/150‑1 (MQ=255)                                                               gCAAAAGATCAAGCAATGACCTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAAc                                                                                                         >  1:1023297/1‑129 (MQ=255)                                                               gCAAAAGATCAAGCAATGACCTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTGCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAAc                                                                                                         <  2:1023297/129‑1 (MQ=255)                                                                aaaaaGATCAAGCAATGACCTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACt                                                                                   <  2:217763/149‑1 (MQ=255)                                                                 aaaaGATCAAGCAATGACCTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTg                                                                                  <  2:704497/150‑1 (MQ=255)                                                                            cAATGACCTGCCACCCGTGGAGGCAAGCAGGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGccac                                                                       >  1:1328783/1‑150 (MQ=255)                                                                                 aCCTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATgg                                                                  <  2:610948/150‑1 (MQ=255)                                                                                   cTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGggc                                                                >  1:817618/1‑150 (MQ=255)                                                                                     gCCACCCGTGGAGGCAAACATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGggcgg                                                              >  2:1029484/1‑150 (MQ=255)                                                                                         cccGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATggg                                                                 <  1:1198603/143‑1 (MQ=255)                                                                                         cccGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATggg                                                                 >  2:1198603/1‑143 (MQ=255)                                                                                         cccGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGaa                                                          <  1:700219/150‑1 (MQ=255)                                                                                         cccGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGaa                                                          >  2:179386/1‑150 (MQ=255)                                                                                           cGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGAc                                                                                                                             <  2:1018954/81‑1 (MQ=255)                                                                                           cGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGAc                                                                                                                             >  1:1018954/1‑81 (MQ=255)                                                                                            gTGGAGGCAAGCATGGCCAGATATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGc                                                       >  1:474679/1‑150 (MQ=255)                                                                                                    aaGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATgg                                               >  2:1197144/1‑150 (MQ=255)                                                                                                      gCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGt                                             <  1:692901/150‑1 (MQ=255)                                                                                                        aTGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGtt                                            <  2:102429/149‑1 (MQ=255)                                                                                                        aTGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGtt                                            >  1:102429/1‑149 (MQ=255)                                                                                                          ggCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTAc                                         <  2:506885/150‑1 (MQ=255)                                                                                                           gCCAGACATATCATCCACTTCACCGGACCGATCAATTCATCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGggcggc                                                             <  2:112742/129‑1 (MQ=255)                                                                                                           gCCAGACATATCATCCACTTCACCGGACCGATCAATTCATCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGggcggc                                                             <  2:112468/129‑1 (MQ=255)                                                                                                           gCCAGACATATCATCCACTTCACCGGACCGATCAATTCATCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGggcggc                                                             >  1:112468/1‑129 (MQ=255)                                                                                                           gCCAGACATATCATCCACTTCACCGGACCGATCAATTCATCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGggcggc                                                             >  1:112742/1‑129 (MQ=255)                                                                                                                   tatCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTACCCTGTACaa                                >  2:436766/1‑150 (MQ=255)                                                                                                                        tcCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGc                                                                                                  >  1:472089/1‑79 (MQ=255)                                                                                                                        tcCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGc                                                                                                  <  2:472089/79‑1 (MQ=255)                                                                                                                        tcCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTACCCTGTACAACTTTc                           >  1:176316/1‑150 (MQ=255)                                                                                                                             ttCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTACCCTGTACAACTTTCTGCgg                      <  2:817618/150‑1 (MQ=255)                                                                                                                              tCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTACCCTGTACAACTTTCTGCggg                     <  2:1328783/150‑1 (MQ=255)                                                                                                                                         aTCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTACCCTGTACAACTTTCTGCgg                      >  1:1311809/1‑138 (MQ=255)                                                                                                                                         aTCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTACCCTGTACAACTTTCTGCgg                      <  2:1311809/138‑1 (MQ=255)                                                                                                                                         aTCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTACCCTGTACAACTTTCTGCGGGCATTACCCGTg          >  1:181465/1‑150 (MQ=255)                                                                                                                                         aTCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTACCCTGTACAACTTTCTGCGGGCATTACCCGTg          >  2:735832/1‑150 (MQ=255)                                                                                                                                         aTCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTACCCTGTACAACTTTCTGCGGGCATTACCCGTg          >  2:1238014/1‑150 (MQ=255)                                                                                                                                          tCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTACCCTGTACAACTTTCTGCGGGCATTACCCGTgg         >  2:1340967/1‑150 (MQ=255)                                                                                                                                            aaTTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTACCCTGTACAACTTTCTGCGGGCATTACCCGTGGa        <  1:237495/149‑1 (MQ=255)                                                                                                                                            aaTTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTACCCTGTACAACTTTCTGCGGGCATTACCCGTGGa        >  2:237495/1‑149 (MQ=255)                                                                                                                                             aTTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTACCCTGTACAACTTTCTGCGGGCATTACCCg            <  1:1424582/144‑1 (MQ=255)                                                                                                                                             aTTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTACCCTGTACAACTTTCTGCGGGCATTACCCg            >  2:1424582/1‑144 (MQ=255)                                                                                                                                             aTTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTACCCTGTACAACTTTCTGCGGGCATTACCCGTGGAgg      >  2:599718/1‑150 (MQ=255)                                                                                                                                             aTTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTACCCTGTACAACTTTCTGCGGGCATTACCCGTGGAgg      >  1:358114/1‑150 (MQ=255)                                                                                                                                             aTTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATATGAACCATGGGCGGCGAATGCAGCTATGGGTTTACACTGTCAAACATTCTGCGGTCATTACCCTTGGAgg      >  1:1022469/1‑150 (MQ=255)                                                                                                                                              ttCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTACCCTGTACAACTTTCTGCGGGCATTACCCGTGGAGGt     <  1:1197144/150‑1 (MQ=255)                                                                                                                                                 gTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTACCCTGTACAACTTTCTGCGGGCATTACCCGTGGAGGTGCa  >  2:529418/1‑150 (MQ=255)                                                                                                                                                 gTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTACCCTGTACAACTTTCTGCGGGCATTACCCGTGGAGGTGCa  >  2:301897/1‑150 (MQ=255)                                                                                                                                                 gTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTACCCTGTACAACTTTCTGCGGGCATTACCCGTGGAGGTGCa  >  1:584098/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                      GGTGCCAGGATGACATTCGAGCGTTTCGCCCACCGCCAGAGGGCCGGGTTGGGCGACGATGGGCAAAAGATCAAGCAATGACCTGCCACCCGTGGAGGCAAGCATGGCCAGACATATCATCCACTTCACCGGACCGATCAATTCGTCAACCTGCGGCAACCTGATCACGACCTGCTCGCGGGCATTGCAACAGGGTGCCGAGCTTCTGCAACTGAACATCGCCACCATGGGCGGCGAATGCAGCTATGGGTTTACCCTGTACAACTTTCTGCGGGCATTACCCGTGGAGGTGCA  >  NC_002947/3701556‑3701849

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