Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,350,476	G→C	5.6%	intergenic (+13/+69)	PP_4704 → / ← PP_4705	hypothetical protein/glyoxalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,350,476	0	G	C	5.6%	225.6 / 3.1	71	intergenic (+13/+69)	PP_4704/PP_4705	hypothetical protein/glyoxalase
Reads supporting (aligned to +/- strand):  ref base G (38/29);  new base C (2/2);  total (40/31)
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.77e-01

CCAGCAGCATCAAGGACTATCTGGTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGG  >  NC_002947/5350330‑5350617                                                                                                                                                   |                                                                                                                                               ccAGCAGCATCAAGGACTATCTGGTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGGCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGgc                                                                                                                                              >  8:173500/1‑149 (MQ=255)      agcaTCAAGGACTATCTGGTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGgcccgcc                                                                                                                                         >  3:289136/1‑149 (MQ=255)          acAAGGACTATCTGGTCGACAAGGGCAAGCTGGAAGACGACAGGGTGTAGTTCATCGATACCAACCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAAGCCCCTCTGGCAGgcccgcccgcc                                                                                                                                     >  8:338727/2‑146 (MQ=255)           cAAGGACTATCTGGTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCAGCCCCca                                                                                                                                     >  3:98237/1‑148 (MQ=255)                 cTATCTGGTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCt                                                                                                                             >  7:55108/1‑149 (MQ=255)                    tCTGGTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTgcg                                                                                                                          >  1:184516/1‑149 (MQ=255)                       ggTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAgac                                                                                                                                >  5:254998/1‑139 (MQ=255)                       ggTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAgac                                                                                                                                <  6:254998/141‑3 (MQ=255)                        gTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGGCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTgcg                                                                                                                          <  7:173500/145‑1 (MQ=255)                        gTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTgcg                                                                                                                          <  1:211596/145‑1 (MQ=255)                        gTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTgcg                                                                                                                          <  6:180824/145‑1 (MQ=255)                        gTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTgcg                                                                                                                          <  7:25978/145‑1 (MQ=255)                           gACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTAcc                                                                                                                   >  3:44293/1‑149 (MQ=255)                              aaGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTCCGCACTACCtgt                                                                                                                >  1:242857/1‑149 (MQ=255)                                gggCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTgg                                                                                                              <  4:44293/149‑1 (MQ=255)                                   cAAGCTGGAAGACGACCGGGTGTACTTCATCGATACAAGCCTGGGGAAGCATCAAAGCGATGGCCGGGTGGTGACACCGATGCACATGGATGCCGAGTAATCCCCTCTGGCAGGCACGCCCACACAGG‑CCTGCGCACTACCTGTGGGAg                                                                                                           >  5:331945/1‑149 (MQ=255)                                            aGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTCTGGGAGCGGCCGTGc                                                                                                  >  5:390077/1‑149 (MQ=255)                                            aGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCACGCCCGCTCCCACAGGTAGTGCGCA‑GGCCTGTGGGGGCGGGCCTGc                                                                                                  <  1:248106/149‑1 (MQ=25)                                              acgacCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGccc                                                                                                >  5:236278/1‑149 (MQ=255)                                                   cGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                           <  4:229713/149‑1 (MQ=255)                                                     ggTGTACTTCATCGATACCAGCCTGGGGCAGACTCAAATCGATGTCCGGGTGGTGACAACGATGCACCTGGTTGCCGAGTAATCACCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGTGCGTGCCCGCGAaga                                                                                         <  6:331945/149‑1 (MQ=255)                                                       tgtACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAgg                                                                                       <  7:274360/149‑1 (MQ=255)                                                         tACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGc                                                                                      >  7:160108/1‑148 (MQ=255)                                                         tACTGCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGc                                                                                      <  4:52321/148‑1 (MQ=255)                                                            ttCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGT‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAg                                                                                   >  8:171845/1‑148 (MQ=255)                                                             tCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGCTGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGACGAGGCCCGcc                                                                                 >  6:361551/1‑148 (MQ=255)                                                             tCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCa                                                                                 >  7:3074/1‑149 (MQ=255)                                                              cATCGATACCATCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACTCCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCa                                                                                 <  6:242929/148‑1 (MQ=255)                                                               aTCGATACCAGCCTGGGGCAGCCTAAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGAATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACg                                                                               >  7:59333/1‑149 (MQ=255)                                                                   aTACCAGCCTGGGGCAGCCTCAAAGCGATGGTCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCCGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCt                                                                           >  6:207593/1‑149 (MQ=255)                                                                   aTACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCt                                                                           >  1:79292/1‑149 (MQ=255)                                                                       cAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCa                                                                        >  6:280849/1‑148 (MQ=255)                                                                       cAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCa                                                                        >  8:97512/1‑148 (MQ=255)                                                                          ccTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAaac                                                                     >  3:206653/1‑148 (MQ=255)                                                                            tGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAaaccaa                                                                  >  4:358330/1‑149 (MQ=255)                                                                              gggCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAaaccaa                                                                  <  3:358330/147‑1 (MQ=255)                                                                                     cTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGccc                                                         <  7:142812/149‑1 (MQ=255)                                                                                      tCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGCCGTGCCCGCCAAGAGCCCAGCACGGGCTTCAAACCAACCCCTCcccc                                                        >  5:310727/1‑149 (MQ=255)                                                                                      tCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCACGCCCGCTCCCACAGGTAGTGCGCA‑GGCCTGTGGgggcgggcct                                                                                                    <  8:139308/105‑3 (MQ=17)                                                                                      tCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCACGCCCGCTCCCACAGGTAGTGCGCA‑GGCCTGTGGgggcgggcct                                                                                                    >  7:139308/1‑103 (MQ=17)                                                                                       cAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGAGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGATGcccc                                                        >  5:246960/1‑148 (MQ=255)                                                                                        aaaGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGcccc                                                        <  2:79292/147‑1 (MQ=255)                                                                                        aaaGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCAGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGa                                                      >  7:254422/1‑149 (MQ=255)                                                                                        aaaGCGATGGCCGGGGGGTGACACCGATGCAACTGGATTCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGcccc                                                        <  6:310727/147‑1 (MQ=255)                                                                                          aGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGcca                                                    >  7:306763/1‑149 (MQ=255)                                                                                              aTGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACAccgc                                                <  2:184516/149‑1 (MQ=255)                                                                                                 gCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCgg                                             <  6:236278/149‑1 (MQ=255)                                                                                                    gggTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGTCGGGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTg                                          <  5:361551/149‑1 (MQ=255)                                                                                                          tgACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGt                                           <  3:75838/142‑1 (MQ=255)                                                                                                          tgACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGt                                           >  4:75838/1‑142 (MQ=255)                                                                                                            acacCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGc                                   <  8:3074/148‑1 (MQ=255)                                                                                                             cacCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCgg                                 <  4:206653/149‑1 (MQ=255)                                                                                                              acCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCGCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGt                                <  7:171845/149‑1 (MQ=255)                                                                                                              acCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGAt                                        <  5:67954/141‑1 (MQ=255)                                                                                                              acCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGAt                                        >  6:67954/1‑141 (MQ=255)                                                                                                              acCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGt                                >  1:304050/1‑149 (MQ=255)                                                                                                              acCGATGCACCTGGATGCAGAGTAATCCCCTCTTGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTGAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGt                                >  8:287502/1‑149 (MQ=255)                                                                                                               ccGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGtt                               <  2:304050/149‑1 (MQ=255)                                                                                                                   tGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAAc                           <  7:97512/149‑1 (MQ=255)                                                                                                                      aCCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTAGCTGTGGGAGAGGGCGTGCCCGCTAATAGGCCAGCAGGGGCTTCAAACCAACCGCTGCCCCGACACCGACGGGTGATACAGCTGTTGAACCtt                        <  2:242857/149‑1 (MQ=255)                                                                                                                      aCCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCtt                        >  7:173598/1‑149 (MQ=255)                                                                                                                          ggATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTATGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTc                    >  4:327795/1‑149 (MQ=255)                                                                                                                               ccGAGTAATCCCCTCTGGCAGGCGCGGCCCCACAGG‑AGTGCGCACTAGCTGTGGGAGCTGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAATCCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTGccgg                <  8:306763/148‑1 (MQ=255)                                                                                                                                  aGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCccggccgg            >  5:37679/1‑149 (MQ=255)                                                                                                                                       tCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTc         >  6:141851/1‑147 (MQ=255)                                                                                                                                       tCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTc         >  2:135618/1‑147 (MQ=255)                                                                                                                                        ccccGCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCcag      <  4:289136/149‑1 (MQ=255)                                                                                                                                          ccTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCcag      >  1:42766/1‑147 (MQ=255)                                                                                                                                            gcgggCACGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAgg  >  4:317628/4‑149 (MQ=255)                                                                                                                                                  aGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGc                                                                                                      <  8:143100/43‑1 (MQ=255)                                                                                                                                                  aGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGc                                                                                                      >  7:143100/1‑43 (MQ=255)                                                                                                                                                   |                                                                                                                                               CCAGCAGCATCAAGGACTATCTGGTCGACAAGGGCAAGCTGGAAGACGACCGGGTGTACTTCATCGATACCAGCCTGGGGCAGCCTCAAAGCGATGGCCGGGTGGTGACACCGATGCACCTGGATGCCGAGTAATCCCCTCTGGCAGGCCCGCCCCCACAGG‑CCTGCGCACTACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGCCAGCACGGGCTTCAAACCAACCGCTGCCCCGACACCGCCGGGTGATACAGCGGTTGAACCTTGTTCCCGGCCGGGTCCAGCAGG  >  NC_002947/5350330‑5350617

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