Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,106,716 G→T 5.1% A148A (GCC→GCA PP_1880 ← outer membrane autotransporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,106,7160GT 5.1% 325.2 / 4.8 99A148A (GCC→GCAPP_1880outer membrane autotransporter
Reads supporting (aligned to +/- strand):  ref base G (40/53);  new base T (2/3);  total (43/56)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

GAGCCGCGGCCGTCGTCCTCAAGAATGCGCGCGCCATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGCCACGGGCGTCACTGCGGATGGTCGAGCGGTGTACGGTGGC  >  NC_002947/2106570‑2106862
                                                                                                                                                  |                                                                                                                                                   
gagCCGCGGCCGTCGTCCTCAAGAATGCGCGCGCCATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCa                                                                                                                                                  <  1:1032074/150‑1 (MQ=255)
gagCCGCGGCCGTCGTCCTCAAGAATGCGCGCGCCATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGAAGTTCGGCa                                                                                                                                                  <  1:954087/150‑1 (MQ=255)
       ggCCGTCGTCCTCAAGAATGCGCGCGCCATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTgg                                                                                                                                           <  1:1129048/150‑1 (MQ=255)
       ggCCGTCGTCCTCAAGAATGCGCGCGCCATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTgg                                                                                                                                           <  1:1509920/150‑1 (MQ=255)
          cgtcgtCCTCAAGAATGCGCGCGCCATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCTGCCGAACTgcagc                                                                                                                                        <  2:2136962/150‑5 (MQ=255)
           gtcgtcCTCAAGAATGCGCGCGCCATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTggacaa                                                                                                                                       <  1:84197/150‑5 (MQ=255)
           gtcgtcCTCAAGAATGCGCGCGCCATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCGCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCAtt                                                                                                                                       <  1:536056/150‑1 (MQ=255)
             cgtcCTCAAGAATGCGCGCGCCATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTAGCTGCTCTGCAGTTCGTCAGATCTGGCATTGa                                                                                                                                     <  2:529428/150‑1 (MQ=255)
               tcCTCAAGAATGCGCGCGCTATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAgg                                                                                                                                   <  2:1313421/150‑1 (MQ=255)
               tcCTCAAGAATGCGCGCGCCATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAgg                                                                                                                                   >  1:1880190/1‑150 (MQ=255)
               tcCTCAAGAATGCGCGCGCCATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAgg                                                                                                                                   <  2:1313984/150‑1 (MQ=255)
                  tCAAGAATGCGCGCGCCATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTgg                                                                                                                                           <  2:2090261/139‑1 (MQ=255)
                  tCAAGAATGCGCGCGCCATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGATGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTgg                                                                                                                                           >  1:2090261/1‑139 (MQ=255)
                    aaGAATGCGCGCGCCATAATTGCCGCCGCTGATCCAGCTGCCGTCCCTCGCAATCAGGGTGCCGGCAACCAGGCGTGCGCTCCCGCTGCTGGAACCGCGCCGGCGTAGCTCGCTGCCCCGCCGTTCAGCAGAACGGGCATTAAGGGTcgc                                                                                                                              >  2:271505/1‑150 (MQ=255)
                        aTGCGCGCGCTATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCt                                                                                                                          <  1:105888/150‑1 (MQ=255)
                        aTGCGCGCGCCATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCt                                                                                                                          <  1:15445/150‑1 (MQ=255)
                        aTGCGCGCGCCATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCt                                                                                                                          <  1:2354974/150‑1 (MQ=255)
                         tGCGCGCGCCATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCtt                                                                                                                         <  1:672487/150‑1 (MQ=255)
                         tGCGCGCGCCATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCtt                                                                                                                         <  2:917415/150‑1 (MQ=255)
                           cgcgcgcCATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGt                                                                                                                       >  2:836753/1‑150 (MQ=255)
                                gcCATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTAcgcc                                                                                                                  >  1:2096430/1‑150 (MQ=255)
                                  cATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTAcgccgc                                                                                                                >  2:628449/1‑150 (MQ=255)
                                   aTAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGTGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCgcgc                                                                                                                            >  2:215269/1‑137 (MQ=255)
                                   aTAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGTGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCgcgc                                                                                                                            <  1:215269/137‑1 (MQ=255)
                                    tAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCt                                                                                                             <  1:1647982/150‑1 (MQ=255)
                                     aaTTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCtt                                                                                                            <  1:2006577/150‑1 (MQ=255)
                                     aaTTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCtt                                                                                                            <  2:507369/150‑1 (MQ=255)
                                        tGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGa                                                                                                         <  1:1266793/150‑1 (MQ=255)
                                            gccgcTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCTGCCGAACTgca                                                                                                                                          >  1:1657394/1‑112 (MQ=255)
                                            gccgcTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCTGCCGAACTgca                                                                                                                                          <  2:1657394/114‑3 (MQ=255)
                                               gcTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCt                                                                                                  <  2:1278335/150‑1 (MQ=255)
                                                  gATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTc                                                                                               >  1:216761/1‑150 (MQ=255)
                                                    tGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGGGGCATCGCTGCTTGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGt                                                                                             <  2:1269864/150‑1 (MQ=255)
                                                       tgctgcCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGa                                                                                          >  2:199114/1‑150 (MQ=255)
                                                       tgctgcCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGa                                                                                          <  2:2096430/150‑1 (MQ=255)
                                                            ccGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAg                                                                                     >  2:257803/1‑150 (MQ=255)
                                                             cGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAgg                                                                                    <  2:2376006/150‑1 (MQ=255)
                                                              gTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGgc                                                                                   <  1:628449/150‑1 (MQ=255)
                                                                  cTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCTGCCGAACTgca                                                                                                                                          >  2:559883/1‑90 (MQ=255)
                                                                  cTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCTGCCGAACTgca                                                                                                                                          <  1:559883/92‑3 (MQ=255)
                                                                   tGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGcgcgcg                                                                              >  2:1116747/1‑150 (MQ=255)
                                                                    ggCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGa                                                                             >  2:1276795/1‑150 (MQ=255)
                                                                       aaTCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGcc                                                                          >  1:1222037/1‑150 (MQ=255)
                                                                        aTCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCg                                                                         >  1:406340/1‑150 (MQ=255)
                                                                         tCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGt                                                                        >  2:1711621/1‑150 (MQ=255)
                                                                           aGGGTGCCGGCAACCAGGCGTGGGAAATCGCTTCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGttt                                                                      <  1:611283/150‑1 (MQ=255)
                                                                            gggTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTc                                                                     >  1:363706/1‑150 (MQ=255)
                                                                            gggTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTc                                                                     >  2:1206451/1‑150 (MQ=255)
                                                                               tGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCg                                                                  >  1:1192424/1‑150 (MQ=255)
                                                                                gCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGc                                                                 <  2:557946/150‑1 (MQ=255)
                                                                                 ccGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGcc                                                                >  1:1116135/1‑150 (MQ=255)
                                                                                    gCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGa                                                                                                                                     >  1:504337/1‑79 (MQ=255)
                                                                                    gCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGa                                                                                                                                     <  2:504337/79‑1 (MQ=255)
                                                                                    gCAACCAGGCCTGCGCCCTCGCTGCTGGCCCCCGTCCCGCGTAGTTCGCTGCGCTGAAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGcccat                                                             <  1:1552684/150‑3 (MQ=255)
                                                                                     cAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTg                                                            <  2:405706/150‑1 (MQ=255)
                                                                                        ccAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCGTTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTgcag                                                         >  1:576000/1‑150 (MQ=255)
                                                                                        ccAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGa                                                                                                                                     <  1:1039263/75‑1 (MQ=255)
                                                                                        ccAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGa                                                                                                                                     >  2:1039263/1‑75 (MQ=255)
                                                                                        ccAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTgcag                                                         <  2:1697682/150‑1 (MQ=255)
                                                                                         cAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTgcagc                                                        >  2:1594428/1‑150 (MQ=255)
                                                                                          aGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACCAGCACATCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTgcagca                                                       >  1:1834908/1‑150 (MQ=255)
                                                                                            gCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACt                                                     >  1:754818/1‑150 (MQ=255)
                                                                                             cGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTc                                                    <  1:95928/150‑1 (MQ=255)
                                                                                                     cTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGa                                                                                                                                     <  2:1492092/62‑1 (MQ=255)
                                                                                                     cTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGa                                                                                                                                     >  1:1492092/1‑62 (MQ=255)
                                                                                                     cTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAg                                            >  2:1392766/1‑150 (MQ=255)
                                                                                                          tgctgGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGCCAc                                       <  1:1594428/150‑1 (MQ=255)
                                                                                                             tgGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGCCACggg                                    >  2:1348924/1‑150 (MQ=255)
                                                                                                              ggCGCACGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGCCACGGGc                                   <  1:1273084/150‑1 (MQ=255)
                                                                                                               gcgcCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCAc                                                                                                     <  2:1491859/83‑1 (MQ=255)
                                                                                                               gcgcCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCAc                                                                                                     >  1:1491859/1‑83 (MQ=255)
                                                                                                                 gcCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCa                                             >  1:1849055/1‑137 (MQ=255)
                                                                                                                 gcCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCa                                             <  2:1849055/137‑1 (MQ=255)
                                                                                                                 gcCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGCCACgg                                     <  1:743952/145‑1 (MQ=255)
                                                                                                                 gcCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGCCACgg                                     >  2:743952/1‑145 (MQ=255)
                                                                                                                 gcCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGCCACGGGCGTc                                <  2:216761/150‑1 (MQ=255)
                                                                                                                   ccGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGCCACGGGCGTCAc                              >  1:90612/1‑150 (MQ=255)
                                                                                                                     gTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGCCACGGGCgccgcgg                            <  1:1116747/150‑7 (MQ=255)
                                                                                                                        ccGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTc                                                                                                           >  2:339752/1‑68 (MQ=255)
                                                                                                                        ccGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTc                                                                                                           <  1:339752/68‑1 (MQ=255)
                                                                                                                         cgcgTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGCCACGGGCGTCACTGCGGa                        >  2:2261436/1‑150 (MQ=255)
                                                                                                                         cgcgTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGCCACGGGCGTCACTGCGGa                        <  1:625139/150‑1 (MQ=255)
                                                                                                                          gcgTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGCCACGGGCGTCACTGCGGAt                       <  1:1206451/150‑1 (MQ=255)
                                                                                                                          gcgTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGCCACGGGCGTCACTGCGGAt                       >  1:575890/1‑150 (MQ=255)
                                                                                                                              aGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGCCACGGGCGTCACTGCGGATGGTc                   <  2:575890/150‑1 (MQ=255)
                                                                                                                              aGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGCCACGGGCGTCACTGCGGATGGTc                   >  2:1570126/1‑150 (MQ=255)
                                                                                                                                ttCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGCCACGGGCGTCACTGCGGATGGTCGa                 <  1:199114/150‑1 (MQ=255)
                                                                                                                                 tCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGTCACGGGCGTCACTGCGGATGGTCGAg                <  2:806472/150‑1 (MQ=255)
                                                                                                                                   gctgcgctgcAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCTGCCCCTGCAGCACTCAGCGCCAGGCCACGGGCGTCACTGCGGATGGTCGAGCg              >  1:1565092/1‑150 (MQ=255)
                                                                                                                                   gctgcgctgcAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGCCACGGGCGTCACTGCGGATGGTCGAGCg              >  1:619330/1‑150 (MQ=255)
                                                                                                                                   gctgcgctgcAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGCCACGGGCGTCACTGCGGATGGTCGAGCg              >  2:1246726/1‑150 (MQ=255)
                                                                                                                                     tgcgctgcAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCTGCCCCTGCAGCACTCAGCGCCAGGCCACGGGCGTCACTGCGGATGGTCGAGCGgt            <  2:1565092/150‑1 (MQ=255)
                                                                                                                                     tgcgctgcAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAg                                            >  1:1537130/1‑118 (MQ=255)
                                                                                                                                     tgcgctgcAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAg                                            <  2:1537130/118‑1 (MQ=255)
                                                                                                                                     tgcgctgcAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGCCACGGGCGTCACTGCGGATGGTCGAGCGgt            <  2:1140573/150‑1 (MQ=255)
                                                                                                                                            cAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGCCACGGGCGTCACTGCGGATGGTCGAGCGGTGTACGGt     <  2:1880190/150‑1 (MQ=255)
                                                                                                                                               ttCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGt                                                                                             <  1:2341239/59‑1 (MQ=255)
                                                                                                                                               ttCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGt                                                                                             >  2:2341239/1‑59 (MQ=255)
                                                                                                                                               ttCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGCCACGGGCGTCACTGCGGATGGTCGAGCGGTGTACGGtggc  <  1:2261436/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GAGCCGCGGCCGTCGTCCTCAAGAATGCGCGCGCCATAATTGCCGCCGCTGATGGTGCTGCCGTCCCTGGCAATCAGGGTGCCGGCAACCAGGCGTGCGCCCTCGCTGCTGGCGCCCGTCCCGCGTAGTTCGCTGCGCTGCAGTTCGGCAGAACTGGCATTGAGGGTCGCGCCTTGTACGCCGCCTTCGATCACGCTGTCGTTGATCAAGGCGCGCGAGCCGTTTCCCGCCCCTGCAGCACTCAGCGCCAGGCCACGGGCGTCACTGCGGATGGTCGAGCGGTGTACGGTGGC  >  NC_002947/2106570‑2106862

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: