Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,762,926	T→C	7.6%	F57L (TTC→CTC)	PP_3326 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,762,926	0	T	C	7.6%	192.2 / 3.6	67	F57L (TTC→CTC)	PP_3326	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/29);  new base C (3/2);  total (35/32)
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.98e-01

GCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGT  >  NC_002947/3762779‑3763068                                                                                                                                                    |                                                                                                                                                  gCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGtt                                                                                                                                                  >  5:218476/1‑149 (MQ=255)       cTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGt                                                                                                                                           >  6:313488/1‑149 (MQ=255)         gTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGa                                                                                                                                         >  5:299835/1‑149 (MQ=255)          ttGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGa                                                                                                                                         <  6:299835/148‑1 (MQ=255)          ttGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGAc                                                                                                                                        >  6:157355/1‑149 (MQ=255)             gCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCgg                                                                                                                                     <  5:313488/149‑1 (MQ=255)             gCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCgg                                                                                                                                     <  4:76445/149‑1 (MQ=255)              cccTGATGGTCAGCGTCCACGCAGGAGCCCACGGCCTGTGGACCGATGATCGACGAGCAAATATCGAGGCCGTGTACGGACACGGAGCCGATGATAACGCCGACAAGTCGGAAAAGATCAAGCGCGCCTGCGCGTTCGA‑CCTGACCgg                                                                                                                                     >  1:194606/1‑148 (MQ=255)               ccTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCggg                                                                                                                                    >  8:308463/1‑148 (MQ=255)               ccTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCggg                                                                                                                                    >  5:26367/1‑148 (MQ=255)               ccTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCgg                                                                                                                                     >  3:227686/1‑148 (MQ=255)                   aTGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAAc                                                                                                                                 <  6:26367/147‑1 (MQ=255)                   aTGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACa                                                                                                                                <  6:303119/148‑1 (MQ=255)                      gTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGAt                                                                                                                            <  1:166914/149‑1 (MQ=255)                      gTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGAt                                                                                                                            >  6:36057/1‑149 (MQ=255)                         aGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATc‑‑cc                                                                                                                       <  4:109928/149‑1 (MQ=255)                           cGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCgg                                                                                                                     >  3:61459/1‑149 (MQ=255)                           cGTAAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGAAc‑‑cccg                                                                                                                     >  1:344304/1‑147 (MQ=255)                            gTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACCGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATGTTC‑‑CCgg                                                                                                                     <  4:227686/149‑1 (MQ=255)                               aaCGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGCGACCGGGAACATGATC‑‑CCGGTCAc                                                                                                                 >  7:305682/1‑149 (MQ=255)                                    aGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTCTACGGCCACGGCGCCGAGGATAACGACTACAAGGCGGAAAAGATCAAGCGCTCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTc                                                                                                             >  7:256404/1‑148 (MQ=255)                                                 cTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGa                                                                                               >  1:220716/1‑149 (MQ=255)                                                  tgtgGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATGTTC‑‑CCgg                                                                                                                     <  2:326333/127‑1 (MQ=39)                                                  tgtgGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATGTTC‑‑CCgg                                                                                                                     >  1:326333/1‑127 (MQ=39)                                                  tgtgGACAGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATGTTC‑‑CCgct                                                                                                                    >  5:199724/1‑126 (MQ=38)                                                      gACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATg                                                                                          <  2:220716/149‑1 (MQ=255)                                                          gaggagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGACGGATCATGCCGt                                                                                      >  1:191266/1‑149 (MQ=255)                                                           aggagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGtg                                                                                     <  5:157355/149‑1 (MQ=255)                                                             gagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGCACATGATC‑‑CCGGTCACGGGCGAGCGCCTGCCGGATCATGCCGagct                                                                                   >  2:177124/1‑149 (MQ=255)                                                              agCGACGTGGAAATATCGAGGCCGGGTACGGCCTCGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGtgctg                                                                                  <  2:344304/149‑1 (MQ=255)                                                                cGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGtgctgc                                                                                 >  1:44794/1‑148 (MQ=255)                                                                  aCGTGGCAATATCGAGGCCGTGTACGGCCAAGGCGCCGAGGATAACGCCTACAAGGAGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGAACGGGAAc                                                                                                                                 >  7:265049/1‑100 (MQ=255)                                                                    gTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCa                                                                            <  5:36057/149‑1 (MQ=255)                                                                     tGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCAc                                                                           >  7:213594/1‑149 (MQ=255)                                                                          atatCGAGGCCGTGTACGGCCTCGGCGCCGAGGTTAACGGCTAGAAGGCGGAAAAGATCAAGCGCGCCTGGGAGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGACCATGCCGTGCTGCAACGACTTAcc                                                                      <  5:234523/149‑1 (MQ=255)                                                                              cGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGc                                                                  >  3:144085/1‑149 (MQ=255)                                                                               gAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGc                                                                  <  4:277695/148‑1 (MQ=255)                                                                                           ggCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGc                                                     >  3:336397/1‑149 (MQ=255)                                                                                           ggCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGc                                                     <  8:124526/149‑1 (MQ=255)                                                                                             ccACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCg                                                   >  5:387725/1‑149 (MQ=255)                                                                                             ccACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCg                                                   >  2:23909/1‑149 (MQ=255)                                                                                              cACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGt                                                  <  4:88985/149‑1 (MQ=255)                                                                                               aCGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtg                                                 <  2:194986/149‑1 (MQ=255)                                                                                                     ccGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGctgctg                                           <  4:144085/149‑1 (MQ=255)                                                                                                       gAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTg                                                      >  6:219600/1‑136 (MQ=255)                                                                                                       gAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTg                                                      <  5:219600/136‑1 (MQ=255)                                                                                                            tAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCg                                                   <  7:167677/134‑1 (MQ=255)                                                                                                            tAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCg                                                   >  8:167677/1‑134 (MQ=255)                                                                                                             aaCGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCAGGCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGCCAAc                                     >  1:113378/1‑149 (MQ=255)                                                                                                              aCGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACgg                                   <  2:44794/148‑1 (MQ=255)                                                                                                                   tACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTg                             >  1:295581/1‑149 (MQ=255)                                                                                                                     cAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTg                             <  7:308463/147‑1 (MQ=255)                                                                                                                               aaGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGccc                 <  4:61459/149‑1 (MQ=255)                                                                                                                               aaGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGccc                 >  2:96966/1‑149 (MQ=255)                                                                                                                               aaGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGACGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGccc                 <  2:191266/149‑1 (MQ=255)                                                                                                                                  aTCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGa              >  1:229539/1‑149 (MQ=255)                                                                                                                                     aaGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACaa           <  6:218476/149‑1 (MQ=255)                                                                                                                                     aaGCGCGCCTGGGCGTTCGA‑CGCGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACaa           <  8:305682/149‑1 (MQ=255)                                                                                                                                         gcgcCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACg                                    >  4:218623/1‑120 (MQ=255)                                                                                                                                         gcgcCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACg                                    <  3:337707/120‑1 (MQ=255)                                                                                                                                         gcgcCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACg                                    <  3:218623/120‑1 (MQ=255)                                                                                                                                         gcgcCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACg                                    >  4:337707/1‑120 (MQ=255)                                                                                                                                            ccTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCt                                                                <  4:24949/89‑1 (MQ=255)                                                                                                                                            ccTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCt                                                                >  3:24949/1‑89 (MQ=255)                                                                                                                                            ccTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg    <  8:213594/149‑1 (MQ=255)                                                                                                                                             cTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATggg   >  3:267441/1‑149 (MQ=255)                                                                                                                                             cTGGGCGGTAGA‑CGGTACCGGGAACATGATC‑‑CCGGGCACGGGCGAGCGCCTGGCGGAGCATGCCGTGCTGCAACGACTTACCCCGCCTGCGGCAATTGCCGGGCGGCTGGACAACGGTCCGTGGTAACAGCGCCCGGACAATCAATggg   <  2:194606/149‑1 (MQ=255)                                                                                                                                              tGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGt  >  4:323501/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GCTTCCCTGTTGGCCCTGATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGT  >  NC_002947/3762779‑3763068

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑