Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	23,510	C→T	6.4%	A112V (GCA→GTA)	PP_0018 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	23,510	0	C	T	6.4%	180.8 / 2.1	63	A112V (GCA→GTA)	PP_0018	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/30);  new base T (2/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-01

GAATTCTTGGGCGTTCCTTAATCGGGATGCTGACTACGCGTGCGTGCTGGCCATGGATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTAGCCGACTACCTTCCGGAAGCCATGCTCCCACTTAAATGCGATGAGCTAGCGGCGCATGGGAAAGTCCCCAAGGAGC  >  NC_002947/23364‑23659                                                                                                                                                   |                                                                                                                                                      gAATTCTTGGGCGTTCCTTAATCGGGATGCTGACTACGCGTGCGTGCTGGCCATGGATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAg                                                                                                                                                     >  8:170151/1‑149 (MQ=255)       ttGGGCGTTCCTTAATCGGGATGCTGACTACGCGTGCGTGCTGGCCATGGATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCAt                                                                                                                                               >  6:42437/1‑149 (MQ=255)          ggCGTTCCTTAATCGGGATGCTGACTACGCGTGCGTGCTGGCCATGGATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTgg                                                                                                                                            <  7:170151/149‑1 (MQ=255)          ggCGTTCCTTAATCGGGATGCTGACTACGCGTGCGTGCTGGCCATGGATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGAAAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTg                                                                                                                                             >  2:353724/1‑148 (MQ=255)            cGTTCCTTAATCGGGATGCTGACTACGCGTGCGTGCTGGCCATGGATCACATGATCCTTCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCATGTGCATTGGat                                                                                                                                          <  4:125676/149‑1 (MQ=255)                    aaTCGGGATGCTGACTACGCGTGCGTGCTGGCCATGGATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCAcc                                                                                                                                    <  4:99456/147‑1 (MQ=255)                     aTCGGGATGCTGACTACGCGTGCGTGCTGGCCATGGATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGc                                                                                                                                  <  5:160242/148‑1 (MQ=255)                             gCTGACTACGCGTGCGTGCTGGCCATGGATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTCAGTGAATTGAGGCGCTGGTAGGTCAAGCAGGTGAATTGGATATCACCGCTAGCGGGAt                                                                                                                         >  8:353824/1‑149 (MQ=255)                               tGACTACGCGTGCGTGCTGGCCATGGATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGc                                                                                                                       >  3:272899/1‑149 (MQ=255)                                  cTACGCGTGCGTGCTGGCCATGGATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGt                                                                                                                                                   <  1:154273/118‑1 (MQ=255)                                  cTACGCGTGCGTGCTGGCCATGGATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGt                                                                                                                                                   >  2:154273/1‑118 (MQ=255)                                   tACGCGTGCGTGCTGGCCATGGATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTa                                                                                                                   >  2:267085/1‑149 (MQ=255)                                     cgcgTGCGTGCTGGCCATGGATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACg                                                                                                                 >  4:233537/1‑149 (MQ=255)                                      gcgtgcgtGCTGGCCATGGATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACg                                                                                                                 >  8:131590/1‑148 (MQ=255)                                      gcgtgcgtGCTGGCCATGGATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGt                                                                                                                <  7:185149/149‑1 (MQ=255)                                           cgtgcTGGCCATGGATCACATGATTCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCg                                                                                                                      >  7:235514/1‑138 (MQ=255)                                           cgtgcTGGCCATGGATCACATGATTCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCg                                                                                                                      <  8:235514/138‑1 (MQ=255)                                           cgtgcTGGCCATGGATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGTAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGAt                                                                                                            >  3:258612/1‑148 (MQ=255)                                           cgtgcTGGCCATGGATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGTAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGAt                                                                                                            >  5:134885/1‑148 (MQ=255)                                             tgcTGGCCATGGATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGt                                                                                                                <  7:105138/142‑1 (MQ=255)                                             tgcTGGCCATGGATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGt                                                                                                                >  8:105138/1‑142 (MQ=255)                                                    cATGGATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAAc                                                                                                  <  4:272286/149‑1 (MQ=255)                                                    cATGGATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAAc                                                                                                  >  4:299364/1‑149 (MQ=255)                                                       ggATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAg                                                                                               >  2:301340/1‑149 (MQ=255)                                                       ggATAACATGATCCTGCAAAGCACAAATTATCTTCATGTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAg                                                                                               >  6:163116/1‑149 (MQ=255)                                                        gATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAg                                                                                               <  4:165049/148‑1 (MQ=255)                                                         aTCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGGAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCt                                                                                             <  4:272899/149‑1 (MQ=255)                                                         aTCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCt                                                                                             >  3:13026/1‑149 (MQ=255)                                                             caTGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGaa                                                                                                   <  8:269789/139‑1 (MQ=255)                                                             caTGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGaa                                                                                                   >  7:269789/1‑139 (MQ=255)                                                               tGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGt                                                                                       <  4:353671/149‑1 (MQ=255)                                                                 aTCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTg                                                                                      <  7:131590/148‑1 (MQ=255)                                                                        aaaGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTa                                                                              <  7:323334/149‑1 (MQ=255)                                                                            cacaAATTATCTTCATTTTGATAACCGCTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTAGcc                                                                           >  3:237449/1‑148 (MQ=255)                                                                              caAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTAGCCGAc                                                                        <  3:233537/149‑1 (MQ=255)                                                                                 aTTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGa                                                                                                    >  6:325174/1‑118 (MQ=255)                                                                                 aTTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGa                                                                                                    <  5:325174/118‑1 (MQ=255)                                                                                    aTCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTg                                                                                      >  7:314377/1‑129 (MQ=255)                                                                                    aTCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTg                                                                                      >  7:308977/1‑129 (MQ=255)                                                                                    aTCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTg                                                                                      <  8:308977/129‑1 (MQ=255)                                                                                    aTCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTg                                                                                      <  8:314377/129‑1 (MQ=255)                                                                                         cATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTAGCCGACTACCTTc                                                                 <  3:299364/145‑1 (MQ=255)                                                                                          aTTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTAGCCGACTACCTTCCGGaa                                                            <  4:13026/149‑1 (MQ=255)                                                                                          aTTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGAACCTAGCAGACTACCTTCCgg                                                              >  4:237850/1‑147 (MQ=255)                                                                                                aTAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAgg                                                                                        >  4:237934/1‑115 (MQ=255)                                                                                                aTAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAgg                                                                                        <  3:237934/115‑1 (MQ=255)                                                                                                    ccACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATCTCACCGCTATCGGTATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTAGCCGACTACCTTCCGGAAGCCATGCTcc                                                  <  5:163116/149‑1 (MQ=255)                                                                                                    ccACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTAGCCGACTACCTTCCGGAAGCCATGCTcc                                                  <  1:301340/149‑1 (MQ=255)                                                                                                          aGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAATATTGAACCAGCTGAAGGTGAGACCCTAGCCGACTACCTTCCGGAAGCCATGCTCCCACtt                                             >  8:2389/1‑148 (MQ=255)                                                                                                            cgcTTTTCCGTGATTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAAGCAGCTGAAGGTGAGACCCTAGCCGACTACCTTCCGGAAGCCATGCTCCCACTTaaa                                          <  3:237850/149‑1 (MQ=255)                                                                                                                  tCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGTAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTAGCCGACTACCTTCCGGAAGCCATGCTCCCACTTAAATGCGa                                     <  6:134885/148‑1 (MQ=255)                                                                                                                  tCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGTAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTAGCCGACTACCTTCCGGAAGCCATGCTCCCACTTAAATGCGa                                     <  4:258612/148‑1 (MQ=255)                                                                                                                  tCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTAGCCGACTACCTTCCGGAAGCCATGCTCCCACTTAAATGCGa                                     <  1:267085/148‑1 (MQ=255)                                                                                                                   ccGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTAGCCGACTACCTTCCGGAAGCCATGCTCCCACTTAAATGCGAt                                    <  7:2389/148‑1 (MQ=255)                                                                                                                   ccGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGGACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTAGCCGACTACCTTCCGGAAGCCATGCTCCCTCTTAAATGCGAt                                    >  2:291138/1‑148 (MQ=255)                                                                                                                    cGTGAGTGCATTGAGGCGCTGTTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTAGCCGACTACCTTCCGGAAGCCATGCTCCCACTTAAATGCGATGa                                  <  7:353824/149‑1 (MQ=255)                                                                                                                    cGTGAGTGCATTGAGGCGCTGGTAGGTCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTAGCCGACTACCTTCCGGAAGCCATGCTCCCACTTAAATGCGATGa                                  <  1:353724/149‑1 (MQ=255)                                                                                                                    cGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTAGCCGACTACCTTCCGGAAGCCATGCTCCCACTTAAATGCGAt                                    >  4:311813/1‑147 (MQ=255)                                                                                                                         gTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTAGCCGACTACCTTCCGGAAGCCATGCTCCCACTTAAATGCGATGAGCTAg                             >  1:346810/1‑149 (MQ=255)                                                                                                                                       tGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTAGCCGACTACCTTCCGGAAGCCATGCTCCCACTTAAATGCGATGAGCTAGCGGCGCATGGGa                 >  5:144842/1‑147 (MQ=255)                                                                                                                                         gTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTAGCCGACTACCTTCCGGAAGCCATGCTCCCACTTAAATGCGATGAGCTAGCGGCGCATGGGAAAGt             >  4:54059/1‑149 (MQ=255)                                                                                                                                               cAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTAGCCGACTACCTTCCGGAAGCCATGCTCCCACTTAAATGCGATGAGCTAGCGGCGCATGGGAAAGTCCCCaa       <  2:346810/149‑1 (MQ=255)                                                                                                                                                aaGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTAGCCGACTACCTTCCGGAAGCCATGCTCCCACTTAAATGCGATGAGCTAGCGGCGCATGGGAAAGTCCCCAAg      <  3:311813/149‑1 (MQ=255)                                                                                                                                                   cAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTAGCCGACTACCTTCCGGAAGCCATGCTCCCACTTAAATGCGATGAGCTAGCGGCGCATGGGAAAGT‑CCCAAGGagc  <  5:42437/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      GAATTCTTGGGCGTTCCTTAATCGGGATGCTGACTACGCGTGCGTGCTGGCCATGGATCACATGATCCTGCAAAGCACAAATTATCTTCATTTTGATAACCACTTAGCGCTTTTCCGTGAGTGCATTGAGGCGCTGGTAGGGCAAGCAGGTGCATTGGATATCACCGCTATCGGGATGCGGTACGTCGATAAAATTGAACCAGCTGAAGGTGAGACCCTAGCCGACTACCTTCCGGAAGCCATGCTCCCACTTAAATGCGATGAGCTAGCGGCGCATGGGAAAGTCCCCAAGGAGC  >  NC_002947/23364‑23659

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