Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,118,628 C→G 6.0% Y134* (TAC→TAG)  PP_4505 → sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,118,6280CG 6.0% 218.5 / 4.1 67Y134* (TAC→TAG) PP_4505sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (30/33);  new base G (2/2);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

GCGGTGAAGCCGTGGTGCCCGGTACCTTCCCACGCCGTGCGGCAGCCCTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCAGCCTGCTGGTGACGCTGTCCATTACCCGCCCGCTAAGCCGCCTGCGCAGTG  >  NC_002947/5118480‑5118776
                                                                                                                                                    |                                                                                                                                                    
gCGGTGAAGCCGTGGTGCCCGGTACCTTCCCACGCCGTGCGGCAGCCCTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTAc                                                                                                                                                      >  3:369357/1‑149 (MQ=255)
      aaGCCGTGGTGCCCGGTACCTTCCCACGCCGTGCGGCAGCCCTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCAt                                                                                                                                                 <  2:94756/148‑1 (MQ=255)
         ccGTGGTGCCCGGTACCTTCCCACGCCGTGCGGCAGCCCTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATccc                                                                                                                                              >  6:101383/1‑148 (MQ=255)
          cGTGGTGCCCGGTACCTTCCCACGCCGTGCGGCAGCCCTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATcc                                                                                                                                               >  4:261515/1‑146 (MQ=255)
          cGTGGTGCCCGGTACCTTCCCACGCCGTGCGGCAGCCCTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATcc                                                                                                                                               <  3:261515/146‑1 (MQ=255)
            tggtgCCCGGTACCTTCCCACGCCGTGCGGCAGCCCTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATcc                                                                                                                                               <  6:312784/144‑1 (MQ=255)
            tggtgCCCGGTACCTTCCCACGCCGTGCGGCAGCCCTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATcc                                                                                                                                               >  5:312784/1‑144 (MQ=255)
                   cGGTACCTTCCCACGCCGTGCGGCAGCCCTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCg                                                                                                                                                    <  3:57424/132‑1 (MQ=255)
                   cGGTACCTTCCCACGCCGTGCGGCAGCCCTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCg                                                                                                                                                    >  4:57424/1‑132 (MQ=255)
                      tACCTTCCCACGCCGTGCGGCAGCCCTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCTCTgg                                                                                                                                <  6:65488/149‑1 (MQ=255)
                      tACCTTCCCACGCCGTGCGGCAGCCCTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTgg                                                                                                                                <  5:53671/149‑1 (MQ=255)
                        ccTTCCCACGCCGTGCGGCAGCCCTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGAt                                                                                                                              <  8:148543/149‑1 (MQ=255)
                             ccACGCCGTGCGGCAGCCCTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTg                                                                                                                         >  2:13557/1‑149 (MQ=255)
                              cACGCCGTGCGGCAGCCCTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTgg                                                                                                                        <  4:369357/149‑1 (MQ=255)
                                           aGCCCTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCga                                                                                                                <  1:323561/144‑1 (MQ=255)
                                           aGCCCTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCga                                                                                                                >  2:323561/1‑144 (MQ=255)
                                              ccTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTAGCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCt                                                                                                          <  6:307951/147‑1 (MQ=255)
                                              ccTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTAGCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCt                                                                                                          >  3:350348/1‑147 (MQ=255)
                                              ccTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTAGCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCt                                                                                                          <  4:350348/147‑1 (MQ=255)
                                              ccTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTAGCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCt                                                                                                          >  5:307951/1‑147 (MQ=255)
                                              ccTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGc                                                                                                        >  2:23046/1‑149 (MQ=255)
                                              ccTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGc                                                                                                        >  6:304416/1‑149 (MQ=255)
                                                  gAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCttt                                                                                                     <  2:58772/148‑1 (MQ=255)
                                                      cgcgcCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCg                                                                                                >  2:224149/1‑149 (MQ=255)
                                                       gcgcCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGc                                                                                               >  5:270156/1‑149 (MQ=255)
                                                        cgcCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCt                                                                                              <  5:304416/149‑1 (MQ=255)
                                                         gcCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCt                                                                                              >  7:248642/1‑148 (MQ=255)
                                                            aGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAgc                                                                                          <  2:29972/149‑1 (MQ=255)
                                                             gCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAgcg                                                                                         <  8:248642/149‑1 (MQ=255)
                                                               acaATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCa                                                                                       >  3:209292/1‑149 (MQ=255)
                                                                   tGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTgg                                                                                   >  5:134553/1‑149 (MQ=255)
                                                                       cAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCAt                                                                               >  3:144342/1‑149 (MQ=255)
                                                                        aGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGTTTTGGCCGCTCAGCGCACTGGGCAtt                                                                              <  3:9101/149‑1 (MQ=255)
                                                                         gcggcgCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCTCTGGATGCCTGGCACAGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTg                                                                             >  7:223572/1‑149 (MQ=255)
                                                                            gcgcTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCt                                                                          >  2:242395/1‑149 (MQ=255)
                                                                              gcTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTg                                                                        <  5:101383/149‑1 (MQ=255)
                                                                               cTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTgg                                                                       <  1:13557/149‑1 (MQ=255)
                                                                                  gCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCAc                                                                                                                     <  5:314294/100‑1 (MQ=255)
                                                                                  gCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCAc                                                                                                                     >  6:314294/1‑100 (MQ=255)
                                                                                   ccATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCga                                                                                                                >  1:310662/1‑104 (MQ=255)
                                                                                   ccATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCga                                                                                                                <  2:310662/104‑1 (MQ=255)
                                                                                   ccATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGt                                                                   <  6:134553/149‑1 (MQ=255)
                                                                                     aTGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGc                                                                 >  7:12112/1‑149 (MQ=255)
                                                                                        gcgcCGGCTGACCGAGGAATACACCAGCCCCGACACGGGCGAAACCTACTTGCTGATCTACCGCTTCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGa                                                              <  3:331477/149‑1 (MQ=255)
                                                                                        gcgcCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTg                                                               <  1:11222/148‑1 (MQ=255)
                                                                                                tGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGtt                                                       >  3:5925/1‑148 (MQ=255)
                                                                                                  aCCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCa                                                     <  4:144342/148‑1 (MQ=255)
                                                                                                   ccGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCAGc                                                   <  1:242395/149‑1 (MQ=255)
                                                                                                    cGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTg                                                                                                                                 <  4:342170/70‑1 (MQ=255)
                                                                                                    cGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTg                                                                                                                                 >  3:342170/1‑70 (MQ=255)
                                                                                                    cGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCa                                                     >  3:309533/1‑146 (MQ=255)
                                                                                                    cGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCa                                                     <  4:309533/146‑1 (MQ=255)
                                                                                                     gAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCAGCct                                                 <  1:224149/149‑1 (MQ=255)
                                                                                                      aGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCAGCctg                                                >  8:215603/1‑149 (MQ=255)
                                                                                                       ggAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCAGCctgc                                               <  6:270156/149‑1 (MQ=255)
                                                                                                       ggAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCAGCctgc                                               <  1:23046/149‑1 (MQ=255)
                                                                                                         aaTACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCAGCctgctg                                             <  8:122141/149‑1 (MQ=255)
                                                                                                          aTACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCAGCCTGCTgg                                            >  6:165168/1‑149 (MQ=255)
                                                                                                                cAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCa                                                     >  8:253749/1‑134 (MQ=255)
                                                                                                                cAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCa                                                     <  7:253749/134‑1 (MQ=255)
                                                                                                                cAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCAGCCTGCTGGTGACGc                                      >  3:8270/1‑149 (MQ=255)
                                                                                                                          acCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCAGCCTGCTGGTGACGCTGTCCATTAc                            <  4:5925/149‑1 (MQ=255)
                                                                                                                             ggCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCAGCCTGCTGGTGACGCTGTCCATTAcccg                         <  2:49211/149‑1 (MQ=255)
                                                                                                                                     cTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCAGCCTGCTGGTGACGCTGTCCATTACCCGCCCGCt                   >  6:96453/1‑147 (MQ=255)
                                                                                                                                         ttGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCAGCCTGCTGGTGa                                         <  8:330449/121‑1 (MQ=255)
                                                                                                                                         ttGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCAGCCTGCTGGTGa                                         >  7:330449/1‑121 (MQ=255)
                                                                                                                                           gCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCAGCCTGCTGGTGACGCTGTCCATTACCCGCCCGCTAAgccgcc           >  6:383692/1‑149 (MQ=255)
                                                                                                                                                   aCCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCAGCCTGCTGGTGACGCTGTCCATTACCCGCCCGCTAAGCCGCCTGCGCAGt   >  5:197404/1‑149 (MQ=255)
                                                                                                                                                    ccGCATCCCTCACCCGGCGCTGGATGGCTGGCACCGCGAGAGCCTGCTTTGGCCGCGCAGCGCACTGTGCATTGCCTTGGTCGTGCTGACCATGTTCAGCCTGCTGTTGACGCTGTCCATTACCCGCCCGATAAGCCGCCTGCGCa     <  8:223572/146‑1 (MQ=255)
                                                                                                                                                    ccGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCAGCCTGCTGGTGACGCTGTCCATTACCCGCCCGCTAAGCCGCCTGCGCAGTg  <  7:215603/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGGTGAAGCCGTGGTGCCCGGTACCTTCCCACGCCGTGCGGCAGCCCTCGAAGCGCGCCAGCACAATGACCAGCGGCGCTTGCCATGGCGCCGGCTGACCGAGGAATACACCAGCCCCGACACCGGCGAAACCTACTTGCTGATCTACCGCATCCCTCACCCGGCGCTGGATGCCTGGCACCGCGAGAGCCTGCTTTGGCCGCTCAGCGCACTGGGCATTGCCTTGGTCGTGCTGACCCTGTTCAGCCTGCTGGTGACGCTGTCCATTACCCGCCCGCTAAGCCGCCTGCGCAGTG  >  NC_002947/5118480‑5118776

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