Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	623,513	G→A	8.5%	H56H (CAC→CAT)	PP_0536 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	623,513	0	G	A	8.5%	202.4 / 5.5	70	H56H (CAC→CAT)	PP_0536	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (29/35);  new base A (4/2);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

CCTCTTGCGTTGACTGCTCCAACATCTTGAGCTGTCCTGCAATTCGCGTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAATTTTGGTAGCCACCTCATGATGATAGGCGT  >  NC_002947/623365‑623655                                                                                                                                                      |                                                                                                                                               ccTCTTGCGTTGACTGCTCCAACATCTTGAGCTGTCCTGCAATTCGCGTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCAtggtgg                                                                                                                                                <  5:99115/149‑1 (MQ=255)   tctTGCGTTGACTGCTCCAACATCTTGAGCTGTCCTGCAATTCGCGTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGgtggtg                                                                                                                                              >  2:44575/1‑149 (MQ=255)          ttGACTGCTCCAACATCTTGAGCTGTCCTGCAATTCGCGTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGctct                                                                                                                                       <  4:253097/149‑1 (MQ=255)          ttGACTGCTCCAACATCTTGAGCTGTCCTGCAATTCGCGTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGctct                                                                                                                                       <  6:23200/149‑1 (MQ=255)          ttGACTGCTCCAACATCTTGAGCTGTCCTGCAATTCGCGTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGGCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGctct                                                                                                                                       <  7:118940/149‑1 (MQ=255)             aCTGCTCCAACATCTTGAGCTGTCCTGCAATTCGCGTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGc                                                                                                                                     <  1:66185/148‑1 (MQ=255)               tGCTCCAACATCTTGAGCTGTCCTGCAATTCGCGTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGTTGAATCGGTGCGTCGGTTGTCAGCAAGGCAAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGATCTGCAt                                                                                                                                   >  7:295992/1‑148 (MQ=255)                 cTCCAACATCTTGAGCTGTCCTGCAATTCGCGTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCAc                                                                                                                                >  2:1082/1‑149 (MQ=255)                    cAACATCTTGAGCTGTCCTGCAATTCGCGTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCAt                                                                                                                             >  6:111959/1‑149 (MQ=255)                      aCATCTTGAGCTGTCCTGCAATTCGCGTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTTACATTTTCATGGTGGTGCCGCTCTGCATCACCATg                                                                                                                            <  1:44575/149‑1 (MQ=255)                                cTGTCCTGCAATTCGCGTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCa                                                                                                                 <  2:102129/149‑1 (MQ=255)                                     cTGCAATTCGCGTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGc                                                                                                            >  7:11168/1‑149 (MQ=255)                                      tGCAATTCGCGTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTCGCt                                                                                                           >  5:194782/1‑149 (MQ=255)                                       gCAATTCGCGTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCt                                                                                                           >  7:98950/1‑148 (MQ=255)                                          aTTCGCGTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTaaa                                                                                                       <  1:29181/149‑1 (MQ=255)                                                gTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCAt                                                                                                                                   <  1:253450/115‑1 (MQ=255)                                                gTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCAt                                                                                                                                   >  2:253450/1‑115 (MQ=255)                                                    aGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCtt                                                                                              <  8:297537/148‑1 (MQ=255)                                                      gatttATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCa                                                                                           <  6:325265/147‑1 (MQ=255)                                                        tttATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAg                                                                                          >  1:302165/1‑148 (MQ=255)                                                        tttATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGt                                                                                         >  5:121368/1‑149 (MQ=255)                                                            tGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCAt                                                                                      <  3:85511/148‑1 (MQ=255)                                                               tCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGATGCAGAGCGGCAccacca                                                                                                                              >  5:295607/1‑105 (MQ=255)                                                               tCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGATGCAGAGCGGCAccacca                                                                                                                              <  6:295607/105‑1 (MQ=255)                                                                cTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCa                                                                                 >  3:25151/1‑149 (MQ=255)                                                                    cTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCg                                                                             <  6:121368/149‑1 (MQ=255)                                                                     tCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCg                                                                             >  2:308368/1‑148 (MQ=255)                                                                        aTGGCCATGACCCGTTGGGGAATCGGTGCGTCGGTTGTCACCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCaga                                                                         <  8:263731/149‑1 (MQ=255)                                                                             cATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCgt                                                                    <  1:37790/149‑1 (MQ=255)                                                                               tGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCgtgt                                                                  <  4:25151/149‑1 (MQ=255)                                                                                 aCCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTc                                                                 <  2:302165/148‑1 (MQ=255)                                                                                  cccGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATTGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCtt                                                               <  7:192295/149‑1 (MQ=255)                                                                                   ccGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTc                                                                 >  7:17271/1‑146 (MQ=255)                                                                                   ccGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTc                                                                 <  8:17271/146‑1 (MQ=255)                                                                                   ccGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTa                                                              <  6:129258/149‑1 (MQ=255)                                                                                   ccGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTa                                                              >  5:272167/1‑149 (MQ=255)                                                                                   ccGTTGGTGAATCGGTGCGTCGGTTGTCAACAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTa                                                              >  5:97981/1‑149 (MQ=255)                                                                                     gTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATc                                                            <  4:302040/149‑1 (MQ=255)                                                                                     gTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATc                                                            <  8:98950/149‑1 (MQ=255)                                                                                     gTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATc                                                            <  3:167475/149‑1 (MQ=255)                                                                                         gTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCctct                                                        >  1:283035/1‑149 (MQ=255)                                                                                         gTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCctct                                                        >  4:45044/1‑149 (MQ=255)                                                                                         gTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCctct                                                        >  1:226521/1‑149 (MQ=255)                                                                                            aaTCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCtctc                                                       >  4:118500/1‑147 (MQ=255)                                                                                                 gTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAAt                                                >  2:6147/1‑149 (MQ=255)                                                                                                    cGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCt                                             <  2:138890/149‑1 (MQ=255)                                                                                                      tCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCgtgt                                                                  >  1:79505/1‑126 (MQ=255)                                                                                                      tCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCgtgt                                                                  <  2:79505/126‑1 (MQ=255)                                                                                                      tCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTCGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCttt                                           <  6:194782/149‑1 (MQ=255)                                                                                                         gTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGAt                                        >  2:283431/1‑149 (MQ=255)                                                                                                             tCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGa                                    <  1:265813/149‑1 (MQ=255)                                                                                                             tCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGa                                    <  4:262020/149‑1 (MQ=255)                                                                                                              cAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAg                                   >  8:90559/1‑149 (MQ=255)                                                                                                                gCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCT‑CATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTaa                                >  4:136365/1‑149 (MQ=255)                                                                                                                     gCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAAtttt                            >  4:202290/1‑149 (MQ=255)                                                                                                                     gCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCAGAGTAAtttt                            <  3:118500/149‑1 (MQ=255)                                                                                                                         gCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAATTTTgg                          >  4:270452/1‑147 (MQ=255)                                                                                                                           cGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAATTTTGGTAGc                      >  1:177734/1‑149 (MQ=255)                                                                                                                              aTAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTTCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAATTTTGGTAGCCAc                   <  5:319319/149‑1 (MQ=255)                                                                                                                                aaGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGCCACCCTCAGAGCCGTGTCTTATACTATCCAGTAATTATTTGAACTGAGTAATTTTGGTAGCCAcc                  >  6:175461/1‑148 (MQ=255)                                                                                                                                  gCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAATTTTGGTAGCCACCTCa               <  7:90559/149‑1 (MQ=255)                                                                                                                                   ctctT‑CATTTTCATGGTTGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAATTTTGGTAGCCACCTCat              <  8:295992/149‑1 (MQ=255)                                                                                                                                    tctT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGACGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAATTTTGGTAGCAAcatcatc             >  6:191362/1‑148 (MQ=255)                                                                                                                                           cttgcATGGTGATGCAGAGCGGCACCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAAtttt                            <  2:199517/124‑1 (MQ=38)                                                                                                                                           cttgcATGGTGATGCAGAGCGGCACCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAAtttt                            >  1:199517/5‑128 (MQ=38)                                                                                                                                             ttCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAATTTTGGTAGCCACCTCATGATGATAg      <  6:272167/148‑1 (MQ=255)                                                                                                                                              tCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAATTTTGGTAGCCACCTCATGATGATAGGc    >  6:195299/1‑149 (MQ=255)                                                                                                                                               cATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAATTTTGGTAGCCACCTCATGATGATAGGc    <  3:202290/148‑1 (MQ=255)                                                                                                                                               cATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAATTTTGGTAGCCACCTCATGATGATAGGc    <  5:195299/148‑1 (MQ=255)                                                                                                                                                aTGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAATTTTGGTAGCCACCTCATGATGATAGGCGt  <  8:11168/149‑1 (MQ=255)                                                                                                                                                aTGGTGATGCAGAGCGGCACCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAATTTTGGTAGCCACCTCATGATGATAGGCGt  >  5:247781/1‑149 (MQ=255)                                                                                                                                                aTGGTGATGCAGAGCGGCACCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATAGGCACCCGCAGAGCCGTCTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAATTTTGGTAGCAACCTCATGATGATAGGCGt  >  2:177158/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                               CCTCTTGCGTTGACTGCTCCAACATCTTGAGCTGTCCTGCAATTCGCGTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAATTTTGGTAGCCACCTCATGATGATAGGCGT  >  NC_002947/623365‑623655

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