Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,443,376:1 +C 6.0% coding (1452/1989 nt) PP_4781 ← sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,443,3761.C 6.0% 341.1 / 14.9 84coding (1452/1989 nt)PP_4781sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base . (41/38);  new base C (3/2);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01

TTCTGCCCGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAATCACTCCGTT  >  NC_002947/5443227‑5443525
                                                                                                                                                      |                                                                                                                                                     
ttCTGCCCGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAgg                                                                                                                                                        >  1:266321/1‑150 (MQ=255)
ttCTGCCCGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAgg                                                                                                                                                        >  1:1308584/1‑150 (MQ=255)
 tCTGCCCGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAgg‑c                                                                                                                                                      <  1:701510/150‑1 (MQ=255)
     cccGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAgg‑c                                                                                                                                                      <  1:1175207/146‑1 (MQ=255)
     cccGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAgg‑c                                                                                                                                                      >  2:1175207/1‑146 (MQ=255)
      ccGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑cggccg                                                                                                                                                 <  2:266321/150‑1 (MQ=255)
        gTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑cggccggc                                                                                                                                               >  2:526272/1‑150 (MQ=255)
        gTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑cggccggc                                                                                                                                               <  1:620370/150‑1 (MQ=255)
        gTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑cggccggc                                                                                                                                               <  1:1655420/150‑1 (MQ=255)
           aCCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTg                                                                                                                                            >  1:1712134/1‑150 (MQ=255)
             cGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCa                                                                                                                                          >  2:1241989/1‑150 (MQ=255)
                gCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAgg‑cg                                                                                                                                                     <  1:612296/136‑1 (MQ=255)
                gCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAgg‑cg                                                                                                                                                     >  2:612296/1‑136 (MQ=255)
                       gcccGTCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAg                                                                                                                                >  2:1425986/3‑150 (MQ=255)
                        tCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAgg                                                                                                                               <  2:1632607/150‑1 (MQ=255)
                         ccGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTTCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAggg                                                                                                                              <  2:1374541/150‑1 (MQ=255)
                           ggCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTa                                                                                                                            >  2:1432994/1‑150 (MQ=255)
                                gTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGT‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAAt                                                                                                                       <  1:1240639/150‑1 (MQ=255)
                                 tCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATg                                                                                                                      <  2:1308584/150‑1 (MQ=255)
                                 tCGGCCGGAATGTCCAGGCTCAGGAACAGGCCTTTGCCCAGCGCCGGCCCGGCACAACCAGCGCGCAAATCCTGCAACCGGCCGCCCAGGCTAAACGGGGCCGCCGGCGCACGCAGG‑CGGCCGGCCTGCAGAGCGAAGAGGGTAaaaagg                                                                                                                      >  2:802769/1‑148 (MQ=255)
                                   ggCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGTGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGCGTAAGAATggg                                                                                                                    <  2:318052/150‑2 (MQ=255)
                                               cAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATc                                                                                                        >  1:1527252/1‑150 (MQ=255)
                                               cAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATc                                                                                                        >  1:880323/1‑150 (MQ=255)
                                                aGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATcc                                                                                                       >  2:1251688/1‑150 (MQ=255)
                                                aGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATcc                                                                                                       >  2:545495/1‑150 (MQ=255)
                                                aGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCAGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATcc                                                                                                       >  1:1260388/1‑150 (MQ=255)
                                                 ggCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCca                                                                                                      <  2:447585/150‑1 (MQ=255)
                                                 ggCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCca                                                                                                      <  1:690878/150‑1 (MQ=255)
                                                   cTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCTTGCAACAGGGCGCACAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCcatc                                                                                                    <  2:1038603/150‑1 (MQ=255)
                                                          aCAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCtt                                                                                             >  1:667225/1‑150 (MQ=255)
                                                              gCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCACAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATc                                                                                                                 <  2:610454/126‑1 (MQ=255)
                                                              gCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCACAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATc                                                                                                                 >  1:610454/1‑126 (MQ=255)
                                                               ccTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCGCCCGCTGGGGCGGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGACGGCCTGCAGTTCGGACAGGGTAAGAATGGCCTCGACCATATCCATCATGCCTTGGGCa                                                                                        <  2:880323/150‑1 (MQ=255)
                                                                cTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAg                                                                                       >  2:1206144/1‑150 (MQ=255)
                                                                 tttGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTg                                                                                                                                            >  2:1730242/1‑96 (MQ=255)
                                                                 tttGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTg                                                                                                                                            <  1:1730242/96‑1 (MQ=255)
                                                                  ttGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCTCGCAGG‑CGGCCGGCCTCAATTGCGGACAGGGGAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGAc                                                                                     <  1:1179547/150‑1 (MQ=255)
                                                                  ttGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTTGGCAGAc                                                                                     <  1:1206144/150‑1 (MQ=255)
                                                                      ccAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCca                                                                                 >  2:956958/1‑150 (MQ=255)
                                                                      ccAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCca                                                                                 >  1:1678410/1‑150 (MQ=255)
                                                                      ccAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCca                                                                                 <  1:1358512/150‑1 (MQ=255)
                                                                       cAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCcac                                                                                >  2:292536/1‑150 (MQ=255)
                                                                       cAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGAACACTGAACGGTGCCGGCTACGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCcac                                                                                >  1:987510/1‑150 (MQ=255)
                                                                             cTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCgt                                                                          <  1:707838/150‑1 (MQ=255)
                                                                              tGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCgtg                                                                         <  2:2647/150‑1 (MQ=255)
                                                                                 ccGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGAcc                                                                                                             >  2:1287135/1‑111 (MQ=255)
                                                                                 ccGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGAcc                                                                                                             <  1:1287135/111‑1 (MQ=255)
                                                                                 ccGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGa                                                                      <  1:1282316/150‑1 (MQ=255)
                                                                                  cGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTggg                                                                                          >  2:26230/1‑129 (MQ=255)
                                                                                  cGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTggg                                                                                          <  1:26230/129‑1 (MQ=255)
                                                                                     cATAACCGGCGCGCAAACTCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGt                                                                  <  1:96415/150‑1 (MQ=255)
                                                                                      aTAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTa                                                                 <  1:526272/150‑1 (MQ=255)
                                                                                          ccGGCGCGCAAACCCTGAAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGa                                                             <  1:1677561/150‑1 (MQ=255)
                                                                                            ggCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGa                                                                      >  2:1288298/1‑139 (MQ=255)
                                                                                            ggCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGa                                                                      <  1:1288298/139‑1 (MQ=255)
                                                                                                 actgaaccTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCAttt                                                      >  2:1380154/7‑150 (MQ=255)
                                                                                                   aaaCCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCa                                                    <  2:1442730/150‑1 (MQ=255)
                                                                                                       ccTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGccc                                                >  2:1829483/1‑150 (MQ=255)
                                                                                                       ccTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGccc                                                >  1:489868/1‑150 (MQ=255)
                                                                                                       ccTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGccc                                                >  1:166054/1‑150 (MQ=255)
                                                                                                       ccTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGAACATAGGACAGGGTAACGACGGCGTCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGccc                                                >  1:1176959/1‑150 (MQ=255)
                                                                                                         tGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGccccc                                              <  1:1425986/150‑1 (MQ=255)
                                                                                                         tGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGccccc                                              >  2:1717455/1‑150 (MQ=255)
                                                                                                         tGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCAAAGCCGTGTGATGGTACTGAGCCATTTCAGccccc                                              >  1:771893/1‑150 (MQ=255)
                                                                                                                  gcgcgcAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAgg                                     <  2:676366/150‑1 (MQ=255)
                                                                                                                  gcgcgcAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAgg                                     <  2:232073/150‑1 (MQ=255)
                                                                                                                     cgcAGGCTGAACGGTGCCGGCTGCGCCTGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGgtg                                  >  2:437208/1‑150 (MQ=255)
                                                                                                                       cAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGc                                                  >  1:184202/1‑132 (MQ=255)
                                                                                                                       cAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGc                                                  <  2:184202/132‑1 (MQ=255)
                                                                                                                         ggCTGAACGGTGCCGGCTGCGCGCGCCGG‑CGGCCGGCCTGCAGTTCGGAGAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCa                              <  2:576996/150‑1 (MQ=255)
                                                                                                                          gCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCAt                             >  2:1396951/1‑150 (MQ=255)
                                                                                                                          gCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCAt                             >  2:1396897/1‑150 (MQ=255)
                                                                                                                           cTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATc                            <  2:1305910/150‑1 (MQ=255)
                                                                                                                            tGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCa                           <  1:783053/150‑1 (MQ=255)
                                                                                                                            tGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCa                           <  2:113745/150‑1 (MQ=255)
                                                                                                                                   tGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTcgag                    >  1:1524076/1‑150 (MQ=255)
                                                                                                                                      cGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTcgagcga                 >  1:31684/1‑150 (MQ=255)
                                                                                                                                      cGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTcgagcga                 <  2:489868/150‑1 (MQ=255)
                                                                                                                                       ggCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGgtgt                                 >  2:1896361/1‑133 (MQ=255)
                                                                                                                                       ggCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGgtgt                                 <  1:1896361/133‑1 (MQ=255)
                                                                                                                                       ggCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTcgagcgag                >  2:547409/1‑150 (MQ=255)
                                                                                                                                         cTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGcc              <  2:1176959/150‑1 (MQ=255)
                                                                                                                                           gcgcCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACt                                                               >  2:91465/1‑99 (MQ=255)
                                                                                                                                           gcgcCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACt                                                               <  1:91465/99‑1 (MQ=255)
                                                                                                                                            gaactgCAGGCCGGCC‑GCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAAt           >  2:1720770/6‑150 (MQ=255)
                                                                                                                                            cgcCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAAt           <  1:1174235/150‑1 (MQ=255)
                                                                                                                                               ctgCAGGCCGGCC‑GCCTGCGGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCgt                                                                          <  2:1112534/82‑1 (MQ=255)
                                                                                                                                               ctgCAGGCCGGCC‑GCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTcgagcga                 <  2:394950/139‑1 (MQ=255)
                                                                                                                                               ctgCAGGCCGGCC‑GCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTcgagcga                 >  1:394950/3‑141 (MQ=255)
                                                                                                                                               ctgCAGG‑CCGGCCGCCTGCGGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCgt                                                                          >  1:1112534/3‑84 (MQ=255)
                                                                                                                                                tgCAGGCCGGCC‑GCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAATCACt       >  1:767703/2‑150 (MQ=255)
                                                                                                                                                cGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAATCACt       <  1:1713384/150‑1 (MQ=255)
                                                                                                                                                 gCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAATCACTc      >  2:1510745/1‑150 (MQ=255)
                                                                                                                                                 gCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGCTGTGCATCAGCTCGACCGAGCCAACCACTc      >  2:330050/1‑150 (MQ=255)
                                                                                                                                                   aGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAATCACTCCg    >  2:992188/1‑150 (MQ=255)
                                                                                                                                                     g‑cggcCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAATCACTCCGtt  <  2:31684/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TTCTGCCCGTCACCGAGCAACCCGTCCGGCAGGTCGGCCGGAATGTCCAGGCTCAGGTACAGGCCTTTGCCCAGCGCCTGCCCGGCATAACCGGCGCGCAAACCCTGCAACAGGGCGCGCAGGCTGAACGGTGCCGGCTGCGCCCGCAGG‑CGGCCGGCCTGCAGTTCGGACAGGGTAAGAATGGCATCGACCATATCCATCATGCCTTGGGCAGACCCCACAGCCGTGTGATGGTACTGAGCCATTTCAGCCCCCATCGGCAGGGTGTGCATCAGCTCGAGCGAGCCAATCACTCCGTT  >  NC_002947/5443227‑5443525

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑