Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	943,994	T→A	6.0%	N543I (AAC→ATC)	PP_0806 ←	surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	943,994	0	T	A	6.0%	239.7 / 4.9	83	N543I (AAC→ATC)	PP_0806	surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (37/41);  new base A (3/2);  total (40/43)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

CGGTCGGTGCGGCAGGGGCGGTGATGTCGGGCGCGGTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCGGTGGTCAGTGGTACATC  >  NC_002947/943848‑944139                                                                                                                                                    |                                                                                                                                                   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                                                                                      <  8:192554/149‑1 (MQ=39)                                                            aCGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGc                                                                                      >  1:152993/1‑149 (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                                                                              >  3:94135/1‑144 (MQ=39)                                                                         cGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCgg                                                                              <  4:94135/144‑1 (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                                                                                      >  1:98399/1‑123 (MQ=38)                                                                                      gTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGc                                                                                      <  2:98399/123‑1 (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                                                       <  2:36995/106‑1 (MQ=40)                                                                                                                                       gCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCa                                                       >  1:36995/1‑106 (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)                                                                                                                                                    |                                                                                                                                                   CGGTCGGTGCGGCAGGGGCGGTGATGTCGGGCGCGGTGATCTGTGCGGGCAGCGAGGAAACGCCGCCGTCATCGATGGCAACCACGTCCAGCAGT‑TCGCCGTTGGCTTGTGCGGGATTCAGTTCGATGCTGAAGCTGCCATCGGCGTTGGCAATGGCACTGCCAATCAGC‑GTGCCGTTGGCGTCATGCACCTCCACACGGCTGCCAGCCGGCGCGGTACCGGTGAGGGTGGTGCCATCAGCATCGATGGCCAGATTGCTCGGGCTGTCGGGCGCGGTGGTCAGTGGTACATC  >  NC_002947/943848‑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: ‑