Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,301,579	G→C	8.6%	intergenic (+69/‑1231)	PP_3774 → / → PP_5627	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,301,579	0	G	C	8.6%	178.1 / 6.4	58	intergenic (+69/‑1231)	PP_3774/PP_5627	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/27);  new base C (2/3);  total (28/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 = 7.68e-01

GGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTGATCTCTGCCTGCCT  >  NC_002947/4301432‑4301727                                                                                                                                                    |                                                                                                                                                        ggTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCgggg                                                                                                                                                        >  7:157273/1‑149 (MQ=255)     gACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGgggctgg                                                                                                                                                    >  7:81183/1‑149 (MQ=255)     gACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGgggctgg                                                                                                                                                    <  3:184978/149‑1 (MQ=255)           acCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCt                                                                                                                                              <  4:46638/149‑1 (MQ=255)            ccGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCCCCGCGCTGGGCttttt                                                                                                                                             <  6:271725/149‑4 (MQ=255)             cGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCCCCGCGCTGGGCtttttg                                                                                                                                            <  5:187890/149‑5 (MQ=255)              gACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATc                                                                                                                                           <  5:137110/149‑1 (MQ=255)               aCGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCCCCGCGCTGGGCtttttgca                                                                                                                                          <  4:135314/149‑7 (MQ=255)                     cccccGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTggg                                                                                                                                    >  2:111593/1‑149 (MQ=255)                     cccccGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTggg                                                                                                                                    >  3:182301/1‑149 (MQ=255)                      ccccGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCtttttgca                                                                                                                                          >  4:334310/1‑136 (MQ=255)                      ccccGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCtttttgca                                                                                                                                          <  3:334310/142‑7 (MQ=255)                      ccccGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGc                                                                                                                                   >  4:229509/1‑149 (MQ=255)                       cccGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCt                                                                                                                                  >  3:186785/1‑149 (MQ=255)                       cccGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCt                                                                                                                                  <  4:186785/149‑1 (MQ=255)                       cccGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCt                                                                                                                                  <  8:233674/149‑1 (MQ=255)                        ccgccgAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGgggctggg                                                                                                                                                   <  3:122338/131‑1 (MQ=255)                        ccgccgAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGgggctggg                                                                                                                                                   >  4:122338/1‑131 (MQ=255)                        ccgccgAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTg                                                                                                                                      <  7:235158/144‑1 (MQ=255)                        ccgccgAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTg                                                                                                                                      >  8:235158/1‑144 (MQ=255)                                gACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TAc                                                                                                                        >  2:254385/1‑149 (MQ=255)                                   ggCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTc                                                                                                                     >  2:202516/1‑149 (MQ=255)                                    gCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTcc                                                                                                                    <  2:302926/149‑1 (MQ=255)                                     cTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCg                                                                                                                   >  8:4539/1‑149 (MQ=255)                                        gATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGc                                                                                                                  <  6:187174/147‑1 (MQ=255)                                        gATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGc                                                                                                                  >  5:187174/1‑147 (MQ=255)                                           cTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGcc                                                                                                             >  7:188333/1‑149 (MQ=255)                                               atcatcAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTa                                                                                                         <  7:4539/149‑1 (MQ=255)                                                 catcaGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTaca                                                                                                       >  5:334137/1‑149 (MQ=255)                                                    caGCGAGCTCTCGCTCAACCATCGCTGAACGTCGATGCCCAGCACTGCTATCAAAACCACTCAGGGCGTCGCATGAATGCAAAAAGCCAAGCGCGGGGCTCGGCTTCTATCAGCTGGGCTCAACAG‑TACTTCCGCCAGCCTTTacaaaa                                                                                                    >  8:135568/1‑146 (MQ=255)                                                        gAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTa                                                                                                <  3:227979/149‑1 (MQ=255)                                                          gCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACt                                                                                              <  8:188333/149‑1 (MQ=255)                                                                gCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGt                                                                                      >  4:19812/1‑149 (MQ=255)                                                                     tccatcGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTc                                                                                 >  8:14515/1‑149 (MQ=255)                                                                     tccatcGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTc                                                                                 <  7:100344/149‑1 (MQ=255)                                                                     tccatcGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTc                                                                                 <  7:14515/149‑1 (MQ=255)                                                                     tccatcGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTc                                                                                 <  8:81183/149‑1 (MQ=255)                                                                      ccatcGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTc                                                                                 >  7:224128/1‑148 (MQ=255)                                                                       catcGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCa                                                                               <  3:229509/149‑1 (MQ=255)                                                                        atcGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAg                                                                              <  4:182301/149‑1 (MQ=255)                                                                         tcGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGt                                                                             <  1:202516/149‑1 (MQ=255)                                                                             tGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGgc                                                                         >  1:35568/1‑149 (MQ=255)                                                                             tGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGgc                                                                         >  2:26500/1‑149 (MQ=255)                                                                                         cccAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATaa                                                              <  8:157273/148‑1 (MQ=255)                                                                                          ccAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTc                                                                                 >  1:291798/1‑128 (MQ=255)                                                                                          ccAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTc                                                                                 <  2:291798/128‑1 (MQ=255)                                                                                             gCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATAACtata                                                         <  3:19812/149‑1 (MQ=255)                                                                                              cccTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATAACTATaa                                                        <  1:26500/149‑1 (MQ=255)                                                                                                 tGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTtgt                                                                   >  3:51471/1‑135 (MQ=255)                                                                                                 tGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTtgt                                                                   <  4:51471/135‑1 (MQ=255)                                                                                                         aaaCCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTCTGATAACTATAACGCTCGGCAtt                                             >  3:300924/1‑149 (MQ=255)                                                                                                          aaCCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTgg                                                                           <  3:356831/118‑1 (MQ=255)                                                                                                          aaCCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTgg                                                                           >  4:356831/1‑118 (MQ=255)                                                                                                               cTCAGGTCGTCGCATGAATGGAAAAAGCCCAGGGCGGGGCTTGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCGTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGGTTGTGATAACTATAACGCTAGGCATTGGATGc                                       <  4:300924/149‑1 (MQ=255)                                                                                                                     gCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGa                                 >  8:178601/1‑149 (MQ=255)                                                                                                                           cATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCCCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTc                           <  1:254385/149‑1 (MQ=255)                                                                                                                               aaTGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCcggc                       <  6:334137/149‑1 (MQ=255)                                                                                                                                  gatagaaGCCCAGCCCCGCGCTGGGCATCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCggg                    >  1:10945/6‑149 (MQ=255)                                                                                                                                   cAAAAAGCCCAGCGCGGGGCTGCGCTTCTATCAGCTGCGCTGAACAG‑TACAGACGCCAGCCTTTACACACTGTACT‑‑GGAAGGTACACCACTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGc                   >  1:89821/1‑149 (MQ=37)                                                                                                                                     agaaGCCCAGCCCCGCGCTGGGCTTTTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCt                 >  1:204199/3‑149 (MQ=255)                                                                                                                                      aaaaCCCCGAGCGGGGGTGGGGTTATATAACCTGGGGTGAAAAGATAATACCCCCCCCTCTTACCCACTTTACTAATGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGc                   <  1:111593/149‑1 (MQ=255)                                                                                                                                      aaaGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCt                 >  3:200508/1‑148 (MQ=255)                                                                                                                                          cccAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTGAtct            >  3:108970/1‑149 (MQ=255)                                                                                                                                                    ggctgggctTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTGATCTCTgcctgcct  <  2:35568/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                        GGTCGACTACACCGACGCAACCCCCGCCGAGGACGGCTCGATCTGCATCATCAGCGAGCTCTCGCTCATCCATCGCTGAACGTCGATGCCCAGCCCTGCTATCAAAACCCCTCAGGGCGTCGCATGAATGCAAAAAGCCCAGCGCGGGGCTGGGCTTCTATCAGCTGGGCTGAACAG‑TACTTCCGCCAGCCTTTACACACTGTACT‑‑GGAAGTTACTCCAGTGGGCGCTTGTGATAACTATAACGCTCGGCATTGGATGCCCGCGATGGGTCCGGCGGGCCTGATCTCTGCCTGCCT  >  NC_002947/4301432‑4301727

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