Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,640,335	C→T	6.9%	A122V (GCG→GTG)	PP_3207 →	cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,640,335	0	C	T	6.9%	168.7 / 2.7	58	A122V (GCG→GTG)	PP_3207	cyclase
Reads supporting (aligned to +/- strand):  ref base C (28/26);  new base T (2/2);  total (30/28)
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.75e-01

TCACATGGACGCGCCCTGGCACTACGCATCCACCACCGACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAAGCCGGAGTGGGGATGGGCCGCGAGGCGACCTTG  >  NC_002947/3640189‑3640480                                                                                                                                                    |                                                                                                                                                  tCACATGGACGCGCCCTGGCACTACGCATCCACCACCGACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTgcgg                                                                                                                                                 >  5:37726/1‑149 (MQ=255) tCACATGGACGCGCCCTGGCACTACGCATCCACCACCGACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCACTTGCCGGATGGCCACGTGGTCACTgcgg                                                                                                                                                 >  6:168739/1‑149 (MQ=255)   acaTGGACGCGCCCTGGCACTACGCATCCACCACCGACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTgcggcg                                                                                                                                               >  4:306941/1‑149 (MQ=255)       ggACGCGCCCTGGCACTACGCATCCACCACCGACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAgg                                                                                                                                           >  1:172677/1‑149 (MQ=255)           gcgcCCTGGCACTACGCATCCACCACCGACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTc                                                                                                                                         >  4:59645/1‑147 (MQ=255)                   gCACTACGCATCCACCACCGACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCACTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGc                                                                                                                               <  5:168739/149‑1 (MQ=255)                       tACGCATCCACCACCGACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTgc                                                                                                                           >  4:266903/1‑149 (MQ=255)                         cGCATCCACCACCGACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTgcgc                                                                                                                         <  2:128727/149‑1 (MQ=255)                         cGCATCCACCACCGACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTgcgc                                                                                                                         <  6:37726/149‑1 (MQ=255)                          gCATCCACCACCGACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGcgcg                                                                                                                        <  8:152903/149‑1 (MQ=255)                           cATCCACCACCGACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGcgcg                                                                                                                        >  2:72089/1‑148 (MQ=255)                               caccacCGACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTg                                                                                                                   >  8:298507/1‑149 (MQ=255)                               caccacCGACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTg                                                                                                                   >  2:37902/1‑149 (MQ=255)                               caccacCGACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTg                                                                                                                   >  7:29462/1‑149 (MQ=255)                               caccacCGACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTg                                                                                                                   <  2:172677/149‑1 (MQ=255)                                 ccaccGACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTggg                                                                                                                 <  1:37902/149‑1 (MQ=255)                                   accGACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCa                                                                                                               >  7:160847/1‑149 (MQ=255)                                   accGACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGACGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGGCGAGGCCGAGCTTGCGCGGATTGGGCa                                                                                                               >  4:143215/1‑149 (MQ=255)                                     cGACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATg                                                                                                             <  7:298507/149‑1 (MQ=255)                                      gACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGTATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGa                                                                                                            <  3:59645/149‑1 (MQ=255)                                      gACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGa                                                                                                            <  3:266903/149‑1 (MQ=255)                                               aaGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGcc                                                                                                   >  4:106810/1‑149 (MQ=255)                                               aaGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGcc                                                                                                   >  3:177380/1‑149 (MQ=255)                                                 gCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGc                                                                                                 <  3:306941/149‑1 (MQ=255)                                                  ccGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCCAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCATCAGCt                                                                                                >  8:216062/1‑149 (MQ=255)                                                       cTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACa                                                                                           >  6:154938/1‑149 (MQ=255)                                                          cGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATc                                                                                         <  4:177380/148‑1 (MQ=255)                                                            gCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGTGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGa                                                                                                            >  6:22000/1‑127 (MQ=255)                                                            gCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGTGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGa                                                                                                            <  5:22000/127‑1 (MQ=255)                                                             cATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGa                                                                                                                                 >  7:288171/1‑105 (MQ=255)                                                             cATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGa                                                                                                                                 <  8:288171/105‑1 (MQ=255)                                                             cATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGc                                                                                     <  1:72089/149‑1 (MQ=255)                                                               tCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCtt                                                                                   <  8:160847/149‑1 (MQ=255)                                                                        tGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTg                                                                                      >  5:76489/1‑137 (MQ=255)                                                                        tGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTg                                                                                      <  6:76489/137‑1 (MQ=255)                                                                          ccGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCcg                                                                        >  5:118547/1‑149 (MQ=255)                                                                                  gtggtgCCTGCAGCCGGGGGTGAAGTTGGATTTTGCGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTgg                                                                                                                  >  1:202300/1‑100 (MQ=255)                                                                                  gtggtgCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTgg                                                                                                                  <  2:202300/99‑1 (MQ=255)                                                                                     gtgCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGc                                                              <  8:29462/148‑1 (MQ=255)                                                                                     gtgCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCt                                                             >  5:136294/1‑149 (MQ=255)                                                                                      tgCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGTGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCtt                                                            >  7:234795/1‑149 (MQ=255)                                                                                         cTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGtt                                                         >  2:279875/1‑149 (MQ=255)                                                                                          tGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTtgt                                                                                 <  1:69385/124‑1 (MQ=255)                                                                                          tGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTtgt                                                                                 >  2:69385/1‑124 (MQ=255)                                                                                             aGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGc                                                     <  7:216062/149‑1 (MQ=255)                                                                                             aGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGc                                                     <  3:106810/149‑1 (MQ=255)                                                                                             aGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGc                                                     <  6:118547/149‑1 (MQ=255)                                                                                                cGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAg                                                  >  4:57109/1‑149 (MQ=255)                                                                                                     gTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTgg                                             >  4:165221/1‑149 (MQ=255)                                                                                                          gTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTAtc                                        <  6:136294/149‑1 (MQ=255)                                                                                                              gATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCg                                     >  7:128379/1‑148 (MQ=255)                                                                                                                 ttt‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGc                                                 >  7:182639/1‑133 (MQ=255)                                                                                                                 ttt‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGc                                                 <  8:182639/133‑1 (MQ=255)                                                                                                                       tCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAAGCCGGAg                            <  5:154938/149‑1 (MQ=255)                                                                                                                       tCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAAGCCGGAg                            >  1:311016/1‑149 (MQ=255)                                                                                                                         aTTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAAGCCGGAGTg                          >  3:10692/1‑149 (MQ=255)                                                                                                                          tttGCCGGATGGCCACGTGGTCACTGTGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAAGCCGGAGTgg                         <  8:234795/149‑1 (MQ=255)                                                                                                                          tttGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAAGCCGGAGTgg                         >  7:167405/1‑149 (MQ=255)                                                                                                                             gCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAAGCCGGAGTgggg                       <  1:279875/148‑1 (MQ=255)                                                                                                                              ccGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAAGCCGGAGTGGGGAt                     <  3:57109/149‑1 (MQ=255)                                                                                                                                                  tGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAAGCCGGAGTGGGGATGGGCAGCGAGGCGACCTTg  >  1:54315/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                  TCACATGGACGCGCCCTGGCACTACGCATCCACCACCGACGGTGGCAAGCCGGCCTTCGGCATCGACGAGCTGCCGCTGGAGTGGTGCCTGCAGCCGGGGGTGAAGTTGGATTTT‑CGTCATTTGCCGGATGGCCACGTGGTCACTGCGGCGCAGGTCGAGGCCGAGCTTGCGCGTATTGGGCATGAGTTGCAGCCGCTGGACATCGTGCTTGTGAACACCCGCGCGGGCAGCTTGTTCGGCCAGCCTGGCTATCTCGAAGCCGGAGTGGGGATGGGCCGCGAGGCGACCTTG  >  NC_002947/3640189‑3640480

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