Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,252,594 T→A 7.2% intergenic (+81/+162) ureG → / ← PP_2851 urease accessory protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,252,5940TA 7.2% 313.0 / 11.4 108intergenic (+81/+162)ureG/PP_2851urease accessory protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (55/45);  new base A (4/4);  total (59/49)
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.22e-01

Only 100 of 110 total aligned reads displayed.

AGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGC  >  NC_002947/3252452‑3252743
                                                                                                                                              |                                                                                                                                                     
tgtggTTGCAGGCCGTGACCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCcgc                                                                                                                                                <  2:2388700/147‑1 (MQ=255)
aGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCcgc                                                                                                                                                <  1:396493/150‑1 (MQ=255)
   ggTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGaa                                                                                                                                             <  1:1825213/150‑1 (MQ=255)
   ggTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTAAACCCGCGaa                                                                                                                                             <  2:1300440/150‑1 (MQ=255)
   ggTTGCAGGACGTGATCGATTTTATGGTCGAGGGGGTGATGCTGGGGGTGCGAGGCTGACAGATCTTAGATTGTTTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGaa                                                                                                                                             <  1:413717/150‑1 (MQ=37)
    gTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAAt                                                                                                                                            >  2:2152588/1‑150 (MQ=255)
    gTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAAt                                                                                                                                            >  2:1808420/1‑150 (MQ=255)
     ttGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGaaca                                                                                                                                           >  2:1503255/1‑148 (MQ=255)
     ttGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATa                                                                                                                                           <  2:1060111/150‑1 (MQ=255)
        cAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGaacaccg                                                                                                                                        <  1:1503255/150‑6 (MQ=255)
        cAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACgg                                                                                                                                        <  2:729765/150‑1 (MQ=255)
        cAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACgg                                                                                                                                        <  2:1961237/150‑1 (MQ=255)
         aGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGt                                                                                                                                       <  1:1604624/150‑1 (MQ=255)
          ggCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTa                                                                                                                                      >  1:2433342/1‑150 (MQ=255)
           gCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAg                                                                                                                                     >  1:1036955/1‑150 (MQ=255)
           gCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAg                                                                                                                                     >  2:1588848/1‑150 (MQ=255)
                 aTCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTATTCGCGGGTAAACCCGCGaacacc                                                                                                                                         >  2:2449921/1‑136 (MQ=38)
                 aTCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTATTCGCGGGTAAACCCGCGAacacc                                                                                                                                         <  1:2449921/140‑2 (MQ=39)
                  tCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGc                                                                                                                              >  2:424454/1‑150 (MQ=255)
                  tCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTAAACCcgc                                                                                                                                                >  2:2225910/1‑132 (MQ=255)
                  tCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTAAACCcgc                                                                                                                                                <  1:2225910/132‑1 (MQ=255)
                    gATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGcaa                                                                                                                            >  2:1919206/1‑150 (MQ=255)
                    gATTTTACCGGGTAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGcaa                                                                                                                            <  1:1284743/150‑1 (MQ=255)
                     aTTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTATTCGCGGGTAAACCCGCGaa                                                                                                                                             >  2:1789511/1‑132 (MQ=255)
                     aTTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTATTCGCGGGTAAACCCGCGaa                                                                                                                                             <  1:1789511/132‑1 (MQ=255)
                      ttttATCGTGGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGcaaca                                                                                                                          >  2:692643/1‑150 (MQ=255)
                      ttttATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGcaaca                                                                                                                          >  2:2408204/1‑150 (MQ=255)
                       tttATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGGTTACCCGCGAATACGGTAGCGGCCGcaacaa                                                                                                                         >  1:1451551/1‑150 (MQ=255)
                        ttATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAAt                                                                                                                        <  1:1588848/150‑1 (MQ=255)
                          aTCGTCGAGCGCGGGATGCTGGGGGTGCGAGGGTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACAAAACCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGa                                                                                                                      >  1:436912/1‑150 (MQ=255)
                           tcgtcgAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGaa                                                                                                                     <  2:195596/150‑1 (MQ=255)
                                 aGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAAcggcgg                                                                                                               <  1:1808420/150‑1 (MQ=255)
                                   cgcgGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGt                                                                                                             <  1:2414239/150‑1 (MQ=255)
                                        gATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGaaa                                                                                                        <  2:1824881/150‑1 (MQ=255)
                                           gCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGCGATCAACAATGAACGGCGGGTGGAAATTg                                                                                                     <  1:1242105/150‑1 (MQ=255)
                                             tGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTgg                                                                                                           <  2:134435/142‑1 (MQ=255)
                                             tGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGc                                                                                                   <  2:465516/150‑1 (MQ=255)
                                               ggggTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAAtt                                                                                                      >  1:1422054/1‑145 (MQ=255)
                                               ggggTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAAtt                                                                                                      <  2:1422054/145‑1 (MQ=255)
                                                         gCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTAc                                                                                                                                                    >  2:368923/1‑89 (MQ=255)
                                                         gCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTAc                                                                                                                                                    <  1:368923/89‑1 (MQ=255)
                                                           tGACAGATCGAGGACTCTGTGCGGGCATGTCAGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGggccggcc                                                                                     >  1:310342/1‑150 (MQ=255)
                                                            gACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAAcggc                                                                                                                 >  1:1038578/1‑121 (MQ=255)
                                                            gACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAAcggc                                                                                                                 <  2:1038578/121‑1 (MQ=255)
                                                                 aTCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCcagca                                                                                              <  1:6734/135‑1 (MQ=255)
                                                                 aTCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCcagca                                                                                              >  2:6734/1‑135 (MQ=255)
                                                                  tCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTcgc                                                                              >  1:1014373/1‑150 (MQ=255)
                                                                   cGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTcgcg                                                                             <  1:1919206/150‑1 (MQ=255)
                                                                   cGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTcgcg                                                                             <  1:2408204/150‑1 (MQ=255)
                                                                    gAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCgg                                                                            <  2:1036955/150‑1 (MQ=255)
                                                                       gACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGACCACCAGCGCCGGCCCCTTCGCGGGTa                                                                         >  2:135363/1‑150 (MQ=255)
                                                                         ctctGTGCGGACATGCCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTaaa                                                                       >  1:2108721/1‑150 (MQ=255)
                                                                          tctGTGCGGGCATGTCCGCACAGCGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGCaaaa                                                                      >  2:1181623/1‑149 (MQ=255)
                                                                          tctGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAc                                                                      >  2:1537153/1‑150 (MQ=255)
                                                                           ctGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAcc                                                                     >  2:1561877/1‑150 (MQ=255)
                                                                            tgtgCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAccc                                                                    >  1:2102421/1‑150 (MQ=255)
                                                                            tgtgCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGGGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAccc                                                                    >  1:1715405/1‑150 (MQ=255)
                                                                            tgtgCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAAGGCGGTAGCGGCCGCAACAATGGACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTACAccc                                                                    >  1:1697582/1‑150 (MQ=255)
                                                                             gtgCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAGGTGTTATCGGCTGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                   <  2:2433342/150‑1 (MQ=37)
                                                                                   gCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCgg                                                                            <  2:2226178/135‑1 (MQ=255)
                                                                                   gCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCgg                                                                            >  1:2226178/1‑135 (MQ=255)
                                                                                     aTGTCCGCACAGAGACCAACGCGGGGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACg                                                           >  2:2041749/1‑150 (MQ=255)
                                                                                        tCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACgggg                                                        >  2:659065/1‑150 (MQ=255)
                                                                                         ccGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGc                                                       <  1:1561877/150‑1 (MQ=255)
                                                                                         ccGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGc                                                       >  1:1153496/1‑150 (MQ=255)
                                                                                             acaGAGACCAACGCGGTGGCTGGCACTGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGc                                                                  <  1:1866847/135‑1 (MQ=255)
                                                                                             acaGAGACCAACGCGGTGGCTGGCACTGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGc                                                                  >  2:1866847/1‑135 (MQ=255)
                                                                                               agagACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGggccggc                                                                                      >  2:1770711/1‑113 (MQ=255)
                                                                                                   aCCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGc                                                                                                                                    <  2:1496916/63‑1 (MQ=255)
                                                                                                   aCCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGc                                                                                                                                    >  1:1496916/1‑63 (MQ=255)
                                                                                                   aCCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAg                                             >  1:739118/1‑150 (MQ=255)
                                                                                                   aCCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAg                                             <  1:1983570/150‑1 (MQ=255)
                                                                                                    ccAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAgg                                            >  1:1167778/1‑150 (MQ=255)
                                                                                                     cAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAGCAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGc                                           >  2:1787178/1‑150 (MQ=255)
                                                                                                     cAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTcg                                                                               <  1:1135817/114‑1 (MQ=255)
                                                                                                     cAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTcg                                                                               >  2:1135817/1‑114 (MQ=255)
                                                                                                     cAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTcg                                                                               <  1:1135839/114‑1 (MQ=255)
                                                                                                     cAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTcg                                                                               >  2:1135839/1‑114 (MQ=255)
                                                                                                     cAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGc                                           >  2:261304/1‑150 (MQ=255)
                                                                                                      aaCGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGcc                                          >  1:1639705/1‑150 (MQ=255)
                                                                                                           ggTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACggg                                                         >  1:251554/1‑130 (MQ=255)
                                                                                                           ggTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACggg                                                         <  2:251554/130‑1 (MQ=255)
                                                                                                            gTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGaa                                    <  1:424454/150‑1 (MQ=255)
                                                                                                             tggctggcACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAAt                                   <  1:692643/150‑1 (MQ=255)
                                                                                                                 tggcACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTg                               >  2:548336/1‑150 (MQ=255)
                                                                                                                   gcACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAg                             >  2:694408/1‑150 (MQ=255)
                                                                                                                     aCCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCa                           <  1:29227/150‑1 (MQ=255)
                                                                                                                      ccGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGaa                                    <  1:104280/140‑1 (MQ=255)
                                                                                                                      ccGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGaa                                    >  2:104280/1‑140 (MQ=255)
                                                                                                                       cGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAgg                         >  1:52333/1‑150 (MQ=255)
                                                                                                                        ggCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGa                        <  2:349697/150‑1 (MQ=255)
                                                                                                                          cagagCCGGGGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACa                      <  1:1852331/146‑1 (MQ=255)
                                                                                                                           tctGCCGGTGTTCGCGGGTTTACCCGCTAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAg                     >  1:54824/1‑150 (MQ=255)
                                                                                                                           tctGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAg                     >  1:2392277/1‑150 (MQ=255)
                                                                                                                              taccgtattCGCGGGTAAACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTg                  >  1:1987462/8‑150 (MQ=25)
                                                                                                                                    gtTTGCGTGTTTTCCCGCGAATACGGTAGCGGCCGCAACAGTGAACGGCGGGTGGAAACTGGCCAGCAGGGCCGGCCTCCTCGCGGGTACACCCGCTCCTACGGGCCGACCGCGCAGGCCGTCGAATCTTGATCGGGTCAGTTGcccccc            <  2:1113447/150‑1 (MQ=18)
                                                                                                                                     ttCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCca           >  1:852829/1‑150 (MQ=255)
                                                                                                                                       cgcgGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGAACGAACCATACAACTAATGTTGAGCAGGACAGTTGCCCCCcaca         >  2:190673/1‑150 (MQ=21)
                                                                                                                                          gggTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGgt      <  2:1451551/150‑1 (MQ=255)
                                                                                                                                              ttACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGc  <  2:1167778/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                     
AGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACAGGTGTGC  >  NC_002947/3252452‑3252743

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: