Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,683,827	A→G	7.2%	G246G (GGT→GGC)	PP_2352 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,683,827	0	A	G	7.2%	157.7 / 2.5	57	G246G (GGT→GGC)	PP_2352	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/26);  new base G (2/2);  total (28/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 = 4.95e-01

TGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCACGTGCATCGAGCTGTTGATCAGGCGAATGACCGGCAGAAAACGTTCGTGATACAGCG  >  NC_002947/2683679‑2683975                                                                                                                                                      |                                                                                                                                                     tGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCa                                                                                                                                                      <  4:48600/149‑1 (MQ=255)    ccGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCAccc                                                                                                                                                   >  7:175257/1‑149 (MQ=255)     cGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCg                                                                                                                                                  >  7:283171/1‑149 (MQ=255)       cAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGc                                                                                                                                                <  2:270474/149‑1 (MQ=255)       cAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGc                                                                                                                                                <  7:163220/149‑1 (MQ=255)            gCCTTCTCGCACAGATCGTCAAAGCCCATCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGa                                                                                                                                            <  6:52510/148‑1 (MQ=255)            gCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGAc                                                                                                                                           <  1:275455/149‑1 (MQ=255)                              tCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTaa                                                                                                                          >  6:167728/1‑148 (MQ=255)                              tCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTaa                                                                                                                          >  7:299490/1‑148 (MQ=255)                                   gCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCt                                                                                                                    <  7:26995/149‑1 (MQ=255)                                     ccAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGTTTGTTCGTTCTGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCTCCCGGCCAGACATAATGGc                                                                                                                                   <  5:224405/132‑1 (MQ=255)                                      cAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGAt                                                                                                                 >  4:148084/1‑149 (MQ=255)                                         caccacCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTgg                                                                                                                              >  5:314976/1‑133 (MQ=255)                                         caccacCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTgg                                                                                                                              <  6:314976/133‑1 (MQ=255)                                            caccCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTcgc                                                                                                           >  2:127201/1‑149 (MQ=255)                                             accCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCt                                                                                                                    >  3:190367/1‑139 (MQ=255)                                             accCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCt                                                                                                                    <  4:190367/139‑1 (MQ=255)                                               ccgccgGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGt                                                                                                        >  6:5813/1‑149 (MQ=255)                                                      ccATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTtggctgg                                                                                                 >  7:188344/1‑149 (MQ=255)                                                              ggCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTg                                                                                          <  8:299490/148‑1 (MQ=255)                                                               gCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGc                                                                                         <  8:132404/148‑1 (MQ=255)                                                                 cAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGa                                                                                       >  4:131076/1‑148 (MQ=255)                                                                         gCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCgt                                                                               <  8:175257/148‑1 (MQ=255)                                                                                    ttGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGAc                                                                   >  2:48942/1‑149 (MQ=255)                                                                                       ttCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCt                                                                                                                    <  3:286765/97‑1 (MQ=255)                                                                                       ttCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCt                                                                                                                    >  4:286765/1‑97 (MQ=255)                                                                                        tCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCGTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATcc                                                                                 <  7:17634/131‑1 (MQ=255)                                                                                        tCACGTCCAGTATCACCGGTTCTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATcc                                                                                 >  8:17634/1‑131 (MQ=255)                                                                                          aCGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGa                                                              <  6:19123/148‑1 (MQ=255)                                                                                           cGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGc                                                            <  5:167728/149‑1 (MQ=255)                                                                                           cGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGc                                                            <  8:188344/149‑1 (MQ=255)                                                                                            gTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCa                                                           >  5:257404/1‑149 (MQ=255)                                                                                              ccAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCACg                                                         >  5:243292/1‑149 (MQ=255)                                                                                              ccAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCACg                                                         >  3:274812/1‑149 (MQ=255)                                                                                                 gTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTc                                                                <  1:172847/139‑1 (MQ=255)                                                                                                 gTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTc                                                                >  2:172847/1‑139 (MQ=255)                                                                                                   aTCACCGGTTGTTCGTTCGGCACGTTCAGCGCCGTGGCGCTGCGCCGCATCGCCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCAc                                                          >  1:197868/1‑144 (MQ=255)                                                                                                   aTCACCGGTTGTTCGTTCGGCACGTTCAGCGCCGTGGCGCTGCGCCGCATCGCCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCAc                                                          >  1:202365/1‑144 (MQ=255)                                                                                                   aTCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCGCCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCAc                                                          <  2:197868/143‑1 (MQ=255)                                                                                                   aTCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCGCCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCAc                                                          <  2:202365/143‑1 (MQ=255)                                                                                                       ccGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACt                                                                  >  4:210699/1‑131 (MQ=255)                                                                                                       ccGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACt                                                                  <  3:210699/131‑1 (MQ=255)                                                                                                           ttgttCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCACGTGCATCGAGCTGt                                            <  7:234112/149‑1 (MQ=255)                                                                                                              ttcgttcgGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCACGTGCATCGAGCTGTTGa                                         >  5:192836/1‑149 (MQ=255)                                                                                                              ttcgttcgGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCACGTGCATCGAGCTGTTGa                                         >  3:87416/1‑149 (MQ=255)                                                                                                                 gttcgGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCACGTGCATCGAGCTGTTGATc                                       <  1:292642/148‑1 (MQ=255)                                                                                                                     ggCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCACGTGCATCGAGCTGTTGATCAGGc                                   <  3:148084/148‑1 (MQ=255)                                                                                                                      gCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCACGTGCATCGAGCTGTTGATCAGGCg                                  >  6:44064/1‑148 (MQ=255)                                                                                                                        aCGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCACGTGCATCGAGCTGTTGATCAGGCGAAt                               <  5:5813/149‑1 (MQ=255)                                                                                                                        aCGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCACGTGCATCGAGCTGTTGATCAGGCGAAt                               <  6:257404/149‑1 (MQ=255)                                                                                                                          gTTCAGCGCC‑CGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCAc                                                                       <  7:209433/107‑1 (MQ=255)                                                                                                                          gTTCAGCGCC‑CGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCAc                                                                       >  8:209433/1‑107 (MQ=255)                                                                                                                                cgcC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCACGTGCATCGAGCTGTTGATCAGGCGAATGACCGGCa                       <  4:87416/149‑1 (MQ=255)                                                                                                                                cgcC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCACGTGCATCGAGCTGTTGATCAGGCGAATGACCGGCa                       <  6:192836/149‑1 (MQ=255)                                                                                                                                 gcC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGATGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCACGTGCATCGAGCTGTTGATCAGGCGAATGACCgg                         >  2:292370/1‑146 (MQ=255)                                                                                                                                     tGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCACGTGCATCGAGCTGTTGATCAGGCGAATGACCGGCGGaa                    >  1:111100/1‑148 (MQ=255)                                                                                                                                     tGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCACGTGCATCGAGCTGTTGATCAGGCGAATGACCGGCAGaa                    >  7:205118/1‑148 (MQ=255)                                                                                                                                             cgccgcATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCACGTGCATCGAGCTGTTGATCAGGCGAATGACCGGCAGAAAACGTTc             <  1:127201/147‑1 (MQ=255)                                                                                                                                                      aCCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCACGTGCATCGAGCTGTTGATCAGGCGAATGACCGGCAGAAAACGTTCGTGATACAGCg  >  6:50336/1‑149 (MQ=255)                                                                                                                                                      aCCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCACGTGCATCGAGCTGTTGATCAGGCGAATGACCGGCAGAAAACGTTCGTGATACAGCg  >  6:304168/1‑149 (MQ=255)                                                                                                                                                      aCCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCACGTGCATCGAGCTGTTGATCAGGCGAATGACCGGCAGAAAACGTTCGTGATACAGCg  >  5:266907/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     TGACCGCAGTGGCCTTCTCGCACAGATCGTCAAAGCCCAGCACCACCCGCCGGCCATCGATGGCCAGCGCTCGCAGCGGGCGCTTGTTCACGTCCAGTATCACCGGTTGTTCGTTCGGCACGTTCAGCGCC‑TGGCGCTGCGCCGCATCACCCGGCCAGACATAATGGCCCTGGGTAAAACGCTGATGCTCGCGGTTGGCTGGCAAACTGCGAAAATCCGTGTCGCCACCGACTTCGAGCACGTGCATCGAGCTGTTGATCAGGCGAATGACCGGCAGAAAACGTTCGTGATACAGCG  >  NC_002947/2683679‑2683975

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