Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,191,853	A→G	11.8%	Q907R (CAG→CGG)	PP_1937 →	helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,191,853	0	A	G	11.8%	150.9 / 10.2	59	Q907R (CAG→CGG)	PP_1937	helicase
Reads supporting (aligned to +/- strand):  ref base A (26/26);  new base G (2/5);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

TGTATGAGCTTTATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAATGCTACTGCC  >  NC_002947/2191709‑2192000                                                                                                                                                 |                                                                                                                                                    tgtATGAGCTTTATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGcc                                                                                                                                                  <  7:66473/148‑1 (MQ=255)       aGCTTTATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGtcat                                                                                                                                            >  8:59840/1‑148 (MQ=255)        gCTTTATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCGGCTGGtcatca                                                                                                                                          <  5:157764/149‑1 (MQ=255)        gCTTTATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCGGCTGGtcatca                                                                                                                                          <  3:284462/149‑1 (MQ=255)            tATCGAAACGAAAGTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTTCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCTTGTCCGGCTCTGATGACCGGCTGGtcatcaga                                                                                                                                        <  8:136350/147‑3 (MQ=255)            tATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCGGCTGGtcatcaga                                                                                                                                        <  5:34407/147‑3 (MQ=255)                gaaacgaaacTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCcgc                                                                                                                                  <  2:289293/149‑1 (MQ=255)                 aaacgaaacTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGt                                                                                                                                               <  4:100035/135‑1 (MQ=255)                 aaacgaaacTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGt                                                                                                                                               >  3:100035/1‑135 (MQ=255)                  aacgaaacTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCcgc                                                                                                                                  >  1:135348/1‑147 (MQ=255)                    cgaaacTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCgaga                                                                                                                              <  6:239760/149‑1 (MQ=255)                     gaaacTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGAt                                                                                                                             >  1:101577/1‑149 (MQ=255)                      aaacTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGAtt                                                                                                                            <  8:81312/149‑1 (MQ=255)                      aaacTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGAtt                                                                                                                            >  8:247331/1‑149 (MQ=255)                           aaCAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAAc                                                                                                                       >  6:300453/1‑149 (MQ=255)                              aaGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCa                                                                                                                    <  6:3680/149‑1 (MQ=255)                                 gTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCt                                                                                                                 <  6:198269/149‑1 (MQ=255)                                     tGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAgg                                                                                                             <  2:135348/149‑1 (MQ=255)                                        cgcgTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCa                                                                                                          >  1:112043/1‑149 (MQ=255)                                               gAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTa                                                                                                    >  5:284631/1‑148 (MQ=255)                                                 gCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCa                                                                                                 <  2:112043/149‑1 (MQ=255)                                                 gCGTTGTCCCTTCCAGTGAACCTTGGGTGTTT‑ACTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGAtt                                                                                                                            <  6:288691/121‑1 (MQ=255)                                                 gCGTTGTCCCTTCCAGTGAACCTTGGGTGTTT‑ACTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGAtt                                                                                                                            >  5:288691/1‑121 (MQ=255)                                                  cGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAg                                                                                                <  5:300453/149‑1 (MQ=255)                                                     tGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGAc                                                                                              >  3:81826/1‑148 (MQ=255)                                                         ccTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCATGGTTATCAGACGAttt                                                                                         >  3:82315/1‑149 (MQ=255)                                                         ccTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGAttt                                                                                         >  1:168956/1‑149 (MQ=255)                                                           ttCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTg                                                                                        >  2:304304/1‑148 (MQ=255)                                                               aGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGaa                                                                                   <  2:101577/149‑1 (MQ=255)                                                               aGTGAACCTTGGGTGTTTCTCTGCTTGGCCAGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGaa                                                                                   >  1:183527/1‑149 (MQ=255)                                                                     ccTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCaa                                                                             >  8:49871/1‑149 (MQ=255)                                                                              tttCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGGTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGAAGAACCTCAACCCGACaa                                                                     <  4:81826/148‑1 (MQ=255)                                                                                 ctctGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAggg                                                                  <  7:247331/148‑1 (MQ=255)                                                                                  tctGCTTGGACGGTGAGAATGCGATCGCACCAAGGTGGATTTTCATGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACCCAGCTGAGGGCACGGTTATCAGACGATTTGCATAACCTCAACCCGACAAGGGCt                                                                >  4:267293/1‑149 (MQ=255)                                                                                      cTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCATGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCaa                                                            <  4:82315/149‑1 (MQ=255)                                                                                         ggCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATgg                                                         >  3:230393/1‑149 (MQ=255)                                                                                          gCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATggg                                                        <  7:49871/149‑1 (MQ=255)                                                                                          gCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATgg                                                         <  1:32540/148‑1 (MQ=255)                                                                                             ggTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCAt                                                     >  8:275093/1‑149 (MQ=255)                                                                                                 agaATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCAtgctg                                                 <  2:183527/149‑1 (MQ=255)                                                                                                   aaTGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCAtgctgc                                                >  6:127468/1‑148 (MQ=255)                                                                                                      gcgcgcgcACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATc                                            >  2:69186/1‑149 (MQ=255)                                                                                                               ccAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTg                                     >  6:33905/1‑147 (MQ=255)                                                                                                               ccAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTg                                     <  5:33905/147‑1 (MQ=255)                                                                                                                    gTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGc                              >  5:253266/1‑149 (MQ=255)                                                                                                                       gATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTgg                           <  4:230393/149‑1 (MQ=255)                                                                                                                             ccTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCgag                      <  2:168956/148‑1 (MQ=255)                                                                                                                               tGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATg                   >  3:1852/1‑149 (MQ=255)                                                                                                                                 tCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATgg                  <  7:275093/148‑1 (MQ=255)                                                                                                                                 tCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATg                   <  1:304304/147‑1 (MQ=255)                                                                                                                                   cGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGgag                <  6:54044/148‑1 (MQ=255)                                                                                                                                       attgatgaCCGGCTGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGaaa             <  8:190226/145‑1 (MQ=255)                                                                                                                                       attgatgaCCGGCTGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAAt           >  3:105197/3‑149 (MQ=255)                                                                                                                                       attgatgaCCGGCTGGTCATCAATATCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGaaaaa           >  7:190226/3‑148 (MQ=255)                                                                                                                                           atgaCCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCt                             >  2:161926/1‑127 (MQ=255)                                                                                                                                           atgaCCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCt                             <  1:161926/127‑1 (MQ=255)                                                                                                                                           atgaCCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAATGCt        <  6:284631/148‑1 (MQ=255)                                                                                                                                           atgaCCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACCGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGCGCTGCAATGGCCATGCTGCAATCGTTTCTGGCCAGGATGGCCGAGATGGAGAAAATGCTa       >  6:158544/1‑149 (MQ=255)                                                                                                                                                cAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAATGCTACTGcc  >  4:180910/1‑149 (MQ=255)                                                                                                                                                 aGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGAGCGAGATGGAGAAAATGCTACTGcc  >  6:187874/1‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                                    TGTATGAGCTTTATCGAAACGAAACTAACAAGGTCATGTCGCGTGTGAGCGTTGTCCCTTCCAGTGAACCTTGGGTGTTTCTCTGCTTGGCCGGTGAGAATGCGCGCGCACCAAGGTGGATTTTCCTGTCCGGCTCTGATGACCAGCCGGTCATCAATCTCCGCGAGATTGCAACGCAGCTGAGGGCACGGTTATCAGACGATTTGCAGAACCTCAACCCGACAAGGGCTGCAATGGGCATGCTGCAATCGTTTCTGGCCAGGCTGGGCGAGATGGAGAAAATGCTACTGCC  >  NC_002947/2191709‑2192000

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