Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,364,437	A→T	6.0%	S278T (TCG→ACG)	PP_2961 ←	LysR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,364,437	0	A	T	6.0%	202.6 / 4.3	67	S278T (TCG→ACG)	PP_2961	LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (37/26);  new base T (2/2);  total (39/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 = 4.69e-01

TGCGAGCCCGAAGCAATGGCTGCAGCAGGGAACGGTGTGACGCCTGGTCGACCGACTTGCACGCCAGGGAATGAGCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCAAACAGCCACTGCGGAGACAGGGAGATGCCCAGGCCTCCGAGGACCAT  >  NC_002947/3364288‑3364584                                                                                                                                                      |                                                                                                                                                    tgCGAGCCCGAAGCAATGGCTGCAGCAGGGAACGGTGTGACGCCTGGTCGACCGACTTGCACGCCAGGGAATGAGCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGCACCAGAGTGATGAACGCCATCGTCCCTGCCGa                                                                                                                                                     <  1:1533209/150‑1 (MQ=255)       cccGAAGCAATGGCTGCAGCAGGGAACGGTGTGACGCCTGGTCGACCGACTTGCACGCCAGGGAATGAGCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAgcggcg                                                                                                                                               >  1:234647/1‑150 (MQ=255)          gAAGCAATGGCTGCAGCAGGGAACGGTGTGACGCCTGGTCGACCGACTTGCACGCCAGGGAATGAGCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGcgccg                                                                                                                                            >  1:1555144/1‑150 (MQ=255)               aaTGGCTGCAGCAGGGAACGGTGTGACGCCTGGTCGACCGACTTGCACGCCAGGGAATGAGCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCAc                                                                                                                                       >  1:974539/1‑150 (MQ=255)                aTGGCCGCAGCAGGGGACGATGTGACGCCTGGTCGACCGACGTGCAGGCCAGGGAATGAGCAAGGGCCGGCGCCTCATCGGGGTGGCCTACACGTTCTCCTGAAGAAAATCGATGAACGCCATCGGCCTTGCCGAGCGGCGCCGCTCaag                                                                                                                                      >  2:480789/1‑148 (MQ=255)                              gAACGGTGTGACGCCTGGTCGACCGACTTGCACGCCAGGGAATGAGCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGgcgc                                                                                                                        <  2:443915/150‑1 (MQ=255)                                 cGGTGTGACGCCTGGTCGACCGACTTGCACGCCAGGGAATGAGCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGa                                                                                                                     <  1:581197/150‑1 (MQ=255)                                               gTCGACCGACTTGCACGCCAGGGAATGAGCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGgtc                                                                                                       >  1:549725/1‑150 (MQ=255)                                                  gaccgacTTGCACGCCAGGGAATGAGCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGgtcgtc                                                                                                    <  1:480789/150‑1 (MQ=255)                                                     cgacTTGCACGCCAGGGAATGAGCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTgg                                                                                                 >  2:531599/1‑150 (MQ=255)                                                       acTTGCACGCCAGGGAATGAGCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTa                                                                                               >  1:1018620/1‑150 (MQ=255)                                                          tGCACGCCAGGGAATGAGCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATtc                                                                                            <  2:234647/150‑1 (MQ=255)                                                             aCGCCAGGGAATGAGCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATtcgtc                                                                                         <  2:798050/150‑1 (MQ=255)                                                               gCCAGGGAATGAGCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGa                                                                                       >  1:1182566/1‑150 (MQ=255)                                                                 cAGGGAATGAGCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACTGGCAGCGCGGTCGTCTGGTATTCGTCGAGc                                                                                     <  1:501966/150‑1 (MQ=255)                                                                 cAGGGAATGAGCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGc                                                                                     >  2:702853/1‑150 (MQ=255)                                                                 cAGGGAATGAGCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGc                                                                                     >  1:1086800/1‑150 (MQ=255)                                                                         gAGCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTaca                                                                             <  2:549725/150‑1 (MQ=255)                                                                           gCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGAc                                                                                                                    <  1:96626/109‑1 (MQ=255)                                                                           gCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGAc                                                                                                                    >  2:96626/1‑109 (MQ=255)                                                                           gCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATtcgtc                                                                                         >  1:1407745/1‑136 (MQ=255)                                                                           gCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATtcgtc                                                                                         <  2:1407745/136‑1 (MQ=255)                                                                            cAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTacacac                                                                          <  2:1018620/150‑1 (MQ=255)                                                                            aaaGGTCCTGCTCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTGACGGCTGAACACGCCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTacacac                                                                          <  2:531345/149‑1 (MQ=255)                                                                               ggTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGcc                                                                       >  2:192740/1‑150 (MQ=255)                                                                                            aTCGGGGTGGCCTACACGTTCTCCCGAATAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAg                                                          >  1:1278670/1‑150 (MQ=255)                                                                                            aTCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGTGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAg                                                          >  2:45553/1‑150 (MQ=255)                                                                                                  gTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGTGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCgca                                                    <  1:45553/150‑1 (MQ=255)                                                                                                       cTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAgc                                                                                                                                                   >  1:1186433/1‑50 (MQ=255)                                                                                                       cTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAgc                                                                                                                                                   <  2:1186433/50‑1 (MQ=255)                                                                                                       cTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCaaa                                               >  2:1385563/1‑150 (MQ=255)                                                                                                       cTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCaaa                                               >  1:17455/1‑150 (MQ=255)                                                                                                       cTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCaaa                                               >  2:1274544/1‑150 (MQ=255)                                                                                                         acacGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAgg                                                                                  <  1:1490486/113‑1 (MQ=255)                                                                                                         acacGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAgg                                                                                  >  2:1490486/1‑113 (MQ=255)                                                                                                           acGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCAAACAGc                                           <  2:974539/150‑1 (MQ=255)                                                                                                             gTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGAc                                                                                                                    <  2:393331/75‑1 (MQ=255)                                                                                                             gTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGAc                                                                                                                    >  1:393331/1‑75 (MQ=255)                                                                                                             gTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAgc                                                                                                             >  1:1320047/1‑82 (MQ=255)                                                                                                             gTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAgc                                                                                                             <  2:1320047/82‑1 (MQ=255)                                                                                                             gTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCCGCTCGaa                                                           <  2:1011124/132‑1 (MQ=255)                                                                                                             gTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCCGCTCGaa                                                           >  1:1011124/1‑132 (MQ=255)                                                                                                                   ccGAAGAAAATCGATGAACGCCATCGTCCTTGCCGTGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAg                                                                                                              >  2:471628/1‑75 (MQ=255)                                                                                                                   ccGAAGAAAATCGATGAACGCCATCGTCCTTGCCGTGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAg                                                                                                              <  1:471628/75‑1 (MQ=255)                                                                                                                      aagaaAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCAAACAGCCACTGCGgaga                                >  2:1566782/1‑150 (MQ=255)                                                                                                                         aaaaTCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCAAACAGCCACTGCGGAGACAg                             <  2:1278670/150‑1 (MQ=255)                                                                                                                           aaTCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCAAACa                                             >  2:935680/1‑132 (MQ=255)                                                                                                                           aaTCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTATGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCAAACa                                             <  1:935680/132‑1 (MQ=255)                                                                                                                             tCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCAAACAGCCACTGCGGAGACAGGgag                         <  1:1566782/150‑1 (MQ=255)                                                                                                                             tCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCAAACAGCCACTGCGGAGACAGGgag                         >  2:343493/1‑150 (MQ=255)                                                                                                                             tCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCAAACAGCCACTGCGGAGACAGGgag                         >  2:1118488/1‑150 (MQ=255)                                                                                                                               gATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCCCGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCAAACAGCCACTGCGGAGACAGGGAGAt                       <  1:343493/150‑1 (MQ=255)                                                                                                                                 tGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCAAACAGCCACTGCGGAGACAGGGAGATGc                     <  2:1086800/150‑1 (MQ=255)                                                                                                                                      gCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCAAACAGCCACTGCGGAGACAGGGAGATGCCCAgg                >  1:435380/1‑150 (MQ=255)                                                                                                                                        cATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCAAACAGCCACTGCgg                                   <  1:40886/129‑1 (MQ=255)                                                                                                                                        cATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCAAACAGCCACTGCgg                                   >  2:40886/1‑129 (MQ=255)                                                                                                                                        cATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCAAACAGCCACTGCGGAGACAGGGAGATGCCCAGGcc              >  1:525408/1‑150 (MQ=255)                                                                                                                                         aTCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCTAAGTCCGCAGCAAACAGCCACTGCGGAGACAGGGAGATGCCCAGGccc             >  2:1486375/1‑149 (MQ=255)                                                                                                                                          tCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTGCGCAGCAAACAGCCACTGCGGAGACAGGGAGATGCCCAGGCCTc            >  2:1184980/1‑150 (MQ=255)                                                                                                                                             tCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCAAACAGCCACTGCGGAGACAGGGAGATGCCCAGGCCTCCGa         >  2:399428/1‑150 (MQ=255)                                                                                                                                                ttGCCGAGCGGCGCCGCTCACGGCTGAACACAGCGCTGACCGACAGTGCGGTCGTCTGGTATTAGTCGAGCAGGCTACACACGCCACCCTGCTCGACGTCCGCAGCAAACAGCCACTGCGGAGACAGGTAGATGCCCAGGCCTCCGAGGa      >  1:220996/1‑150 (MQ=255)                                                                                                                                                  gCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCAAACAGCCACTGCGGAGACAGGGAGATGCCCAGGCCTCCGAGGAcc    <  1:1274544/150‑1 (MQ=255)                                                                                                                                                  gCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCAAACAGCCACTGCGGAGACAGGGAGATGCCCAGGCCTCCGAGGAc     >  1:1146449/1‑149 (MQ=255)                                                                                                                                                  gCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCAAACAGCCACTGCGGAGACAGGGAGATGCCCAGGCCTCCGAGGAc     <  2:1146449/149‑1 (MQ=255)                                                                                                                                                   ccGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCAAACAGCCACTGCGGAGACAGGGAGATGCCCAGGCCTCCGAGGACCa   >  1:1050930/1‑150 (MQ=255)                                                                                                                                                    cGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCAAACAGCCACTGCGGAGACAGGGAGATGCCCAGGCCTCCGAGGACCAt  >  2:450557/1‑150 (MQ=255)                                                                                                                                                     gAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAg                                                          <  1:219192/93‑1 (MQ=255)                                                                                                                                                     gAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAg                                                          >  2:219192/1‑93 (MQ=255)                                                                                                                                                      |                                                                                                                                                    TGCGAGCCCGAAGCAATGGCTGCAGCAGGGAACGGTGTGACGCCTGGTCGACCGACTTGCACGCCAGGGAATGAGCAAGGTCCTGCGCCTCATCGGGGTGGCCTACACGTTCTCCCGAAGAAAATCGATGAACGCCATCGTCCTTGCCGAGCGGCGCCGCTCACGGCTGAACACGGCGCTGACCGGCAGCGCGGTCGTCTGGTATTCGTCGAGCAGGCTACACACGCCACCCTGCTCGAAGTCCGCAGCAAACAGCCACTGCGGAGACAGGGAGATGCCCAGGCCTCCGAGGACCAT  >  NC_002947/3364288‑3364584

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