Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,581,445	C→A	6.8%	F166L (TTC→TTA)	PP_3160 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,581,445	0	C	A	6.8%	189.1 / 2.7	59	F166L (TTC→TTA)	PP_3160	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/27);  new base A (2/2);  total (30/29)
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.01e-01

TTCAATGCGCCAGACAGCTTGAGCCCGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATGAATTCGACAGTACCCTCCCAGCCAGGATAAAACACGC  >  NC_002947/3581305‑3581594                                                                                                                                             |                                                                                                                                                       ttCAATGCGCCAGACAGCTTGAGCCCGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTTACTATCCGCCCAGCAGCTTCAT‑TGCgggc                                                                                                                                              <  1:807339‑M2/150‑5 (MQ=255)  tCAATGCGCCAGACAGCTTGAGCCCGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGTTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCgggcg                                                                                                                                             >  1:825302‑M2/1‑145 (MQ=255)     aTGCGCCAGACAGCTTGAGCCCGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGtt                                                                                                                                          <  2:1817783/150‑1 (MQ=255)      tGCGCCAGACAGCTTGAGCCCGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCAAGTTAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGttt                                                                                                                                         <  1:233994/150‑1 (MQ=255)           cAGACAGCTTGAGCCCGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGtttattta                                                                                                                                    <  1:648345/150‑1 (MQ=255)            aGACAGCTTGAGCCCGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTAc                                                                                                                                   >  2:798313/1‑150 (MQ=255)               cAGCTTGAGCCCGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCaa                                                                                                                                <  1:658567/150‑1 (MQ=255)                  cTTGAGCCCGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCg                                                                                                                             >  1:986778/1‑150 (MQ=255)                         ccGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGa                                                                                                                      <  2:1524825/150‑1 (MQ=255)                         ccGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGa                                                                                                                      <  1:452757/150‑1 (MQ=255)                         ccGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCC‑GCAA‑TGAAGCTGCtgggcggatagtcaacgtgaccagttgga                                                                                                                     <  2:1536825‑M2/150‑30 (MQ=255)                                 tGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTc                                                                                                              <  2:672731/150‑1 (MQ=255)                                    tGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCC‑GCAATGAAGC‑TGCtgggcggatagtcaacgtgaccagttggaacgtgaccact                                                                                                          >  1:1536825‑M2/1‑110 (MQ=255)                                        ggggCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGca                                                                                                       <  1:518296/150‑1 (MQ=255)                                             cAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCata                                                                                                  >  1:1310357/1‑150 (MQ=255)                                             cAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCC‑GCAA‑TGAAGCTGCtgggcggatagtcaacgtgaccagttggaacgt                                                                                                                 >  2:1449680‑M2/1‑101 (MQ=255)                                             cAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCC‑GCAA‑TGAAGCTGCtgggcggatagtcaacgtgaccagttggaacgt                                                                                                                 <  1:1449680‑M2/134‑34 (MQ=255)                                                    tATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCAGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTgg                                                                                           >  2:495294/1‑150 (MQ=255)                                                     aTTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTggg                                                                                          >  2:244136/1‑150 (MQ=255)                                                        gTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGtt                                                                                                                                          <  1:446433/99‑1 (MQ=255)                                                        gTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGtt                                                                                                                                          >  2:446433/1‑99 (MQ=255)                                                        gTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTAACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATAATATCAGTggggcc                                                                                       >  2:619580/1‑148 (MQ=255)                                                              gCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAg                                                                                 <  2:752980/150‑1 (MQ=255)                                                               cAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCAGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGcat                                                                                                      <  2:96853/128‑1 (MQ=255)                                                               cAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCAGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGcat                                                                                                      >  1:96853/1‑128 (MQ=255)                                                                       aTCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGcatca                                                                                                    <  2:1307133/122‑1 (MQ=255)                                                                       aTCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGcatca                                                                                                    >  1:1307133/1‑122 (MQ=255)                                                                        tCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCtt                                                                       >  2:152964/1‑150 (MQ=255)                                                                            gTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCCTATCAGTGGGGTCAAAAAGCTCCGGCCTTGagc                                                                   >  1:81491/1‑148 (MQ=255)                                                                             tGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTggg                                                                                          <  1:1549821/126‑1 (MQ=255)                                                                             tGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTggg                                                                                          >  2:1549821/1‑126 (MQ=255)                                                                             tGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGaaaa                                                                  <  2:986778/150‑1 (MQ=255)                                                                               aaCAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATggg                                                                                                                        <  1:73996/94‑1 (MQ=255)                                                                               aaCAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATggg                                                                                                                        >  2:73996/1‑94 (MQ=255)                                                                                aCAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATTGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAAccc                                                               >  2:1928877/1‑150 (MQ=255)                                                                                     gTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGgt                                                          <  1:244136/150‑1 (MQ=255)                                                                                          tCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCa                                                     >  1:1132911/1‑150 (MQ=255)                                                                                             gAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGaa                                                  >  2:1162723/1‑150 (MQ=255)                                                                                                      cGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGa                                         <  1:1217802/150‑1 (MQ=255)                                                                                                       gTTCCAACTGGTCACGTTTACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATAAAAAAAAAAAAAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGAt                                        <  2:807339/150‑1 (MQ=255)                                                                                                       gTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGAt                                                                                                                           <  1:252218/67‑1 (MQ=255)                                                                                                       gTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGAt                                                                                                                           >  2:252218/1‑67 (MQ=255)                                                                                                        ttCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATg                                       >  2:583133/1‑150 (MQ=255)                                                                                                         tCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTggg                                                                                          >  1:1652539/1‑98 (MQ=255)                                                                                                         tCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTggg                                                                                          <  2:1652539/98‑1 (MQ=255)                                                                                                            aaCTGGTCACGTTGACTATCCGCCCATCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGCGCAGAAGGCAGTGGATGAAtt                                   >  2:1744985/1‑150 (MQ=255)                                                                                                                       ttGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCtt                                                                       >  1:100984/1‑103 (MQ=255)                                                                                                                       ttGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCtt                                                                       <  2:100984/103‑1 (MQ=255)                                                                                                                       ttGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATGAATTCGACAGTAccc                        <  2:1132911/150‑1 (MQ=255)                                                                                                                         gACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATGAATTCGACAGTACCCTc                      <  1:798313/150‑1 (MQ=255)                                                                                                                              tCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATGAATTCGACAGTACCCTCccagc                 >  2:1430986/1‑150 (MQ=255)                                                                                                                                cGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATGAATTCGACAGTACCCTCccagcca               <  2:1115135/150‑1 (MQ=255)                                                                                                                                cGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATGAATTCGACAGTACCCTCccagcca               >  1:1403039/1‑150 (MQ=255)                                                                                                                                 gCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATGAATTCGACAGTACCCTCccagccag              >  1:811191/1‑150 (MQ=255)                                                                                                                                 gCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATGAATTCGACAGTACCCTCccagccag              >  1:1654495/1‑150 (MQ=255)                                                                                                                                        agcTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCGTGAAAACCCTCGGTGCGCAGAAGGCAGTGGATGAATTCGACAGTACCCTCCCAGCCAGGATaaaa       <  1:1744985/150‑1 (MQ=255)                                                                                                                                           ttCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATGAATTCGACAGTACCCGCGCAGCCAGGATAAAacac    >  2:1116711/1‑150 (MQ=255)                                                                                                                                           ttCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATGAATTAGACAGTACCCTCCCAGCCAGGATAAAacac    <  1:1162723/150‑1 (MQ=255)                                                                                                                                            tCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATGAATTCGACAGTACCCTCCCAGCCAGGATAAAACACg   >  2:419691/1‑150 (MQ=255)                                                                                                                                             cAT‑TGCGGGCGGTTTATTTACCAAGCGATTGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATGAATTCGACAGTACCCTCCCAGCCAGGATAAAACACGc  <  1:1928877/150‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                       TTCAATGCGCCAGACAGCTTGAGCCCGTTCGGTGATGGGGGGGCCAATGCGTATTGTTATTGCAAGGGTGATCCAGTGAACAAAGTGGATCCGAGTGGTCACGTTCCAACTGGTCACGTTGACTATCCGCCCAGCAGCTTCAT‑TGCGGGCGGTTTATTTACCAAGCGATGGGGACTCCGCTCGGCGGCATCATATCAGTGGGGTCAAAAAGCTCCGGCCTTGAAAACCCTCGGTGTGCAGAAGGCAGTGGATGAATTCGACAGTACCCTCCCAGCCAGGATAAAACACGC  >  NC_002947/3581305‑3581594

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