Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,030,972 C→G 6.8% P343A (CCG→GCG)  PP_0887 → sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,030,9720CG 6.8% 421.0 / 15.6 133P343A (CCG→GCG) PP_0887sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (54/70);  new base G (4/5);  total (58/75)
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.33e-01

Only 100 of 137 total aligned reads displayed.

GCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTGGGACGCCCAGGACATCGTCATCAGTATCCACGACCATGGTCCTGGTGTACCGCCAG  >  NC_002947/1030823‑1031116
                                                                                                                                                     |                                                                                                                                                 
gcgcAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGgcgc                                                                                                                                                   <  1:840724/150‑1 (MQ=255)
     gCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGcgccgccg                                                                                                                                              >  1:1693316/1‑150 (MQ=255)
      ccGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGcgccgccgc                                                                                                                                             >  1:1867612/1‑150 (MQ=255)
      ccGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGcgccgccgc                                                                                                                                             <  2:768412/150‑1 (MQ=255)
        gAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCa                                                                                                                                           >  1:2383876/1‑150 (MQ=255)
          aGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGGCGGCGCCAGa                                                                                                                                         <  1:1728325/150‑1 (MQ=255)
           gCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGAc                                                                                                                                        <  1:572271/150‑1 (MQ=255)
             cAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGGCGGCGCCAGAcgc                                                                                                                                      <  1:1143781/150‑3 (MQ=255)
             cAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGGCGGCGCCAGAcgc                                                                                                                                      <  1:775080/150‑3 (MQ=255)
              aaCCGGCGCCTGGCCATCGTGGAGCAGGATGTGAGCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTg                                                                                                                                     >  1:1377943/1‑150 (MQ=255)
              aaCCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTATTGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTg                                                                                                                                     >  1:626520/1‑150 (MQ=255)
                ccGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGAc                                                                                                                                  <  2:573045/150‑1 (MQ=255)
                       cTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCg                                                                                                                           >  1:84677/1‑150 (MQ=255)
                        tGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGt                                                                                                                          >  2:1642262/1‑150 (MQ=255)
                          gCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTa                                                                                                                        >  2:828298/1‑150 (MQ=255)
                           ccATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGAcc                                                                                                                                 >  1:1339123/1‑140 (MQ=255)
                           ccATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGAcc                                                                                                                                 <  2:1339123/140‑1 (MQ=255)
                           ccATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTAc                                                                                                                       <  1:1624140/150‑1 (MQ=255)
                           ccATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTAc                                                                                                                       <  2:334087/150‑1 (MQ=255)
                           ccATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTAc                                                                                                                       >  1:2095598/1‑150 (MQ=255)
                            cATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACGAAATCGTTACt                                                                                                                      >  1:1282679/1‑150 (MQ=255)
                            cATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACt                                                                                                                      <  1:1431708/150‑1 (MQ=255)
                                 tGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATc                                                                                                                 >  1:713255/1‑150 (MQ=255)
                                    aGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCGAACCTGACCCAGGCGTTACTGAATctgc                                                                                                              <  1:1881740/150‑1 (MQ=255)
                                    aGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATctgc                                                                                                              <  2:709564/150‑1 (MQ=255)
                                      cAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATctgctg                                                                                                            >  1:199514/1‑150 (MQ=255)
                                         gATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAAc                                                                                                         >  2:2045517/1‑150 (MQ=255)
                                          aTGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGa                                                                                                                                   <  2:2234912/123‑1 (MQ=255)
                                          aTGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGa                                                                                                                                   >  1:2234912/1‑123 (MQ=255)
                                          aTGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACa                                                                                                        >  1:2081949/1‑150 (MQ=255)
                                            gtgACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGGAGGCTGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAAt                                                                                                      <  2:1640138/150‑1 (MQ=255)
                                              gcccGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATgc                                                                                                    >  2:593548/3‑150 (MQ=255)
                                              gACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATgc                                                                                                    <  2:2292559/150‑1 (MQ=255)
                                               aCCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGGCGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATgcc                                                                                                   >  1:1693728/1‑150 (MQ=255)
                                               aCCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATgcc                                                                                                   <  2:1900938/150‑1 (MQ=255)
                                                    ctggctggACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCTCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGa                                                                                              >  2:1259175/1‑150 (MQ=255)
                                                     tggctggACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCTTAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGAt                                                                                             <  1:624517/150‑1 (MQ=255)
                                                       gctggACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGAtgc                                                                                           <  2:1170088/150‑1 (MQ=255)
                                                       gctggACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAACCAGCTACCTTTTCCTGAGCATGCCCGACCGCAACGTGCCGCTCCTGGCCCCGCCGCCAGACCCGACCCAGGCGTTACTGACTCTGCTGAACAATGCCGCCGAtgc                                                                                           <  2:1084940/150‑1 (MQ=255)
                                                          ggACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGAtgcctg                                                                                        >  2:583667/1‑150 (MQ=255)
                                                              gAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCt                                                                                    >  2:1065636/1‑150 (MQ=255)
                                                               aGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCtg                                                                                   <  2:1867612/150‑1 (MQ=255)
                                                                ggCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACt                                                                                                                      <  2:2092046/114‑1 (MQ=255)
                                                                ggCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATgc                                                                                                    <  2:1975020/132‑1 (MQ=255)
                                                                    cTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACTCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTgatgat                                                                              <  2:2403199/150‑1 (MQ=255)
                                                                    cTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTgatgat                                                                              <  2:149597/150‑1 (MQ=255)
                                                                     tCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATc                                                                             <  2:84677/150‑1 (MQ=255)
                                                                     tCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATc                                                                             >  1:817631/1‑150 (MQ=255)
                                                                      cAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCt                                                                            >  1:2073081/1‑150 (MQ=255)
                                                                      cAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCt                                                                            <  1:1853266/150‑1 (MQ=255)
                                                                        aCCGCTGGCACCTGATGCGCCGGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTgg                                                                          <  2:859959/150‑1 (MQ=255)
                                                                        aCCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTgg                                                                          <  2:1857780/150‑1 (MQ=255)
                                                                          cGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGaa                                                                        >  2:629237/1‑150 (MQ=255)
                                                                            cTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACt                                                                                                                      >  1:619041/1‑102 (MQ=255)
                                                                            cTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGt                                                                      <  1:2240704/150‑1 (MQ=255)
                                                                             gatgtccTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGtt                                                                     >  2:396626/6‑150 (MQ=255)
                                                                                cACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGc                                                                  <  2:817632/150‑1 (MQ=255)
                                                                                cACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGc                                                                  <  2:817631/150‑1 (MQ=255)
                                                                                 aCCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGcc                                                                 >  1:701347/1‑150 (MQ=255)
                                                                                 aCCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGcc                                                                 <  2:1693316/150‑1 (MQ=255)
                                                                                     gATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGCCCTCCCCCAGGCGTTACTGAAGCTGCTGAACCATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGa                                                             <  2:1586533/150‑1 (MQ=255)
                                                                                     gATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGa                                                             <  2:2073081/150‑1 (MQ=255)
                                                                                     gATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGa                                                             >  1:2356363/1‑150 (MQ=255)
                                                                                     gATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACACAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGa                                                             <  2:1377943/150‑1 (MQ=255)
                                                                                         cgcCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTgg                                                         <  1:583667/150‑1 (MQ=255)
                                                                                           cccGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTAc                                                                                                                       <  2:1320855/86‑1 (MQ=255)
                                                                                           cccGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTAc                                                                                                                       >  1:1320855/1‑86 (MQ=255)
                                                                                             cGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCtgat                                                                                 >  1:389179/1‑122 (MQ=255)
                                                                                              ggAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTGGGACGc                                                    >  1:396870/1‑150 (MQ=255)
                                                                                                 aGCCAGCTACCGCTTCAAGCGCCTGCGCGGCGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCTTGGACTGGAACGCCCa                                                 >  2:1223389/1‑150 (MQ=255)
                                                                                                  gCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTGGGACGCCCAg                                                <  2:396870/150‑1 (MQ=255)
                                                                                                   ccAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCtga                                                                                  <  2:1068695/115‑1 (MQ=255)
                                                                                                   ccAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCtga                                                                                  >  1:1068695/1‑115 (MQ=255)
                                                                                                   ccAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGACAAGGCTTGCTCCGCCGATACACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTGGGACGCCCAgg                                               >  2:683820/1‑150 (MQ=255)
                                                                                                      gCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCGGCCGATGCCTGCCCTGATGATCGGGAAGTTCGCCTGGACTGGGACGCCCAGGACa                                            <  2:332943/150‑1 (MQ=255)
                                                                                                        tACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTGGGACGCCCAGGACATc                                          <  1:396626/150‑1 (MQ=255)
                                                                                                          ccGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGaa                                                                        <  2:1579533/118‑1 (MQ=255)
                                                                                                          ccGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGaa                                                                        >  1:1579533/1‑118 (MQ=255)
                                                                                                          ccGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTGGGACGCCCAGGACATCGt                                        <  1:1389825/150‑1 (MQ=255)
                                                                                                          ccGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTGGGACGCCCAGGACATCGt                                        >  1:2191927/1‑150 (MQ=255)
                                                                                                              ttCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCtga                                                                                  <  1:1167179/104‑1 (MQ=255)
                                                                                                              ttCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCtga                                                                                  >  2:1167179/1‑104 (MQ=255)
                                                                                                              ttCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTGGGACGCCCAGGACATCGtcatc                                    >  1:1053212/1‑150 (MQ=255)
                                                                                                              ttCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTGGGACGCCCAGGACATCGtcatc                                    >  1:242868/1‑150 (MQ=255)
                                                                                                               tCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTGGGACGCCCAGGACATCGtcatca                                   <  2:713255/150‑1 (MQ=255)
                                                                                                                         gcgcGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTGGGACGCCCAGGACATCGTCATCAGTATCCACGa                         <  1:1200514/150‑1 (MQ=255)
                                                                                                                             gACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCTCCTGGACTGGGACGCCCAGGACATCGTCATCAGTATCCACGACCAt                     <  1:1065636/150‑1 (MQ=255)
                                                                                                                                ggCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTGGGACGCCCAGGACATCGTCATCAGTATCCACGACCATg                    <  2:2198881/148‑1 (MQ=255)
                                                                                                                                ggCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTGGGACGCCCAGGACATCGTCATCAGTATCCACGACCATg                    >  1:2198881/1‑148 (MQ=255)
                                                                                                                                ggCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTGGGACGCCCAGGACATCGTCATCAGTATCCACGACCATg                    >  1:952980/1‑148 (MQ=255)
                                                                                                                                      gTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTGGGACGCCCAGGACATCGTCATCAGTATCCACGCCCATGGTCCTGgt            <  1:1887667/150‑1 (MQ=255)
                                                                                                                                          cgcgTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTGGGACGCCCAGGACATCGTCATCAGTATCCACGACCATGGTCCTGGTGTAc        <  2:1282679/150‑1 (MQ=255)
                                                                                                                                           gcgTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTg                                                               <  2:228764/94‑1 (MQ=255)
                                                                                                                                           gcgTCTGGCGCCGCCGACCGAGCAGCGCAGAAGCGTACTGGATCGGGTGGTCAAGGCGGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTGGGACGCCCAGGACATCGTCATCAGTATCCACGACCATGGTCCTGGTGTAcc       <  2:1453705/150‑1 (MQ=255)
                                                                                                                                             gTCTGGCGGCGGCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGCTCTGGAAGTTCGCCTGGACTGGGACGCCCAGGACATCGTCATCAGTATCCACGACCATGGTCCTGGTGTACCgc     >  2:395459/1‑150 (MQ=255)
                                                                                                                                             gTCTGGCGGCGGCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTGGGACGCCCAGGACATCGTCATCAGTATCCACGACCATg                    >  1:1824195/1‑135 (MQ=255)
                                                                                                                                             gTCTGGCGGCGGCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTGGGACGCCCAGGACATCGTCATCAGTATCCACGACCATg                    <  2:1824195/135‑1 (MQ=255)
                                                                                                                                             gTCTGGCGGCGGCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGGTCGCCTGGACTGGGACGCCCAGGACAGCGTCATCAGTATCCACGACCATGGTCCTGGTGTACCgc     >  1:695031/1‑150 (MQ=255)
                                                                                                                                              tCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTGGGACGCCCAGGACATCGTCATCAGTATCCACGACCATGGTCCTGGTGTACCgcc    >  2:2240530/1‑150 (MQ=255)
                                                                                                                                                tGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTGGGACGCCCAGGACATCGTCATCAGTATCCACGACCATGGTCCTGGTGTACCgccag  >  2:110839/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCTGCGCGACGGCAAGGTGCCGCGTCTGGCGCCGCCGCCAGACCTGACCCAGGCGTTACTGAATCTGCTGAACAATGCCGCCGATGCCTGCCCTGATGATCTGGAAGTTCGCCTGGACTGGGACGCCCAGGACATCGTCATCAGTATCCACGACCATGGTCCTGGTGTACCGCCAG  >  NC_002947/1030823‑1031116

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