Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,743,847	A→C	7.1%	L673R (CTC→CGC)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,743,847	0	A	C	7.1%	176.1 / 4.3	56	L673R (CTC→CGC)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base A (24/28);  new base C (2/2);  total (26/30)
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.34e-01

GAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAGGAACACCACCTTCAGCAGGC  >  NC_002947/5743700‑5743994                                                                                                                                                             |                                                                                                                                                    gAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGa                                                                                                                                                     <  8:28723/148‑1 (MQ=255) gAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAg                                                                                                                                                    >  3:137381/1‑149 (MQ=255)                        gcCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCAtt                                                                                                                                           <  5:90220/135‑1 (MQ=255)                        gcCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCAtt                                                                                                                                           >  6:90220/1‑135 (MQ=255)                        gcCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCa                                                                                                                             >  2:37456/1‑149 (MQ=255)                         ccAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCAt                                                                                                                            <  1:30794/149‑1 (MQ=255)                         ccAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCAt                                                                                                                            <  1:37456/149‑1 (MQ=255)                                 aCCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTg                                                                                                                     <  5:96312/148‑1 (MQ=255)                                   cGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTa                                                                                                                  >  1:18317/1‑149 (MQ=255)                                     gCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATCGTGAGCGCCCCATTGAGGCCGAACTTCATGTTGATGGTAcc                                                                                                                >  7:114533/1‑149 (MQ=255)                                              cgcgcTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTc                                                                                                       <  4:35865/149‑1 (MQ=255)                                              cgcgcTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTc                                                                                                       <  8:62535/149‑1 (MQ=255)                                               gcgcTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCg                                                                                                      <  7:233227/149‑1 (MQ=255)                                               gcgcTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCg                                                                                                      <  5:202367/149‑1 (MQ=255)                                                  cTTCAACATGCTGCGAGCTCAGGCCAAATATGAACATGTTGTCGGCCCCCACAAGCTCTGACTCCCTATATCATTCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGaa                                                                                                            <  1:260375/149‑1 (MQ=255)                                                   ttCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAAc                                                                                                  >  6:245483/1‑149 (MQ=255)                                                   ttCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAAc                                                                                                  >  8:70058/1‑149 (MQ=255)                                                     cAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAAcc                                                                                                 <  5:245483/148‑1 (MQ=255)                                                      aaCCTGCTGCGAGGTCAGACCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCcg                                                                                               <  2:218725/149‑1 (MQ=255)                                                            ctgcGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGt                                                                                           <  7:70058/147‑1 (MQ=255)                                                            ctgcGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGt                                                                                           <  1:186617/147‑1 (MQ=255)                                                            ctgcGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTg                                                                                          >  4:78268/1‑148 (MQ=255)                                                             tgcGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTc                                                                                                                              <  7:26523/111‑1 (MQ=255)                                                             tgcGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTc                                                                                                                              >  8:26523/1‑111 (MQ=255)                                                                         ccAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCgg                                                                             >  8:139750/1‑148 (MQ=255)                                                                                 gAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCATGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTccg                                                                    >  6:236376/1‑149 (MQ=255)                                                                                     aTGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGc                                                                                                                       <  3:89882/94‑1 (MQ=255)                                                                                     aTGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGc                                                                                                                       >  4:89882/1‑94 (MQ=255)                                                                                       gttgtCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTCCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTccgccgccg                                                              >  2:210500/1‑149 (MQ=255)                                                                                          gtCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGaa                                                                           >  5:76735/1‑133 (MQ=255)                                                                                          gtCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGaa                                                                           <  6:76735/133‑1 (MQ=255)                                                                                           tCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGgat                                                          <  7:139750/149‑1 (MQ=255)                                                                                             ggCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCATGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGgatg                                                         <  5:236376/148‑1 (MQ=255)                                                                                              gCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGgatgat                                                       >  3:244279/1‑149 (MQ=255)                                                                                               cccccAC‑‑‑‑‑CTG‑TTCGG‑ACATC‑‑TCGACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCTCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGatgatg                                                      <  2:101841/149‑1 (MQ=255)                                                                                               cccccAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGatgatg                                                      <  2:18317/149‑1 (MQ=255)                                                                                                ccccAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGc                                                     >  3:307263/1‑149 (MQ=255)                                                                                                           cTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCgg                                               <  4:307263/149‑1 (MQ=255)                                                                                                                  gC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGc                                         >  1:3134/1‑149 (MQ=255)                                                                                                                       aTC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGAc                                     >  6:80233/1‑149 (MQ=255)                                                                                                                       aTC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGAc                                     >  3:28098/1‑149 (MQ=255)                                                                                                                            tcCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGAcatg                                  >  4:186710/1‑146 (MQ=255)                                                                                                                               aCGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACAATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGtgg                                 >  2:256470/1‑145 (MQ=255)                                                                                                                                     gcgcCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTgg                          <  3:83944/148‑1 (MQ=255)                                                                                                                                      cgcCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGc                        <  3:78268/149‑1 (MQ=255)                                                                                                                                      cgcCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGc                        <  4:137381/149‑1 (MQ=255)                                                                                                                                         cATCCAGGGTGCCAATGGTGAGCGCACCATTGAGTCCGAAGTTCATGTTGCTGGTACCGGAAGCGTGGTAACCCGCGGGGTATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGTTGACGTTGTAGTTGGGCAgg                     <  8:114533/149‑1 (MQ=255)                                                                                                                                             cAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTCCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAGGAaca                 <  1:210500/149‑1 (MQ=255)                                                                                                                                                    tcAATGGTGCGCTCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAgg                                                                        <  6:296599/86‑1 (MQ=255)                                                                                                                                                    tcAATGGTGCGCTCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAgg                                                                        >  5:296599/2‑87 (MQ=255)                                                                                                                                                    tcAATGGTGCGCTCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCa                       >  2:304103/2‑136 (MQ=255)                                                                                                                                                    tcAATGGTGCGCTCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCa                       <  1:304103/135‑1 (MQ=255)                                                                                                                                                      aaTGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGa                                      <  7:79905/119‑1 (MQ=255)                                                                                                                                                      aaTGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGa                                      >  8:79905/1‑119 (MQ=255)                                                                                                                                                      aaTGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAGGAACACCACCtt          >  2:11963/1‑147 (MQ=255)                                                                                                                                                      aaTGGTGAGCGCACCATTGAGGACGAACTTCAAGTTGGTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAGGAACACCACCtt          >  1:23737/1‑147 (MQ=255)                                                                                                                                                       aTGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAGGAACACCACCTTcag       <  4:28098/149‑1 (MQ=255)                                                                                                                                                             aGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAGGAACACCACCTTCAGCAGGc  <  7:43950/148‑1 (MQ=255)                                                                                                                                                             |                                                                                                                                                    GAAGCGGCTATGGCAGCATTGGCGCCAAAATCACCGGCACGCTTGCGCGCTTCAACCTGCTGCGAGGTCAGGCCAAAGATGAACATGTTGTCGGCCCCCAC‑‑‑‑‑CTG‑TTCGC‑ACATC‑‑TCCACGTTGGCGCCATCCAGGGTGCCAATGGTGAGCGCACCATTGAGGCCGAACTTCATGTTGCTGGTACCGGAAGCTTCGTAACCCGCGGTGGATATCTGCTCGGAAAGGTCCGCCGCCGGGATGATGCTCTCGGCCAGGCTGACGTTGTAGTTGGGCAGGAACACCACCTTCAGCAGGC  >  NC_002947/5743700‑5743994

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