Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	466,450	G→T	6.3%	G150V (GGC→GTC)	PP_0385 →	serine acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	466,450	0	G	T	6.3%	202.1 / 3.4	63	G150V (GGC→GTC)	PP_0385	serine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (28/31);  new base T (2/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-01

TACGTCAGCATCGGGCGCAACTTCTTCATCCGGCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGGTACATTTGCGGACGTTGGCAGAGGGGATGG  >  NC_002947/466302‑466596                                                                                                                                                     |                                                                                                                                                   tACGTCAGCATCGGGCGCAACTTCTTCATCCGGCAGAACTGCTCGATCGGCATCAAGTCCCTGGGGTTGCAGGAGTACGCACTGTGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGc                                                                                                                                                   <  2:300267/150‑1 (MQ=255)    gTCAGCATCGGGCGCAACTTCTTCATCCGGCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGAt                                                                                                                                                <  2:531650/150‑1 (MQ=255)                 gcAACTTCTTCATCCGGCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCg                                                                                                                                   <  2:158313/150‑1 (MQ=255)                   aaCTTCTTCATCCGGCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGgc                                                                                                                                 >  2:225499/1‑150 (MQ=255)                        cttcATCCGGCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGAt                                                                                                                            >  1:654179/1‑150 (MQ=255)                        cttcATCCGGCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGAGATCGGCGCGAt                                                                                                                            <  2:1126119/150‑1 (MQ=255)                              ccGGCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGtt                                                                                                                      <  2:788840/150‑1 (MQ=255)                               cGGCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTg                                                                                                                     <  2:646942/150‑1 (MQ=255)                               cGGCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTg                                                                                                                     >  2:971838/1‑150 (MQ=255)                                ggCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGtt                                                                                                                      >  1:1443881/1‑148 (MQ=255)                                ggCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGtt                                                                                                                      <  2:1443881/148‑1 (MQ=255)                                   ggAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAg                                                                                                                                         <  1:1332831/125‑1 (MQ=255)                                   ggAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAg                                                                                                                                         >  2:1332831/2‑126 (MQ=255)                                       cTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTCGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTAcgcg                                                                                                             <  1:1567883/150‑1 (MQ=255)                                         gCTCGATCGGCATCAAGCCCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCggg                                                                                                           <  1:563467/150‑1 (MQ=255)                                         gCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGa                                                                                                                                                 >  1:695340/1‑112 (MQ=255)                                         gCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGa                                                                                                                                                 <  2:695340/112‑1 (MQ=255)                                          cTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGa                                                                                                          >  1:95971/1‑150 (MQ=255)                                             gATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGatat                                                                                                       >  1:712797/1‑150 (MQ=255)                                                  gCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGc                                                                                                  >  1:1083912/1‑150 (MQ=255)                                                      cAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTgg                                                                                              >  1:226468/1‑150 (MQ=255)                                                                gggTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTc                                                                                    <  2:1083912/150‑1 (MQ=255)                                                                       aGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAAcc                                                                             <  1:971838/150‑1 (MQ=255)                                                                               gCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGa                                                                                                                             <  2:524911/94‑1 (MQ=255)                                                                               gCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGa                                                                                                                             >  1:524911/1‑94 (MQ=255)                                                                               gCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAg                                                                     >  1:282045/1‑150 (MQ=255)                                                                                 aCTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGtt                                                                                                                      <  1:566197/99‑1 (MQ=255)                                                                                 aCTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGtt                                                                                                                      >  2:566197/1‑99 (MQ=255)                                                                                 aCTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGc                                                                                                  <  2:75405/119‑1 (MQ=255)                                                                                 aCTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGc                                                                                                  >  1:75405/1‑119 (MQ=255)                                                                                  cTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTc                                                                  <  2:654179/150‑1 (MQ=255)                                                                                     cgcATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCgag                                                               >  1:1163703/1‑150 (MQ=255)                                                                                      gcATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCgaga                                                              <  1:225499/150‑1 (MQ=255)                                                                                      gcATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCgaga                                                              <  2:1601998/150‑1 (MQ=255)                                                                                      gcATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCgaga                                                              >  1:586249/1‑150 (MQ=255)                                                                                            gTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTg                                                                                                                     >  2:659441/1‑89 (MQ=255)                                                                                            gTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTg                                                                                                                     <  1:659441/89‑1 (MQ=255)                                                                                              gACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATGTTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCg                                                      <  2:83388/150‑1 (MQ=255)                                                                                                 aaCGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGCCAATGTCGATgc                                                                                                                                              <  2:777009/59‑2 (MQ=37)                                                                                                 aaCGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGCCAATGTCGATgc                                                                                                                                              <  2:1602164/59‑2 (MQ=37)                                                                                                 aaCGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGCCAATGTCGATgc                                                                                                                                              >  1:1602164/1‑58 (MQ=255)                                                                                                 aaCGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGCCAATGTCGATgc                                                                                                                                              >  1:777009/1‑58 (MQ=255)                                                                                                  aCGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAg                                                  <  2:95971/150‑1 (MQ=255)                                                                                                   cGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAgg                                                 <  2:282045/150‑1 (MQ=255)                                                                                                         cATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTgg                                                                                                                    >  2:698192/1‑77 (MQ=255)                                                                                                         cATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTgg                                                                                                                    <  1:698192/77‑1 (MQ=255)                                                                                                         cATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGAGTTCGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTa                                           >  2:1050454/1‑150 (MQ=255)                                                                                                           tGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAAc                                         <  2:586249/150‑1 (MQ=255)                                                                                                              gcgcCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCaa                                      >  2:1498084/1‑150 (MQ=255)                                                                                                                  cAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTgg                                                                                                                    <  1:155322/68‑1 (MQ=255)                                                                                                                  cAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTgg                                                                                                                    >  2:155322/1‑68 (MQ=255)                                                                                                                          gCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGa                                                                                                          >  2:31335/1‑70 (MQ=255)                                                                                                                          gCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGa                                                                                                          <  1:31335/70‑1 (MQ=255)                                                                                                                                    tGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGCTACATTTGCGGACGtt                >  2:1109789/1‑150 (MQ=255)                                                                                                                                     gATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGGTACATTTGc                       <  2:1228698/142‑1 (MQ=255)                                                                                                                                     gATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGGTACATTTGc                       >  1:1228698/1‑142 (MQ=255)                                                                                                                                      aTCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGGTACATTTGCGGACGTTgg              <  1:97580/150‑1 (MQ=255)                                                                                                                                       tCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGGTACATTTGCGGACGTTGGc             <  2:712797/150‑1 (MQ=255)                                                                                                                                         gCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGGTACATTTGCGGACGTTGGCag           >  1:369398/1‑150 (MQ=255)                                                                                                                                          cATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGGTACATTTGCGGACGTTGGCaga          >  1:1029537/1‑150 (MQ=255)                                                                                                                                           aTCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGGTACATTTGCGGACGTTGGCagag         >  2:986415/1‑150 (MQ=255)                                                                                                                                               aCATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGGTATATTTGCGGACGTTGGCAGAGGGGa     <  1:1498084/150‑1 (MQ=255)                                                                                                                                                  ttGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCAACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGGTACATTTGCGGACGTTGGCAGAGGGGAtgg  >  1:521864/1‑150 (MQ=255)                                                                                                                                                     gCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAgggg                                               <  1:1343087/102‑1 (MQ=255)                                                                                                                                                     gCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAgggg                                               >  2:1343087/1‑102 (MQ=255)                                                                                                                                                     gCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCg                                 <  1:817155/116‑1 (MQ=255)                                                                                                                                                     gCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCg                                 >  2:817155/1‑116 (MQ=255)                                                                                                                                                     |                                                                                                                                                   TACGTCAGCATCGGGCGCAACTTCTTCATCCGGCAGAACTGCTCGATCGGCATCAAGACCCTGGGGTTGCAGGAGTACGCACTGCGCATTGGTGACAACGTGAGCATGGGCGCCAACAGTTGCATCGTTGCTGATCGCATCGACATTGGCGATGACGTAGTGATCGGCGCGATGTCGTTGGTTACGCGGGATATTCCGCCTGGCATGACCTTCTACAACCCACGGCAGGTCGAGATGCGTCCGGGAGGGGATTAACCAACAGCGGTACATTTGCGGACGTTGGCAGAGGGGATGG  >  NC_002947/466302‑466596

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