Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,685,015	C→A	6.3%	intergenic (+27/‑112)	PP_3250 → / → PP_5593	transporter/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,685,015	0	C	A	6.3%	200.3 / 3.6	63	intergenic (+27/‑112)	PP_3250/PP_5593	transporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/29);  new base A (2/2);  total (32/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 = 9.00e-01

GATGGGCGAACACCCGCAATTGGCGATGCTCATCGGCGCGTTGCTTGCCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGCCGTGCATTCCGAAGTCATATGAGGTTTGAACACCATGGCCACTAAAGACAACAGCCGCACC  >  NC_002947/3684868‑3685162                                                                                                                                                    |                                                                                                                                                      gATGGGCGAACACCCGCAATTGGCGATGCTCATCGGCGCGTTGCTTGCCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCt                                                                                                                                                      >  1:137517/1‑149 (MQ=255)  aTGGGCGAACACCCGCAATTGGCGATGCTCATCGGCGCGTTGCTTGCCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCtt                                                                                                                                                     <  5:259895/149‑1 (MQ=255)         aaCACCCGCAATTGGCGATGCTCATCGGCGCGTTGCGTGCCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAAc                                                                                                                                              <  8:246158/149‑1 (MQ=255)              ccGCAATTGGCGATGCTCATCGGCGCGTTGCTTGCCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCa                                                                                                                                          <  1:215556/148‑1 (MQ=255)              ccGCAATTGGCGATGCTCATCGGCGCGTTGCTTGCCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCa                                                                                                                                          >  2:215556/1‑148 (MQ=255)                gCAATTGGCGATGCTCATCGGCGCGTTGCTTGCCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCg                                                                                                                                       <  5:208888/149‑1 (MQ=255)                   aTTGGCGATGCTCATCGGCGCGTTGCTTGCCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTc                                                                                                                                    <  2:304888/149‑1 (MQ=255)                   aTTGGCGATGCTCATCGGCGCGTTGCTTGCCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTc                                                                                                                                    <  7:145584/149‑1 (MQ=255)                   aTTGGCGATGCGCATCGGCGCGTTGCTTGCCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTc                                                                                                                                    <  8:34036/149‑1 (MQ=255)                      ggCGATGCTCATCGGCGCGTTGCTTGCCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCg                                                                                                                                  >  1:258870/1‑148 (MQ=255)                        cGATGCTCATCGGCGCGTTGCTTGCCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAAc                                                                                                                               <  2:137517/149‑1 (MQ=255)                            gCTCATCGGCGCGTTGCTTGCCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACctct                                                                                                                           >  7:134407/1‑149 (MQ=255)                              tCATCGGCGCGTTGCTTGCCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCg                                                                                                                         >  1:24768/1‑149 (MQ=255)                              tCATCGGCGCGTTGCTTGCCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATAAGGCGATGCCGCTTATTTAACGCGACCGGTCCGCTCCtctc                                                                                                                          >  7:254383/1‑148 (MQ=255)                                aTCGGCGCGTTGCTTGCCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGc                                                                                                                                            >  4:259681/1‑128 (MQ=255)                                aTCGGCGCGTTGCTTGCCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGc                                                                                                                                            <  3:259681/128‑1 (MQ=255)                                aTCGGCGCGTTGCTTGCCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGt                                                                                                                        >  3:135772/1‑148 (MQ=255)                                 tCGGCGCGTTGCTTGCCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGtaa                                                                                                                      <  3:286889/149‑1 (MQ=255)                                   ggCGCGTTGCGTGCCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGtaata                                                                                                                    <  2:223406/149‑1 (MQ=255)                                         ttgcttgcCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCGACCGGTCCGAACCTCTCGTAATAAAGCGa                                                                                                              <  1:303725/149‑1 (MQ=255)                                         ttgcttgcCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGa                                                                                                              <  6:35263/149‑1 (MQ=255)                                             ttgcCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCt                                                                                                        >  3:42786/1‑149 (MQ=255)                                                ccGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGcc                                                                                                     <  4:135772/149‑1 (MQ=255)                                                   tGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGat                                                                                                  >  2:134025/1‑149 (MQ=255)                                                   tGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGa                                                                                                   >  2:94783/1‑148 (MQ=255)                                                                   aGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAAc                                                                                                                                              >  1:253956/1‑91 (MQ=255)                                                                   aGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAAc                                                                                                                                              <  2:253956/91‑1 (MQ=255)                                                                     gcgcTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATc                                                                                >  3:204045/1‑149 (MQ=255)                                                                     gcgcTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATc                                                                                >  7:98745/1‑149 (MQ=255)                                                                      agcgtGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAATCGA‑‑ATCTGCCTATATCGCCGAACTCCGAATCt                                                                               <  1:178982/145‑1 (MQ=255)                                                                       gcTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCtg                                                                              >  8:122154/1‑149 (MQ=255)                                                                           ggCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTgtgt                                                                           <  6:44277/148‑1 (MQ=255)                                                                            gCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGAAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTAc                                                                         <  1:232320/149‑1 (MQ=255)                                                                              tGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAAc                                                                                        <  4:139315/132‑1 (MQ=255)                                                                              tGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAAc                                                                                        >  3:139315/1‑132 (MQ=255)                                                                               gggCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACgg                                                                       >  2:115673/1‑148 (MQ=255)                                                                                   cgccgGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCt                                                                  >  1:305220/1‑149 (MQ=255)                                                                                       ggCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGATTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGcc                                                              >  2:321581/1‑149 (MQ=255)                                                                                       ggCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGATTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGcc                                                              >  2:326620/1‑149 (MQ=255)                                                                                              gctcATGTATAGCACGCCGAACCACTGACACTCTCACCAAATATGGAGATGCCGCTTATTTAACTCCACCGGTAAGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCTCCGATCTAAGAATCTGTTTACGGCATCTTGCCTTGCAtt                                                       <  8:254383/146‑1 (MQ=255)                                                                                                   tGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGCCGTGCATTCCGaa                                                  <  4:204045/149‑1 (MQ=255)                                                                                                   tGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGCCGTGCATTCCGaa                                                  <  8:134407/149‑1 (MQ=255)                                                                                                   tGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGCCGTGCATTCCGaa                                                  <  2:305220/149‑1 (MQ=255)                                                                                                       tAGCCCGCCGCACCACTGACACGCGCACCAAAAATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGATCATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGCCGTGCATTCCGAAGt                                                <  1:94783/149‑1 (MQ=255)                                                                                                          cccGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGCCGTGCATTCCGAAGTCatat                                           >  5:307567/1‑149 (MQ=255)                                                                                                                gcACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAAc                                                                                        >  2:322251/1‑98 (MQ=255)                                                                                                                gcACCACTGACACGCGCACCAAATATGGCGATGCCGATTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGCCGTGCATTCCGAAGTCATATGAGGtt                                     <  1:321581/149‑1 (MQ=255)                                                                                                                gcACCACTGACACGCGCACCAAATATGGCGATGCCGATTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGCCGTGCATTCCGAAGTCATATGAGGtt                                     <  1:326620/149‑1 (MQ=255)                                                                                                                    cacTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGCCGTGCATTCCGAAGTCATATGAGGTTTGac                                 >  4:267507/1‑148 (MQ=255)                                                                                                                         acacGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGc                                                                                                      <  7:30908/75‑1 (MQ=255)                                                                                                                         acacGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGc                                                                                                      >  8:30908/1‑75 (MQ=255)                                                                                                                         acacGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGCCGTGCATTCCGAAGTCATATGAGGTTTGAACACCa                            <  6:307567/149‑1 (MQ=255)                                                                                                                         acacGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGCCGTGCATTCCGAAGTCATATGAGGTTTGAACACCa                            >  7:222932/1‑149 (MQ=255)                                                                                                                            cgcgCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGCCGTGCATTCCGAAGTCATATGAGGTTTGAACACCATgg                         <  7:122154/149‑1 (MQ=255)                                                                                                                             gcgcACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGCACTCCGAATCTGTGTACGGCATCTTGCCGTGCATTCCGAAGTCATATGAGGTTTGAACACCATGGc                        >  2:193690/1‑149 (MQ=255)                                                                                                                               gcACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGCCGTGCATTCCGAAGTCATATGAGGTTTGAACACCATGGCCa                      <  2:24768/149‑1 (MQ=255)                                                                                                                                  ccAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGCAGTGCATTCCGAAGTCATATGAGGTTTGAACACCATGGCCACt                    >  8:256370/1‑148 (MQ=255)                                                                                                                                   cAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGCCGTGCATTCCGAAGTCATATGAGGTTTGAACACCATGGCCACTa                   >  2:69092/1‑148 (MQ=255)                                                                                                                                    aaaTATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGCCGTGCATTCCGAAGTCATATGAGGTTTGAACACCATGGCCACTa                   >  1:169174/1‑147 (MQ=255)                                                                                                                                        atGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGCCGTGCATTCCGAAGTCATATGAGGTTTg                                   >  1:300894/1‑127 (MQ=255)                                                                                                                                        atGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGCCGTGCATTCCGAAGTCATATGAGGTTTg                                   <  2:300894/127‑1 (MQ=255)                                                                                                                                                 ccGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGCCGTGCATTCCGAAGTCATATGAGGTTTGAACACCATGGCCACTAAAGACAACAGCCGCa    >  7:51326/1‑149 (MQ=255)                                                                                                                                                   gCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGCCGTGCATTCCGAAGTCATATGAGGTTTGAACACCATGGCCACTAAAGACAACAGCCGCAcc  >  3:104705/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GATGGGCGAACACCCGCAATTGGCGATGCTCATCGGCGCGTTGCTTGCCGTGGGCTGGGGCGCGACAGGCGCTTGGCTGGGCCGCCGGCAGGCGCGCATGGATAGCCCGCCGCACCACTGACACGCGCACCAAATATGGCGATGCCGCTTATTTAACGCCACCGGTCCGAACCTCTCGTAATAAAGCGA‑‑ATCTGCCGATATCGCCGAACTCCGAATCTGTGTACGGCATCTTGCCGTGCATTCCGAAGTCATATGAGGTTTGAACACCATGGCCACTAAAGACAACAGCCGCACC  >  NC_002947/3684868‑3685162

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