Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,160,676	T→C	6.1%	*458R (TGA→CGA)	PP_3662 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,160,676	0	T	C	6.1%	187.8 / 2.0	65	*458R (TGA→CGA)	PP_3662	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/36);  new base C (2/2);  total (27/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

GCCTGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCC  >  NC_002947/4160527‑4160824                                                                                                                                                      |                                                                                                                                                     gCCTGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCt                                                                                                                                                      >  2:349314/1‑150 (MQ=255)   cTGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGa                                                                                                                                                    <  1:701994/150‑1 (MQ=255)    tGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGAc                                                                                                                                                   <  1:1290386/150‑1 (MQ=255)     gATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTTGCCACCTGata                                                                                                                                                  <  1:420117/150‑3 (MQ=255)     gATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTTGCCACCTGata                                                                                                                                                  <  2:880450/150‑3 (MQ=255)              acacGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCttt                                                                                                                                         <  1:857832/150‑1 (MQ=255)                acGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGGACCCTTCTATCAGCTGGCAACCTGCCCCGCTGCTTTaa                                                                                                                                       >  2:945105/1‑149 (MQ=255)                 cGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACt                                                                                                                                      <  2:358188/150‑1 (MQ=255)                     ccGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAgg                                                                                                                                  >  2:1332087/1‑150 (MQ=255)                        tATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACg                                                                                                                               >  2:405200/1‑150 (MQ=255)                         aTCAGATCCACGGTGAAAGCGCCGTACTACACCCGCTGGGCCGGCTTCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACgg                                                                                                                              <  2:373102/150‑1 (MQ=255)                          tCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACggg                                                                                                                             >  2:841020/1‑150 (MQ=255)                             gATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCAt                                                                                                                          <  1:559552/150‑1 (MQ=255)                               tCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATgg                                                                                                                        <  1:735145/150‑1 (MQ=255)                                        aCAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATCCATGCCCGTc                                                                                                               <  1:405200/150‑1 (MQ=255)                                            cgcCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAAc                                                                                                           <  2:1409515/150‑1 (MQ=255)                                             gcCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAAcg                                                                                                          >  1:1213452/1‑150 (MQ=255)                                                  aCTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCgggg                                                                                                     <  1:850060/150‑1 (MQ=255)                                                    tGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGc                                                                                                   <  2:1213452/150‑1 (MQ=255)                                                      gACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCCGACCCGCTGCTTTACTGAGGACGGGCCTGGATGCCCGTCGAACGCGCGGGTCg                                                                                                 >  2:1078397/1‑150 (MQ=255)                                                      gACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCCGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCg                                                                                                 >  2:1078384/1‑150 (MQ=255)                                                       aCCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCgg                                                                                                <  1:349314/150‑1 (MQ=255)                                                       aCCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCgg                                                                                                <  2:1227228/150‑1 (MQ=255)                                                        cccGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCggg                                                                                               <  1:908429/150‑1 (MQ=255)                                                           gCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGa                                                                                            >  2:616114/1‑150 (MQ=255)                                                                gCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGTTAAAGCTGGCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCCGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGa                                                                                       <  1:1078397/150‑1 (MQ=255)                                                                gCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCCGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGa                                                                                       <  1:1078384/150‑1 (MQ=255)                                                                          gcAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCCTGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTAc                                                                             <  1:1014305/150‑1 (MQ=255)                                                                           cAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACg                                                                            >  2:228813/1‑150 (MQ=255)                                                                            aGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACgg                                                                           >  2:735633/1‑150 (MQ=255)                                                                            aGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACgg                                                                           <  1:735633/150‑1 (MQ=255)                                                                                aTTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATCCATGCCCGTc                                                                                                               >  1:1019547/1‑110 (MQ=255)                                                                                aTTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATCCATGCCCGTc                                                                                                               <  2:1019547/110‑1 (MQ=255)                                                                                 ttCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCAc                                                                      <  1:475525/150‑1 (MQ=255)                                                                                         ccAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATgg                                                              <  1:1321598/150‑1 (MQ=255)                                                                                          cAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGCATGTACCGTGCACAACGGGGGCAACCTGACCCGCTGCTGGACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATggg                                                             >  1:958433/1‑150 (MQ=255)                                                                                            gCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGtt                                                           <  1:945105/150‑1 (MQ=255)                                                                                              aCCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGa                                                                                       <  2:1249574/120‑1 (MQ=255)                                                                                              aCCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGa                                                                                       >  1:1249574/1‑120 (MQ=255)                                                                                              aCCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATCCATGCCCGTc                                                                                                               >  1:830197/1‑96 (MQ=255)                                                                                              aCCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATCCATGCCCGTc                                                                                                               <  2:830197/96‑1 (MQ=255)                                                                                                cGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACAt                                                       <  2:958433/150‑1 (MQ=255)                                                                                                    tGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGaaa                                                   <  1:616114/150‑1 (MQ=255)                                                                                                        gCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTc                                               >  1:251790/1‑150 (MQ=255)                                                                                                         cTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGGCCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTcc                                              <  1:606912/150‑1 (MQ=255)                                                                                                             ccGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGccc                                                                                                                  >  2:61064/1‑78 (MQ=255)                                                                                                             ccGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGccc                                                                                                                  <  1:61064/78‑1 (MQ=255)                                                                                                                  ggcgcggcgTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAg                                     <  1:841020/150‑1 (MQ=255)                                                                                                                    cgcggcgTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGt                                   >  1:1091952/1‑150 (MQ=255)                                                                                                                            aTGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCt                           >  2:806821/1‑150 (MQ=255)                                                                                                                                  cGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACgggca                                                                                      >  2:83700/1‑83 (MQ=255)                                                                                                                                  cGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACgggca                                                                                      <  1:83700/85‑3 (MQ=255)                                                                                                                                  cGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTa                                                          <  2:434753/113‑1 (MQ=255)                                                                                                                                  cGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGGTa                                                          >  1:434753/1‑113 (MQ=255)                                                                                                                                      cTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGa                 <  1:228813/150‑1 (MQ=255)                                                                                                                                      cTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGa                 <  1:1332087/150‑1 (MQ=255)                                                                                                                                       tATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAg                >  1:499380/1‑150 (MQ=255)                                                                                                                                       tATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAg                >  1:564209/1‑150 (MQ=255)                                                                                                                                             gTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATgg                                                              <  1:242366/98‑1 (MQ=255)                                                                                                                                             gTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATgg                                                              >  2:242366/1‑98 (MQ=255)                                                                                                                                              tGGCAACCTGACCCGCTGCTTTACTGAGGACGGGTATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTAc                                                                             <  1:1069478/82‑1 (MQ=255)                                                                                                                                              tGGCAACCTGACCCGCTGCTTTACTGAGGACGGGTATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTAc                                                                             >  2:1069478/1‑82 (MQ=255)                                                                                                                                                  aaCCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAg                                     >  2:344341/1‑118 (MQ=255)                                                                                                                                                  aaCCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAg                                     <  1:344341/118‑1 (MQ=255)                                                                                                                                                    ccTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGACAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCAccgcc   >  2:452672/1‑150 (MQ=255)                                                                                                                                                     cTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCAccgccc  >  1:1200010/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GCCTGATCGAGGAACACGGGCCGTATCAGATCCACGGTGACAGCGCCGTACTGGACCCGCTGGGCCGGCTGCTGCAGGCATTCGTCGACCAGCACCGGATGAAGCTGCCCGGCGGCGCGGCGTATGTACCGTGCTATCAGGTGGCAACCTGACCCGCTGCTTTACTGAGGACGGGCATGGATGCCCGTCGAACGCGGGGGCCGGGTGACGGGAGCGTCTGTACGGATCACCCACATGGGTTACATGAAAGCTCCCCCACCTAGGTGACTCGCTCACTGCCGGAGTTCAACCACCGCCC  >  NC_002947/4160527‑4160824

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