Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	943,991	G→C	5.6%	A544G (GCC→GGC)	PP_0806 ←	surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	943,991	0	G	C	5.6%	249.7 / 4.0	84	A544G (GCC→GGC)	PP_0806	surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (38/41);  new base C (3/2);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTCGGTCGGTGCGGCAGGGGCGGTGATGTCGGGCGCGGTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCGGTGGTCAGTGGTACATC  >  NC_002947/943844‑944139                                                                                                                                                     |                                                                                                                                                      aGTTCGGTCGGTGCGGCAGGGGCGGTGATGACGGGCGCGGTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGc                                                                                                                                                      >  1:36310/1‑149 (MQ=17)  gTTCGGTCGGTGCGGCAGGGGCGGTGATGTCGGGCGCGGTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCg                                                                                                                                                     >  1:357592/1‑149 (MQ=32)     cggtcggtGCGGCAGGGGCGGTGATGTCGGGCGCGGTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGtt                                                                                                                                                   <  6:65291/148‑1 (MQ=32)     cggtcggtGCGGCAGGGGCGGTGATGTCGGGCGCGGTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTg                                                                                                                                                  >  5:78306/1‑149 (MQ=32)       gtcggtGCGGCAGGGGCGGTGATGTCGGGCGCGGTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGc                                                                                                                                                <  4:118228/149‑1 (MQ=32)          ggtGCGGCAGGGGCGGTGATGTCGGGCGCGGTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAAt                                                                                                                                             >  2:198527/1‑149 (MQ=32)                  aGGGGCGGTGATGTCGGGCGCGGTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGc                                                                                                                                     <  6:217090/149‑1 (MQ=32)                    gggCGGTGATGTCGGGCGCGGTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGcc                                                                                                                                    >  3:350332/1‑148 (MQ=32)                    gggCGGTGATGTCGGGCGCGGTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCa                                                                                                                                   <  6:78306/149‑1 (MQ=32)                        ggTGATGTCGGGCGCGGTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCa                                                                                                                               <  2:106452/149‑1 (MQ=32)                          tGATGTCGGGCGCGGTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGc                                                                                                                             >  1:304161/1‑149 (MQ=33)                           gATGTCGGGCGCGGTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGc                                                                                                                             <  7:1228/148‑1 (MQ=33)                            aTGTCGGGCGCGGTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCGTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑Gt                                                                                                                          >  5:356396/1‑149 (MQ=21)                             tggcGGGCGCGGTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCCCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTg                                                                                                                         <  6:244558/146‑1 (MQ=21)                               tCGGGCGCGGTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAATGTCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGATGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGc                                                                                                                        >  1:120911/1‑149 (MQ=255)                                 gggCGCGGTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGt                                                                                                                     >  7:192554/1‑149 (MQ=34)                                  ggCGCGGTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGtt                                                                                                                    <  7:251494/149‑1 (MQ=34)                                     gcgGTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGc                                                                                                                 <  6:28941/149‑1 (MQ=34)                                       ggTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGt                                                                                                               >  5:313419/1‑149 (MQ=35)                                        gTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGt                                                                                                               >  5:190703/1‑148 (MQ=35)                                           aTCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCAt                                                                                                            <  6:190703/148‑1 (MQ=37)                                           aTCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATg                                                                                                           >  4:330852/1‑149 (MQ=37)                                             cTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCa                                                                                                         >  2:283682/1‑149 (MQ=37)                                               gtgCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCAcc                                                                                                       >  1:362185/1‑149 (MQ=38)                                               gtgCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCAcc                                                                                                       >  8:325085/1‑149 (MQ=38)                                                tgCGGGCAGCGAGGAAACGCCGCCGCCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGc                                                                                                          <  5:342780/145‑1 (MQ=37)                                                tgCGGGCACCGAGGAAACGCAGCCGTCATCGATGGCAACCACGTCCAGCAGG‑TCGCCGTTGGCTTGTGCGGGAATCACATCGATGCTGAAGCTGCCATCGGCGTTGTCAAGGGCACTGCCAATCAGC‑GCGCCGTTGGCGTCATGCACCt                                                                                                      <  6:88556/149‑1 (MQ=11)                                                         cGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGc                                                                                             >  8:140521/1‑149 (MQ=39)                                                          gAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAAt                                                                                                                                 <  8:279340/113‑1 (MQ=32)                                                              aaaCGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCa                                                                                        <  2:226912/149‑1 (MQ=39)                                                                aCGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGc                                                                                      >  1:152993/1‑149 (MQ=39)                                                                aCGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGc                                                                                      <  8:192554/149‑1 (MQ=39)                                                                   ccgccgTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCgg                                                                                   <  7:153476/149‑1 (MQ=39)                                                                      ccgTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGgcgc                                                                                <  4:350332/149‑1 (MQ=39)                                                                       cgTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATAGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGcgcg                                                                               >  7:51388/1‑149 (MQ=255)                                                                            tCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTAAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTAcc                                                                          >  7:362428/1‑149 (MQ=255)                                                                             cGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCgg                                                                              <  4:94135/144‑1 (MQ=39)                                                                             cGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCgg                                                                              >  3:94135/1‑144 (MQ=39)                                                                                tGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTg                                                                      <  2:36310/149‑1 (MQ=39)                                                                                    aaCCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAggg                                                                  >  6:18650/1‑149 (MQ=38)                                                                                     aCCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGggt                                                                 <  4:183227/149‑1 (MQ=37)                                                                                          gTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGc                                                                                      <  2:98399/123‑1 (MQ=38)                                                                                          gTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGc                                                                                      >  1:98399/1‑123 (MQ=38)                                                                                              agcagT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTcc                                                                                                    >  3:42698/1‑105 (MQ=40)                                                                                              agcagT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTcc                                                                                                    <  4:42698/105‑1 (MQ=40)                                                                                              agcagT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATc                                                        <  6:327801/149‑1 (MQ=37)                                                                                                  gT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGc                                                     >  7:319342/1‑148 (MQ=37)                                                                                                   t‑tCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCa                                                    <  2:357592/148‑1 (MQ=37)                                                                                                        ccGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGAt                                               <  1:198527/149‑1 (MQ=37)                                                                                                        ccGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGAt                                               <  2:304161/149‑1 (MQ=37)                                                                                                          gTTGGCTTGTGCGGGATGCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATgg                                             <  3:330852/149‑1 (MQ=37)                                                                                                             ggCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGcc                                           <  2:120911/148‑1 (MQ=37)                                                                                                                   tgCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGc                                    <  2:320657/149‑1 (MQ=37)                                                                                                                      gggATTCAGTTCGATGCTAAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCg                                 <  8:362428/149‑1 (MQ=37)                                                                                                                        gATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCggg                               >  4:283787/1‑149 (MQ=37)                                                                                                                        gATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGCCTGACAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCggg                               >  6:135775/1‑149 (MQ=255)                                                                                                                         aTTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGc                              >  2:205815/1‑149 (MQ=37)                                                                                                                             aGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTc                          >  5:200808/1‑149 (MQ=37)                                                                                                                               ttCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCgg                        >  3:197819/1‑149 (MQ=37)                                                                                                                                tCGTTGCTGAAGCTGCCATGGGCGTTGGCAATGGCACTGCCAATGAGCGGTGCGGTTGG‑GTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGGGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGGCgg                        <  7:27788/148‑1 (MQ=11)                                                                                                                                 cGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGgc                      >  7:24463/1‑149 (MQ=37)                                                                                                                                  gATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGgcg                     >  1:71354/1‑149 (MQ=37)                                                                                                                                       tGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCggtg                >  5:266912/1‑149 (MQ=37)                                                                                                                                         ccattGCCAACGCCGATGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCggtggt              >  1:347701/5‑149 (MQ=255)                                                                                                                                         ccattGCCAACGCCGATGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACATCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCggtggt              >  8:88162/5‑149 (MQ=11)                                                                                                                                         aaGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCCCACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCggtggt              <  6:313419/149‑1 (MQ=37)                                                                                                                                          cattGCCAACGCCGATGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGGGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCGGTGGTc             >  8:79900/4‑149 (MQ=11)                                                                                                                                           gCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCa                                                       >  1:36995/1‑106 (MQ=40)                                                                                                                                           gCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCa                                                       <  2:36995/106‑1 (MQ=40)                                                                                                                                           gCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCgg                        >  7:299020/1‑137 (MQ=37)                                                                                                                                           gCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCgg                        <  8:299020/137‑1 (MQ=37)                                                                                                                                           gCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCggtggt              >  2:310862/1‑147 (MQ=37)                                                                                                                                             tGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCGGTGGTCAGt          >  3:127077/1‑149 (MQ=37)                                                                                                                                              gCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATTCACCTCCACACTGCTACCAGCCGGCGCGGTACCGTTGATGGTTGTGCCATCAGCATCGATTGCCTGATTGGTCGTGATGTCGGGCGCGGTGGTCAGTg         <  6:356396/149‑1 (MQ=11)                                                                                                                                              gCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCGGTGGTCAGTg         >  4:170972/1‑149 (MQ=37)                                                                                                                                              gCCAACGCCGATGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCGGTGGTCAGt          <  7:88162/148‑1 (MQ=255)                                                                                                                                              gCCAACGCCGATGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCGGTGGTCAGt          <  2:347701/148‑1 (MQ=255)                                                                                                                                               ccctCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGGCATGCACCGCCCCACGGCTGCCAGCCTGCGCGCTACCGGGGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCGGTGGTCAGTg         <  2:301801/145‑1 (MQ=11)                                                                                                                                               ccATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCGGTGGTCAGTg         <  5:18650/148‑1 (MQ=37)                                                                                                                                                cATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATTCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCGGTGGTCAGTGGt       <  4:293812/149‑1 (MQ=37)                                                                                                                                                cATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCGGTGGTCAGTGGt       <  8:24463/149‑1 (MQ=37)                                                                                                                                                cATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCGGTGGTCAGTGGt       <  8:319342/149‑1 (MQ=37)                                                                                                                                                 aTCGGCGTTGGCAATGGCACTGCCAATCAGA‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCGGTGGTCAGTGGTa      <  7:325085/149‑1 (MQ=37)                                                                                                                                                  tCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCGGTGGTCAGTGGTa      >  3:75169/1‑148 (MQ=37)                                                                                                                                                    ggCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCGGTGGTCAGTGGTACa    <  3:170972/148‑1 (MQ=37)                                                                                                                                                     gCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCGGTGGTCAGTGGTACATc  >  5:277352/1‑149 (MQ=37)                                                                                                                                                     |                                                                                                                                                      AGTTCGGTCGGTGCGGCAGGGGCGGTGATGTCGGGCGCGGTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCGGTGGTCAGTGGTACATC  >  NC_002947/943844‑944139

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