Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,104,697	A→G	6.6%	F63L (TTC→CTC)	PP_1879 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,104,697	0	A	G	6.6%	172.7 / 2.6	60	F63L (TTC→CTC)	PP_1879	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/27);  new base G (2/2);  total (31/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 = 5.44e-01

CAATAAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGC  >  NC_002947/2104550‑2104843                                                                                                                                                    |                                                                                                                                                      cAATAAGGTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATTGGCTTGCTGAAGGTGTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGaa                                                                                                                                                       <  3:308924/148‑1 (MQ=255)    tAAGCTTGTCGACCGGTACCAATATGGATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCAGGGTGCGAATTTGCCAACGGGTGGAACTc                                                                                                                                                    >  2:166591/1‑148 (MQ=255)    tAAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTc                                                                                                                                                    >  7:208409/1‑148 (MQ=255)      aGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCAc                                                                                                                                                 <  7:96802/149‑1 (MQ=255)       gCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCAcc                                                                                                                                                >  4:339539/1‑149 (MQ=255)       gCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCAcc                                                                                                                                                >  2:120595/1‑149 (MQ=255)           gTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGc                                                                                                                                             <  5:145116/148‑1 (MQ=255)            tCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATACCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAg                                                                                                                                           >  7:73945/1‑149 (MQ=255)            tCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGAAAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAGTTCCACCCgtgg                                                                                                                                           <  6:90637/149‑4 (MQ=255)               aCCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAgggg                                                                                                                                        >  7:25024/1‑149 (MQ=255)                ccGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAggggg                                                                                                                                       <  4:109653/149‑1 (MQ=255)                   gTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTACTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGaa                                                                                                                                    <  7:131974/149‑1 (MQ=255)                       cAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAgg                                                                                                                                <  3:339539/149‑1 (MQ=255)                                   ctctTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑Ct                                                                                                                  >  3:255792/1‑149 (MQ=255)                                   ctctTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑Ct                                                                                                                  <  7:223814/149‑1 (MQ=255)                                    tctTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGc                                                                                                                            >  7:204254/1‑140 (MQ=255)                                    tctTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGc                                                                                                                            <  8:204254/140‑1 (MQ=255)                                      ttCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAggggg                                                                                                                                       <  6:21477/127‑1 (MQ=255)                                      ttCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAggggg                                                                                                                                       >  5:21477/1‑127 (MQ=255)                                          ggTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTACGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCgt                                                                                                          >  8:1449/1‑149 (MQ=255)                                               gATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTgg                                                                                                     >  6:321415/1‑149 (MQ=255)                                               gATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATACCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTg                                                                                                      <  8:73945/148‑1 (MQ=255)                                                aTAGGCTTGCTGAATGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTggg                                                                                                    <  8:208409/149‑1 (MQ=255)                                                   ggCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTgggggt                                                                                                 <  4:146487/149‑4 (MQ=255)                                                    gCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTc                                                                                                >  4:115500/1‑149 (MQ=255)                                                     cTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTACTCATCTAATTACGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGCGGCAAGAGGTTGCCCGGAG‑‑AT‑CGCGCGAGTTGGGTTTcc                                                                                               >  4:74731/1‑149 (MQ=255)                                                         gtGAAGGTCTTCTTGTCGACCAGTAAGTCATTCGAGTCTAGGTATTGGTATTGCTCATCTATTTCCGGGTGCGAATTAGCACACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCGCGGAA‑‑CT‑CGAGCGTGTTGGGTTTcccccc                                                                                           <  5:90538/148‑1 (MQ=255)                                                         cTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTcccccc                                                                                           <  3:115500/149‑1 (MQ=255)                                                           gAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCtt                                                                                         >  8:190250/1‑149 (MQ=255)                                                            aaGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTggg                                                                                                    >  6:65746/1‑137 (MQ=255)                                                            aaGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTggg                                                                                                    <  5:65746/137‑1 (MQ=255)                                                                   tcttGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCg                                                                                 >  1:415/1‑149 (MQ=255)                                                                    cttGTCGACCAGTACGTGATTCCATTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGTTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCg                                                                                 <  1:166591/148‑1 (MQ=255)                                                                    cttGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATc                                                                                  <  5:321415/147‑1 (MQ=255)                                                                     ttGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTc                                                                               >  4:365179/1‑149 (MQ=255)                                                                        tCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTccac                                                                            >  1:300907/1‑149 (MQ=255)                                                                               gTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCtt                                                                     >  2:192927/1‑149 (MQ=255)                                                                                   gTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGCCCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCCCGCTTGGTa                                                                 >  4:345920/1‑149 (MQ=255)                                                                                      aTTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACgg                                                              <  7:30343/149‑1 (MQ=255)                                                                                          cAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTgg                                                          <  2:300907/149‑1 (MQ=255)                                                                                             tAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTggggg                                                       <  7:190250/149‑1 (MQ=255)                                                                                             tAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGGTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTggggg                                                       <  4:255792/149‑1 (MQ=255)                                                                                                             agtCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGCTACGGTTGGGGGGACAAGAGAGGCCGAg                                       <  1:210481/147‑1 (MQ=255)                                                                                                               tCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGTAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACAGTGTTGGGTTTCCCCCCTTCAACATCGTCGACAACGCTTGGTACGGTTGGGGGGACAAGAGAGCCCGAGTc                                     >  4:9535/1‑149 (MQ=255)                                                                                                                 aTCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCa                                   >  7:234909/1‑149 (MQ=255)                                                                                                                  tCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAg                                  >  2:203771/1‑149 (MQ=255)                                                                                                                       ttCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCATGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAg                             >  7:182691/1‑149 (MQ=255)                                                                                                                        tCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAg                             <  2:415/148‑1 (MQ=255)                                                                                                                          cGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAggg                           <  1:192927/148‑1 (MQ=255)                                                                                                                             gTGCGAATTTGCCCACGGGTGGAACTCCACCCGCACGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTg                       <  8:182691/149‑1 (MQ=255)                                                                                                                                   aTTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATATGGTGCCCGGAATCCTAAGACCTCATACAAATACCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAg                    <  1:120595/149‑1 (MQ=255)                                                                                                                                      cccctgcGGGTGGAGTTCCACCCGTAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAg                             >  2:208803/7‑134 (MQ=255)                                                                                                                                      cccctgcGGGTGGAGTTCCACCCGTAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAg                             <  1:208803/128‑1 (MQ=255)                                                                                                                                         ctgcGGGTGGAGTTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGt           >  7:299969/4‑149 (MQ=255)                                                                                                                                           aCGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTa         >  7:90083/1‑149 (MQ=255)                                                                                                                                             gggTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAggg                           <  7:115762/129‑1 (MQ=255)                                                                                                                                             gggTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAggg                           >  8:115762/1‑129 (MQ=255)                                                                                                                                                tGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCg      >  6:296039/1‑147 (MQ=255)                                                                                                                                                  gAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAgc  >  1:252649/1‑149 (MQ=255)                                                                                                                                                  gAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAgc  <  2:252649/149‑1 (MQ=255)                                                                                                                                                  gAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAgc  >  3:358369/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CAATAAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCACCCGCAGGGGGGAATAGGTTGCCCGGAA‑‑CT‑CGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCGGAGC  >  NC_002947/2104550‑2104843

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