Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,613,523	T→A	8.6%	A158A (GCT→GCA)	PP_4933 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,613,523	0	T	A	8.6%	194.3 / 7.7	70	A158A (GCT→GCA)	PP_4933	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/33);  new base A (2/4);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

CACACGGCCGTACAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAACTCAGCT  >  NC_002947/5613380‑5613671                                                                                                                                                 |                                                                                                                                                      cacaCGGCCGTACAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCg                                                                                                                                                    <  8:278724/147‑1 (MQ=255) cacaCGGCCGTACAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGa                                                                                                                                                  <  2:70814/149‑1 (MQ=255)       gCCGTACAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGCAGAAACACCCCACTTCGTGCTGGCTCCGGAGTCCg                                                                                                                                             <  1:54972/149‑1 (MQ=255)           tACAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGcc                                                                                                                                        <  6:26928/149‑1 (MQ=255)             cAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCACCGGACTCCGGAGccag                                                                                                                                      <  5:255938/149‑3 (MQ=255)             cAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCACCGGACTCCGGAGcca                                                                                                                                       <  4:282630/148‑2 (MQ=255)              aGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCa                                                                                                                                    >  1:94057/1‑149 (MQ=255)                gAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCACCGGACTCCGGAGccag                                                                                                                                      <  4:123768/146‑3 (MQ=255)                gAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCACCGGACTCCGGAGccag                                                                                                                                      >  3:123768/1‑144 (MQ=255)                  aGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCATTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGcc                                                                                                                                <  6:419561/149‑1 (MQ=255)                  aGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGC‑GGCACCGGACTCCGGAGCCAGCACgaa                                                                                                                                <  2:294184/149‑3 (MQ=255)                    tGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCg                                                                                                                              <  1:143312/149‑1 (MQ=255)                    tGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCg                                                                                                                              >  6:6204/1‑149 (MQ=255)                      gCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGt                                                                                                                            >  5:365106/1‑149 (MQ=255)                            gagTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTg                                                                                                                      <  4:61056/149‑1 (MQ=255)                               tGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCgg                                                                                                                   >  2:360904/1‑149 (MQ=255)                                     ctgGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACt                                                                                                             <  6:140114/149‑1 (MQ=255)                                     ctcgcacccccGCAAGGTTGGCTGGTCGCTCCCCTCCACGACACCTCTCTGCCCATCCCCTCCCCCGCCGCCCTGCACTACCGC‑GACACACCCCAGATCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCaaa                                                                                                               <  4:304947/140‑1 (MQ=255)                                      tgGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGTAAACt                                                                                                             >  3:297247/1‑148 (MQ=255)                                        gAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGAt                                                                                                          >  5:238423/1‑149 (MQ=255)                                        gAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGAt                                                                                                          <  1:360904/149‑1 (MQ=255)                                         aaaCCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATg                                                                                                         <  7:239573/149‑1 (MQ=255)                                            ccccGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTa                                                                                                      >  6:135457/1‑149 (MQ=255)                                             cccGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTAc                                                                                                     <  6:365106/149‑1 (MQ=255)                                               cGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTg                                                                                                   <  5:6204/149‑1 (MQ=255)                                                            gggcgggcGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACggg                                                                                      >  6:264068/1‑149 (MQ=255)                                                              gcgggcgCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGcc                                                                                    <  6:238423/149‑1 (MQ=255)                                                               cgggcgCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCt                                                                                   <  1:43133/149‑1 (MQ=255)                                                                    gcTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTAtgtg                                                                              <  5:135457/149‑1 (MQ=255)                                                                      ttCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGAACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCa                                                                            >  1:281790/1‑149 (MQ=255)                                                                         aaCGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCCCTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCt                                                                         >  6:21607/1‑149 (MQ=255)                                                                          aCGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTg                                                                        >  2:37170/1‑149 (MQ=255)                                                                          aCGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGCCCTATGTGCACCTg                                                                        >  8:115106/1‑149 (MQ=255)                                                                           cGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTAt                                                                                 <  1:271494/139‑1 (MQ=255)                                                                           cGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTAt                                                                                 >  2:271494/1‑139 (MQ=255)                                                                                 cTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCg                                                                                                                    <  8:151394/98‑1 (MQ=255)                                                                                 cTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCg                                                                                                                    >  7:151394/1‑98 (MQ=255)                                                                                    tcTGCGCATCCCCTACGCCGACGCCCTGCACTACCGA‑GAAAGACCCCACTTCGTGATTGCTCAGGAGTCCGGTGCC‑GCACTCCCGGTGCATGGCGGCAACCTGATGCTACTGGCGCTCCCACGGGCCTATGTGCACCTGCCCACCGCtt                                                              <  1:151154/148‑1 (MQ=255)                                                                                    gCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCGCACGGGCCTATGTGCACCTGCCCACCGCtt                                                              <  1:37170/149‑1 (MQ=255)                                                                                         ccATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCg                                                            >  2:20278/1‑146 (MQ=255)                                                                                         ccATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCg                                                            <  1:20278/146‑1 (MQ=255)                                                                                         ccATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTa                                                         <  7:318909/149‑1 (MQ=255)                                                                                           aTCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTAtgt                                                                               <  6:296052/125‑1 (MQ=255)                                                                                           aTCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTAtgt                                                                               >  5:296052/1‑125 (MQ=255)                                                                                              cccTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAAATGATCATACTGAAGCGCCCACGGGCAAATGTGCACCTGCCCACCCCTTCGCAACAgat                                                    >  8:176918/1‑147 (MQ=255)                                                                                               ccTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCt                                                    >  7:169688/1‑148 (MQ=255)                                                                                                  aCGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCa                                                                                          <  3:359489/107‑1 (MQ=255)                                                                                                  aCGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCa                                                                                          >  4:359489/1‑107 (MQ=255)                                                                                                  aCGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATc                                                >  2:212947/1‑149 (MQ=255)                                                                                                    gccgccgccCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGACAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTa                                              >  5:372898/1‑149 (MQ=255)                                                                                                      cgccgccCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATg                                            <  2:281790/149‑1 (MQ=255)                                                                                                            ccTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGc                                                                                                                  <  4:276899/73‑1 (MQ=255)                                                                                                            ccTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGc                                                                                                                  >  3:276899/1‑73 (MQ=255)                                                                                                            ccTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATgc                                      >  2:58822/1‑149 (MQ=255)                                                                                                            ccTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATgc                                      >  7:190269/1‑149 (MQ=255)                                                                                                              tGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATgcgc                                    >  3:64291/1‑149 (MQ=255)                                                                                                                cACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCt                                  >  5:327048/1‑149 (MQ=255)                                                                                                                  cTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTg                                 >  2:176403/1‑148 (MQ=255)                                                                                                                   tACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGc                                >  8:112027/1‑148 (MQ=255)                                                                                                                     ccGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCACACGGGCTTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCaa                              <  5:264068/148‑1 (MQ=255)                                                                                                                      cgc‑gAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGt                            >  3:188961/1‑149 (MQ=255)                                                                                                                            aaCACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGcccc                       <  4:64291/149‑1 (MQ=255)                                                                                                                                ccgcACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGa                   <  7:112027/146‑1 (MQ=255)                                                                                                                                  ccACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCgg                     <  3:55923/145‑1 (MQ=255)                                                                                                                                  ccACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCgg                     >  4:55923/1‑145 (MQ=255)                                                                                                                                   cACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCATCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACg                 >  2:256254/1‑148 (MQ=255)                                                                                                                                        cGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGTAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCaa             <  4:297247/147‑1 (MQ=255)                                                                                                                                          tGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGaaa         <  1:176403/149‑1 (MQ=255)                                                                                                                                            gcggCACCGGACTCCGGAGCCAGCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAAc        >  2:274201/3‑149 (MQ=255)                                                                                                                                                cTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGACAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAACTCAGc   <  6:372898/149‑1 (MQ=255)                                                                                                                                                 tCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAACCTCAGCt  >  3:322314/1‑149 (MQ=255)                                                                                                                                                 tCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAACTCAGCt  >  5:78535/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                      CACACGGCCGTACAGGAAGTGGCCGAGGAGTGCCTGCTGGAAACCCCGGAAGGTTGGCTGGGCGGGCGCTTCAACGACACCTGGCTGCCCATCCCCTACGCCGCCGCCCTGCACTACCGC‑GAAACACCCCACTTCGTGCTGGCTCCGGAGTCCGGTGCC‑GCACGCCCGGTGCATTGCGGCAAACTGATGCTACTGGAGCGCCCACGGGCCTATGTGCACCTGCCCACCGCTTCGCTACAGCTGATCTATGACATGCGCCTGCAAGTGCCCCGGGACGCCAAGAAACTCAGCT  >  NC_002947/5613380‑5613671

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