Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	623,522	T→G	8.6%	A53A (GCA→GCC)	PP_0536 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	623,522	0	T	G	8.6%	202.1 / 9.5	70	A53A (GCA→GCC)	PP_0536	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/34);  new base G (4/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

TTGACTGCTCCAACATCTTGAGCTGTCCTGCAATTCGCGTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAATTTTGGTAGCCACCTCATGATGATAGGCGTGCAGCTCTTT  >  NC_002947/623374‑623665                                                                                                                                                      |                                                                                                                                                ttGACTGCTCCAACATCTTGAGCTGTCCTGCAATTCGCGTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGctct                                                                                                                                                 <  6:23200/149‑1 (MQ=255) ttGACTGCTCCAACATCTTGAGCTGTCCTGCAATTCGCGTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGctct                                                                                                                                                 <  4:253097/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                                                                                                                                             >  2:253450/1‑115 (MQ=255)                                       gTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCAt                                                                                                                                             <  1:253450/115‑1 (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                                                                                                                                        <  6:295607/105‑1 (MQ=255)                                                      tCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGATGCAGAGCGGCAccacca                                                                                                                                        >  5:295607/1‑105 (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                                                                      <  8:98950/149‑1 (MQ=255)                                                                            gTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATc                                                                      <  4:302040/149‑1 (MQ=255)                                                                            gTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATc                                                                      <  3:167475/149‑1 (MQ=255)                                                                                gTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCctct                                                                  >  4:45044/1‑149 (MQ=255)                                                                                gTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCctct                                                                  >  1:283035/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                                      >  1:199517/5‑128 (MQ=38)                                                                                                                                  cttgcATGGTGATGCAGAGCGGCACCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAAtttt                                      <  2:199517/124‑1 (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)                                                                                                                                                ccGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAATTTTGGTAGCCACCTCATGATGATAGGCGTGCAGCTCtt   >  8:227949/1‑149 (MQ=255)                                                                                                                                                ccGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAATTTTGGTAGCCACCTCATGATGATAGGCGTGCAGCTCtt   <  1:224950/149‑1 (MQ=255)                                                                                                                                                 cGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAATTTTGGTAGCCACCTCATGATGATAGGCGTGCAGCTCttt  <  1:283431/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                TTGACTGCTCCAACATCTTGAGCTGTCCTGCAATTCGCGTGAAGGCTTTATGATCTTGCTCAATGGCCATGACCCGTTGGTGAATCGGTGCGTCGGTTGTCAGCAAGGCCAGCCGAATAAGCTCTT‑CATTTTCATGGTGGTGCCGCTCTGCATCACCATGCAAGGCTTCCAGTTGCTTAAACAACAGCTTGCAGTCATCGGCACCCGCAGAGCCGTGTCTTATCCTCTCCAGTAATTCTTTGATCTGAGTAATTTTGGTAGCCACCTCATGATGATAGGCGTGCAGCTCTTT  >  NC_002947/623374‑623665

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