Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,108,759	G→C	7.4%	A312G (GCG→GGG)	PP_5359 ←	CobW/P47K family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,108,759	0	G	C	7.4%	198.3 / 4.8	64	A312G (GCG→GGG)	PP_5359	CobW/P47K family protein
Reads supporting (aligned to +/- strand):  ref base G (26/32);  new base C (3/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACTCCCTTGGCCCGTCGC  >  NC_002947/6108617‑6108905                                                                                                                                               |                                                                                                                                                   aTCGAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGggttg                                                                                                                                             >  2:1317917/1‑150 (MQ=255)  tCGAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGggttgc                                                                                                                                            <  2:1250261/150‑1 (MQ=255)   cGAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGggttgcg                                                                                                                                           <  1:974451/150‑1 (MQ=255)         aGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGgttgcggttgcg                                                                                                                                     >  1:349736/1‑150 (MQ=255)          gcgcCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGCGCATAACCGCTGGCGCAACTGGTCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGttgcggttgcgg                                                                                                                                    >  1:882858/1‑150 (MQ=255)              ccTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGa                                                                                                                                <  2:250395/150‑1 (MQ=255)              ccTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGCCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGa                                                                                                                                <  1:377752/150‑1 (MQ=255)                tACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATg                                                                                                                              >  1:1058966/1‑150 (MQ=255)                   ccGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAAATGATCCACAGAAGAGCTCGTCCAGGAGCTCGGCGTGCAACTCGGAGTGCGTTTGAGGCTGATGGaa                                                                                                                           >  2:1085360/1‑150 (MQ=255)                    cGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTGCGAGATATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGaaa                                                                                                                          <  1:1085360/150‑1 (MQ=255)                      ttCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAAt                                                                                                                        >  2:258579/1‑150 (MQ=255)                        cGAAATCAAGACCCCACCACGCGCAGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCa                                                                                                                      <  2:472944/150‑1 (MQ=255)                            aTCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACCCGcgttgcaag                                                                                                                                         >  2:1194014/1‑124 (MQ=255)                            aTCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACCCGcgttgcaag                                                                                                                                         <  1:1194014/127‑4 (MQ=255)                                aGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATa                                                                                                              <  1:1221004/150‑1 (MQ=255)                                   ccccACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACgg                                                                                                           <  1:43148/150‑1 (MQ=255)                                    cccACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGc                                                                                                          >  1:1325492/1‑150 (MQ=255)                                       accacGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGt                                                                                                       <  2:623242/150‑1 (MQ=255)                                             cgcCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTgc                                                                                                 <  2:179167/150‑1 (MQ=255)                                                  aCAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAg                                                                                            <  1:717739/150‑1 (MQ=255)                                                       gTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCt                                                                                       <  2:467131/150‑1 (MQ=255)                                                        tCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGGTCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTg                                                                                      <  2:882858/150‑1 (MQ=255)                                                         cagagagATGCAGCAAAGGCGCCAGGCGCCTAACCGCTGGCGCAACTGATCACCTGACGCGACATTCATTCAACCCGTCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTgg                                                                                     <  2:1008142/147‑1 (MQ=255)                                                         cATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTgg                                                                                     <  2:1317414/150‑1 (MQ=255)                                                           tAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGAGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGc                                                                                   >  2:157258/1‑150 (MQ=255)                                                                   gagcaATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAg                                                                                            >  1:205289/2‑133 (MQ=255)                                                                   cagcaATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCt                                                                                                                                  <  1:126546/95‑1 (MQ=255)                                                                   cagcaATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAg                                                                                            <  2:205289/133‑1 (MQ=255)                                                                   cagcaATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGAGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCt                                                                                                                                  >  2:126546/1‑95 (MQ=255)                                                                     gcaATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGCAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCACTGGGCTGCCAAGCgg                                                                         <  2:1058966/150‑1 (MQ=255)                                                                                aGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCa                                                              <  1:1317917/150‑1 (MQ=255)                                                                                    gCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCa                                                          >  2:100205/1‑150 (MQ=255)                                                                                       tAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGcc                                                       >  1:1122785/1‑150 (MQ=255)                                                                                                 cgcAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGa                                             <  2:1203038/150‑1 (MQ=255)                                                                                                  gcAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGac                                            >  1:1050317/1‑150 (MQ=255)                                                                                                   cAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATcc                                                                 >  1:890632/1‑128 (MQ=255)                                                                                                   cAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATcc                                                                 <  2:890632/128‑1 (MQ=255)                                                                                                    aaCTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGacca                                          >  1:220735/1‑150 (MQ=255)                                                                                                    aaCTGATCAACTGACGCGACATTCATTCAACCCGCATTGCAAAGCGGGTCCGATTTGCGTCTGAGGGAAGAACAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGacca                                          >  1:388348/1‑150 (MQ=255)                                                                                                     aCTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGaccac                                         <  1:1408589/150‑1 (MQ=255)                                                                                                             aCTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTgcgc                                                                                               >  2:670017/1‑88 (MQ=255)                                                                                                             aCTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTgcgc                                                                                               <  1:670017/88‑1 (MQ=255)                                                                                                             aCTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTc                                 >  1:203923/1‑150 (MQ=255)                                                                                                                     gACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCtgt                         <  2:1271651/150‑1 (MQ=255)                                                                                                                      aCATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCtgtg                        <  1:258579/150‑1 (MQ=255)                                                                                                                        attcattcAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGACAAGAGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGa                      >  2:1107544/1‑150 (MQ=255)                                                                                                                         ttcattcaACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGAt                     >  1:588382/1‑150 (MQ=255)                                                                                                                         ttcattcaACCCAGCTTGCAACCCGGGCTGCGGTTGCGGCCGACGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCGACGCGGGTGTATCCCCATCCGACCGCCTCAAGAACCACCAACCATCCATGCTGTAGAt                     >  2:373958/1‑150 (MQ=255)                                                                                                                            attcaACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAAc                  >  1:399329/1‑150 (MQ=255)                                                                                                                             ttcaACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACt                 >  1:184871/1‑150 (MQ=255)                                                                                                                             cgcaACCCGCGTTGCAAGCCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCg                                                                          <  2:381248/91‑1 (MQ=38)                                                                                                                             cgcaACCCGCGTTGCAAGCCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCg                                                                          >  1:381248/3‑93 (MQ=38)                                                                                                                               caACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACTcc               <  1:985991/150‑1 (MQ=255)                                                                                                                                  cccGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACTCCCtt            <  2:1325492/150‑1 (MQ=255)                                                                                                                                  cccGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACTCCCtt            >  1:164115/1‑150 (MQ=255)                                                                                                                                   ccGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCCCATGCCCTCCATGCTGTGGATAACTCCCTTg           <  1:1107544/150‑1 (MQ=255)                                                                                                                                    cGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGcc                                                       <  2:135080/105‑1 (MQ=255)                                                                                                                                    cGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGcc                                                       >  1:135080/1‑105 (MQ=255)                                                                                                                                       cTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGAt                     <  1:964312/136‑1 (MQ=255)                                                                                                                                       cTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGAt                     >  2:964312/1‑136 (MQ=255)                                                                                                                                          gCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACTCCCTTTGCCCGTc    <  2:1122785/150‑1 (MQ=255)                                                                                                                                          gCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACTCCCTTGGCCCGTc    <  2:588382/150‑1 (MQ=255)                                                                                                                                            aaCGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACTCCCTTGGCCCGTCGc  >  2:497403/1‑150 (MQ=255)                                                                                                                                            aaCGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACTCCCTTGGCCCGTCGc  >  1:174675/1‑150 (MQ=255)                                                                                                                                               |                                                                                                                                                   ATCGAACTAGCGCCCTACCCGTTCGAAATCAAGACCCCACCACGCGCCGACAACGTCATAGAGATGCAGCAATGGCGCCAGGAGCATAACCGCTGGCGCAACTGATCAACTGACGCGACATTCATTCAACCCGGCTTGCAACGCGGGTTGCGGTTGCGGCTGATGGAAAATCAGCTCGATACGGCTGTCCCTGCGCCAGTCGCTGGGCTGCCAAGCGGGCATATCCCCATTCAGCCCATTGAATGACCACCAGCCCTCCATGCTGTGGATAACTCCCTTGGCCCGTCGC  >  NC_002947/6108617‑6108905

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