Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 832,299 T→C 6.5% L21P (CTG→CCG)  PP_0716 → HlyD family secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947832,2990TC 6.5% 183.5 / 2.6 64L21P (CTG→CCG) PP_0716HlyD family secretion protein
Reads supporting (aligned to +/- strand):  ref base T (28/31);  new base C (2/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

ACCCGCCAGCAAATGCAGCAACGCACCGACTGGGGCGAGTTGCTCGACGAGCCAGCCCCCACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTACGTAGTCCCGGTGGCCGCCGAAGTGGCCGGGCAGATCAAGCGTGTGGCCGTAG  >  NC_002947/832152‑832436
                                                                                                                                                   |                                                                                                                                          
aCCCGCCAGCAAATGCAGCAACGCACCGACTGGGGCGAGTTGCTAGACGAGCCAGCCCCCACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGtgctg                                                                                                                                           <  5:314760/149‑1 (MQ=255)
 cccGCCAGCAAATGCAGCAACGCACCGACTGGGGCGAGTTGCTCGACGAGCCAGCCCCCACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGtgctgc                                                                                                                                          >  6:313683/1‑149 (MQ=255)
      cAGCAAATGCAGCAACGCACCGACTGGGGCGAGTTGCTCGACGAGCCAGCCCCCACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTgctgctgct                                                                                                                                      >  6:4309/1‑148 (MQ=255)
      cAGCAAATGCAGCAACGCACCGACTGGGGCGAGTTGCTCGACGAGCCAGCCCCCACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGctgctgctg                                                                                                                                     <  7:202537/149‑1 (MQ=255)
           aaTGCAGCAACGCACCGACTGGGGCGAGTTGCTCGACGAGCCAGCCCCCACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGc                                                                                                                                >  7:127989/1‑149 (MQ=255)
            aTGCAGCAACGCACCGACTGGGGCGAGTTGCTCGACGAGCCAGCCCCCACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGcc                                                                                                                               <  6:234338/149‑1 (MQ=255)
             tGCAGCAACGCACCGACTGGGGCGAGTTGCTCGACGAGCCAGCCCCCACCGCGCATGCACAAGGTGAAACCGAATGAGAGACGCCACGCCAACGCCGCCCCCATCAGCCCAGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCg                                                                                                                              >  2:171303/1‑149 (MQ=255)
               cagcaACGCACCGACTGGGGCGAGTTGCTCGACGAGCCAGCCCCCACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTc                                                                                                                            <  8:123712/149‑1 (MQ=255)
                   aaCGCACCGACTGGGGCGAGTTGCTCGACGAGCCAGCCCCCACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGcc                                                                                                                        >  4:277812/1‑149 (MQ=255)
                   aaCGCACCGACTGGGGCGAGTTGCTCGACGAGCCAGCCCCCACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGcc                                                                                                                        <  4:197342/149‑1 (MQ=255)
                    aCGCACCGACTGGGGCGAGTTGCTCGACGAGCCAGCCCCCACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCt                                                                                                                       >  1:64287/1‑149 (MQ=255)
                    aCGCACCGACTGGGGCGAGTTGCTCGACGAGCCAGCCCCCACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCt                                                                                                                       >  3:145447/1‑149 (MQ=255)
                       cACCGACTGGGGCGAGTTGCTCGACGAGCCAGCCCCCACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGct                                                                                                                    >  5:193875/1‑149 (MQ=255)
                        aCCGACTGGGGCGAGTTGCTCGACGAGCCAGCCCCCACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGctc                                                                                                                   >  7:206843/1‑149 (MQ=255)
                         ccGACTGGGGCGAGTTGCTCGACGAGCCAGCCCCCACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGctct                                                                                                                  >  4:289886/1‑149 (MQ=255)
                             cTGGGCCGAGTTGCTCGACGAGCCAGCCCCCACCGCGCATGCACAAGGGGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGTGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTa                                                                                                              <  7:297743/149‑1 (MQ=255)
                                gggCGAGTTGCTCGACGAGCCAGCCCCCACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACtt                                                                                                           <  5:313683/149‑1 (MQ=255)
                                    gAGTTGCTCGACGAGCCAGCCCCCACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTgctg                                                                                                       <  5:140010/149‑1 (MQ=255)
                                          cTCGACGAGCCAGCCCCCACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGAc                                                                                                 >  7:253431/1‑149 (MQ=255)
                                          cTCGACGAGCCAGCCCCCACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGAc                                                                                                 >  8:2776/1‑149 (MQ=255)
                                             gacgaGCCAGCCCCCACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCg                                                                                               <  5:4309/148‑1 (MQ=255)
                                                   ccagccCCCACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACg                                                                                        >  4:117833/1‑149 (MQ=255)
                                                        ccccACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTa                                                                                   <  2:64287/149‑1 (MQ=255)
                                                         cccACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTa                                                                                   >  7:177172/1‑148 (MQ=255)
                                                          ccACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTaca                                                                                 >  6:331189/1‑149 (MQ=255)
                                                                 gcATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCCGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGc                                                                           >  3:81060/1‑148 (MQ=255)
                                                                 gcATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCCGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGc                                                                           >  5:42493/1‑148 (MQ=255)
                                                                  cATGCGCAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCaa                                                                         <  2:10115/149‑1 (MQ=255)
                                                                  cATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGc                                                                           >  7:231043/1‑147 (MQ=255)
                                                                   aTGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATATGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAg                                                                        <  5:59534/149‑1 (MQ=255)
                                                                      cacaAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑c                                                                    <  3:289886/149‑1 (MQ=255)
                                                                      cacaAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑c                                                                    <  3:117833/149‑1 (MQ=255)
                                                                        caAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGc                                                                  <  6:48920/149‑1 (MQ=255)
                                                                         aaGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGcc                                                                 >  1:153800/1‑149 (MQ=255)
                                                                         aaGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGcc                                                                 <  5:331189/149‑1 (MQ=255)
                                                                          aGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCt                                                                >  4:4/1‑149 (MQ=255)
                                                                          aGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCt                                                                <  7:73222/149‑1 (MQ=255)
                                                                          aGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCCGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCt                                                                <  4:81060/149‑1 (MQ=255)
                                                                          aGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCCGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCt                                                                <  6:42493/149‑1 (MQ=255)
                                                                           ggTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTg                                                               >  2:159753/1‑149 (MQ=255)
                                                                             tGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTg                                                               >  7:221439/1‑147 (MQ=255)
                                                                             tGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCa                                                             <  8:127989/149‑1 (MQ=255)
                                                                               aaaCCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAgg                                                           >  1:276754/1‑149 (MQ=255)
                                                                                    gAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTAc                                                      <  8:206843/149‑1 (MQ=255)
                                                                                    gAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTAc                                                      <  1:159753/149‑1 (MQ=255)
                                                                                    gAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTAc                                                      <  2:33067/149‑1 (MQ=255)
                                                                                     aaTGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTACg                                                     >  8:267315/1‑149 (MQ=255)
                                                                                     aaTGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTACg                                                     >  2:71989/1‑149 (MQ=255)
                                                                                           cGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTACGTAGTcc                                               <  8:221439/149‑1 (MQ=255)
                                                                                             acgccacgccAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGctct                                                                                                                  >  2:42161/1‑81 (MQ=255)
                                                                                             acgccacgccAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGctct                                                                                                                  <  1:42161/81‑1 (MQ=255)
                                                                                             acgccacgccAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTACGTAGTCCCg                                             <  8:231043/149‑1 (MQ=255)
                                                                                              cgccacgccaACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTACGTAGTCCCgg                                            <  3:4/149‑1 (MQ=255)
                                                                                                ccacgccaACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTACGTAGTCCCGGTg                                          <  4:145447/149‑1 (MQ=255)
                                                                                                ccacgccaACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGGCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTACGTAGTCCCGGTg                                          <  6:190787/149‑1 (MQ=255)
                                                                                                 cacgccaACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTACGTAGTCCCGGTgg                                         >  2:263706/1‑149 (MQ=255)
                                                                                                         cgccgcCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTACGTAGTCCCGGTGGCCGCCGaa                                 <  2:153800/149‑1 (MQ=255)
                                                                                                         cgccgcCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTACGTAGTCCCGGTGGCCGCCGaa                                 <  8:253431/149‑1 (MQ=255)
                                                                                                          gccgccCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTACGTAGTCCCGGTGGCCGCCGAAg                                <  7:267315/149‑1 (MQ=255)
                                                                                                                 ccATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTACGTAGTCCCGGTGGCCGCCGAAGTGGCCgg                         >  4:132269/1‑149 (MQ=255)
                                                                                                                 ccATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTACGTAGTCCCGGTGGCCGCCGAAGTGGCCgg                         <  6:193875/149‑1 (MQ=255)
                                                                                                                      aGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCTCCGCCTGCAGGCTTACGTAGTCCCGGTGGCCGCCGAAGTGGCCGGGCa                      <  1:71989/148‑1 (MQ=255)
                                                                                                                        ccccGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTACGTAGTCCCGGTGGCCGCCGAAGTTGCCGGGCAGAt                   >  4:72304/1‑148 (MQ=255)
                                                                                                                        ccccGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTACGTAGTCCCGGTGGCCGCCGAAGTGGCCGGGCAGAt                   >  6:40574/1‑148 (MQ=255)
                                                                                                                                 gcgGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTACGTAGTCCCGGTGGCCGCCGAAGTGGCCGGGCAGATCAAGCgtgt          <  3:132269/148‑1 (MQ=255)
                                                                                                                                      aTACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTACGTAGTCCCGGTGGCCGCCGAAGTGGCCGGGCAGATCAAGCGTGTGGCCg     >  3:37010/1‑148 (MQ=255)
                                                                                                                                      aTACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTAAGTAGTACCCGTGGCCGCCGAAGTCGCCGGGCAGATAAAGCGTGTGGCCg     >  4:158334/1‑148 (MQ=255)
                                                                                                                                         cGCTGGGTGCAGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTACGTAGTCCCGGTGGCCGCCGAAGTGGCCGGGCAGATCAAGCGTGTGGCCGTAg  <  6:173630/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                          
ACCCGCCAGCAAATGCAGCAACGCACCGACTGGGGCGAGTTGCTCGACGAGCCAGCCCCCACCGCGCATGCACAAGGTGAAACCGAATGAGCGACGCCACGCCAACGCCGCCCCCATCAGCCCCGGATCGCGGCATACGCTGGGTGCTGCTGCTGATTGCCGTCAGCCTGCTCTGGTACTTGCTGGCTGACCGCTATACGCCGTACACCCAGCAAGCT‑CGCCTGCAGGCTTACGTAGTCCCGGTGGCCGCCGAAGTGGCCGGGCAGATCAAGCGTGTGGCCGTAG  >  NC_002947/832152‑832436

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