Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,673,595 G→C 7.8% P107P (CCG→CCC PP_4137 → outer membrane siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,673,5950GC 7.8% 313.8 / 14.9 103P107P (CCG→CCCPP_4137outer membrane siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (48/47);  new base C (4/4);  total (52/51)
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.26e-01

Only 100 of 109 total aligned reads displayed.

AGCAACCTGGTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTCCCTGGC  >  NC_002947/4673446‑4673744
                                                                                                                                                     |                                                                                                                                                      
aGCAACCTGGTGGACATGCAGCGCGTGGGGTAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATGAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCg                                                                                                                                                        >  2:1116348/1‑150 (MQ=255)
aGCAACCTGGTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCg                                                                                                                                                        >  2:1591876/1‑150 (MQ=255)
aGCAACCTGGTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCg                                                                                                                                                        <  1:883720/150‑1 (MQ=255)
  cAACCTGGTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTTCGCCgtg                                                                                                                                                      <  2:696605/150‑3 (MQ=255)
    aCCTGGTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTTTCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCggggg                                                                                                                                                    <  1:962209/150‑1 (MQ=255)
        tatggACATGAAGGGCGGGCGGGTGGGCCAAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCa                                                                                                                                                <  2:1852215/148‑1 (MQ=255)
        ggtggACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCa                                                                                                                                                >  2:868529/1‑150 (MQ=255)
         gtggACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACTATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCAc                                                                                                                                               <  1:982710/150‑1 (MQ=255)
         gtggACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCAc                                                                                                                                               >  2:1065427/1‑150 (MQ=255)
          tggACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACg                                                                                                                                              <  2:1395207/150‑1 (MQ=255)
            gACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGTGCCGGGGGCGCACGGt                                                                                                                                            >  2:1744523/1‑150 (MQ=255)
            gACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCAcaga                                                                                                                                            >  1:1852215/1‑147 (MQ=255)
               aTGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTTGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTcc                                                                                                                                         >  1:1884161/1‑150 (MQ=255)
                 gCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCgg                                                                                                                                       <  1:1065427/150‑1 (MQ=255)
                  cAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCACGCAGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGgt                                                                                                                                      >  2:816168/1‑150 (MQ=255)
                  cAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGgt                                                                                                                                      <  2:490730/150‑1 (MQ=255)
                   aGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGgtg                                                                                                                                     >  1:6601/1‑150 (MQ=255)
                   aGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGgtg                                                                                                                                     >  2:641709/1‑150 (MQ=255)
                       cgTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTa                                                                                                                                 <  2:6601/150‑1 (MQ=255)
                          ggggCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTacac                                                                                                                              <  2:54466/150‑1 (MQ=255)
                          ggggCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTacac                                                                                                                              >  2:1310868/1‑150 (MQ=255)
                           gggCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACAcc                                                                                                                             <  1:1766885/150‑1 (MQ=255)
                             gCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCAt                                                                                                                           >  2:1779928/1‑150 (MQ=255)
                             gCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCAt                                                                                                                           >  1:374139/1‑150 (MQ=255)
                             gCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACgg                                                                                                                                             >  2:249865/1‑132 (MQ=255)
                             gCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACgg                                                                                                                                             <  1:249865/132‑1 (MQ=255)
                             gCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCGCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCAt                                                                                                                           >  1:881360/1‑150 (MQ=255)
                                  gCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCg                                                                                                                      >  1:1770820/1‑150 (MQ=255)
                                     gAGTGGCCCGCACCCAGGACGTGCGGGCCTACCAGCCAGGGGTGTTCGCCCATTCGGGATGTAACGATGGCAGCAAGCTGTCGATACGCGGCTCGGGCCTCAATCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACg                                                                                                                   <  2:913632/150‑1 (MQ=255)
                                     gAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACg                                                                                                                   <  2:1676921/150‑1 (MQ=255)
                                     gAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACg                                                                                                                   <  1:1597103/150‑1 (MQ=255)
                                     gAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACg                                                                                                                   <  1:751766/150‑1 (MQ=255)
                                     gAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACGGttga                                                                                                                                         >  1:1371543/1‑126 (MQ=255)
                                     gAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACGGttga                                                                                                                                         <  2:1371543/128‑3 (MQ=255)
                                      aGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACgg                                                                                                                  <  1:445426/150‑1 (MQ=255)
                                          gCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTg                                                                                                              >  2:1298945/1‑150 (MQ=255)
                                            cAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGGGCACGGGTCCGGTGTCTAGACGATGTTCGACGGCCTGcc                                                                                                            >  1:1221152/1‑150 (MQ=255)
                                            cAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGcc                                                                                                            >  1:722404/1‑150 (MQ=255)
                                               cAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGtt                                                                                                         <  1:1298945/150‑1 (MQ=255)
                                                 aCCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGa                                                                                                       <  2:1063386/150‑1 (MQ=255)
                                                  ccAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGAc                                                                                                      <  2:1405687/150‑1 (MQ=255)
                                                  ccAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGAc                                                                                                      >  2:1670491/1‑150 (MQ=255)
                                                  ccAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGAc                                                                                                      <  1:1227654/150‑1 (MQ=255)
                                                          tGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGTCGTTGACCGGCCCgg                                                                                              >  1:489490/1‑150 (MQ=255)
                                                           gCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCggg                                                                                             <  2:722404/150‑1 (MQ=255)
                                                               gCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGggcggc                                                                                         >  1:1437902/1‑150 (MQ=255)
                                                                       gccagGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTa                                                                                >  2:718578/1‑150 (MQ=255)
                                                                          agGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACGGttgag                                                                                                                                        <  2:1735251/92‑4 (MQ=38)
                                                                          agGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACGGttgag                                                                                                                                        >  1:1735251/1‑89 (MQ=38)
                                                                           ggggTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCGCGGTTc                                                                                                                                          >  2:1063521/1‑89 (MQ=255)
                                                                           ggggTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCGCGGTTc                                                                                                                                          <  1:1063521/89‑1 (MQ=255)
                                                                           ggggTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGc                                                                                     <  1:1577023/142‑1 (MQ=255)
                                                                           ggggTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGc                                                                                     >  2:1577023/1‑142 (MQ=255)
                                                                           ggggTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGaa                                                                            >  2:1042259/1‑150 (MQ=255)
                                                                              gtgtTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTg                                                                         <  1:1591876/150‑1 (MQ=255)
                                                                                gtTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGtt                                                                       <  1:1779928/150‑1 (MQ=255)
                                                                                 ttCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCg                                                                                                                                        >  2:300564/1‑85 (MQ=255)
                                                                                 ttCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCg                                                                                                                                        <  1:300564/85‑1 (MQ=255)
                                                                                 ttCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTc                                                                      >  1:771757/1‑150 (MQ=255)
                                                                                    gCCCAATCGGCAGTTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCCACACCATGTCCGACGGCCTGCCGTTGACCCGCCCGGGCGGCACGCCTTACGAACAGTTCGaa                                                                   >  2:319309/1‑150 (MQ=255)
                                                                                      ccAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAAcc                                                                 >  2:570919/1‑150 (MQ=255)
                                                                                       cAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCg                                                                >  2:1594134/1‑150 (MQ=255)
                                                                                       cAATCGGCAGGTAACGATGGCACCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCg                                                                <  1:1213449/150‑1 (MQ=255)
                                                                                        aaTCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGc                                                               >  1:965203/1‑150 (MQ=255)
                                                                                           cGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCtgt                                                            <  1:1596984/150‑1 (MQ=255)
                                                                                           cGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCtgt                                                            <  2:1924972/150‑1 (MQ=255)
                                                                                           cGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCtgt                                                            <  2:1854967/150‑1 (MQ=255)
                                                                                                      gATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGt                                                 >  1:61988/1‑150 (MQ=255)
                                                                                                      gATAGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGCCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGt                                                 >  2:634872/1‑150 (MQ=255)
                                                                                                        tGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGc                                               <  1:1116348/150‑1 (MQ=255)
                                                                                                            aTCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGaa                                           >  2:1838154/1‑150 (MQ=255)
                                                                                                              cAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGt                                         <  2:374139/150‑1 (MQ=255)
                                                                                                               aaGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGtg                                        <  2:1437902/150‑1 (MQ=255)
                                                                                                                   tGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGAc                                                                                                      <  2:1736902/85‑1 (MQ=255)
                                                                                                                   tGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGAc                                                                                                      >  1:1736902/1‑85 (MQ=255)
                                                                                                                      cGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGTAGTGCTGCGTg                                 <  2:783663/150‑1 (MQ=255)
                                                                                                                        aTACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGgc                               >  1:1346396/1‑150 (MQ=255)
                                                                                                                           cgcgGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGcc                            <  1:816168/150‑1 (MQ=255)
                                                                                                                            gcgGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGTTGTCTAGACTATTTTTGCCGGCCTGATGGAGCACGGCCGGAGGTGAGGCCTTGATGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCCAAGTTCTGGGTGGCGCCa                           >  1:1876418/1‑150 (MQ=255)
                                                                                                                            gcgGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCa                           >  1:350405/1‑150 (MQ=255)
                                                                                                                             cgGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCaa                          <  1:854186/150‑1 (MQ=255)
                                                                                                                                cTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATgg                       <  1:1042259/150‑1 (MQ=255)
                                                                                                                                      ccTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTAc                                                                               >  1:1290338/1‑88 (MQ=255)
                                                                                                                                      ccTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTAc                                                                               <  2:1290338/88‑1 (MQ=255)
                                                                                                                                       cTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGAc                >  2:271639/1‑150 (MQ=255)
                                                                                                                                        tCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGAcc               >  2:28873/1‑150 (MQ=255)
                                                                                                                                         cAACCGTGCGCCGGGGGCGCACGGTTCAGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTAc                                                                               >  1:1706767/1‑85 (MQ=255)
                                                                                                                                         cAACCGTGCGCCGGGGGCGCACGGTTCAGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTAc                                                                               <  2:1706767/85‑1 (MQ=255)
                                                                                                                                          aaCCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGCATCGTGCTCGCCGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCaa             >  2:527698/1‑150 (MQ=255)
                                                                                                                                          aaCCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCaa             <  1:1168834/150‑1 (MQ=255)
                                                                                                                                          aaCCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCaa             <  1:1171328/150‑1 (MQ=255)
                                                                                                                                           aCCGTGCGCCCCCGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAg            >  1:1395234/1‑150 (MQ=255)
                                                                                                                                            ccGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGCAGATCCACGTGGCGCCAATGGCTTTGACCAAgg           <  1:527698/150‑1 (MQ=255)
                                                                                                                                              gTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGa                                                      >  2:1188924/1‑105 (MQ=255)
                                                                                                                                              gTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGa                                                      <  1:1188924/105‑1 (MQ=255)
                                                                                                                                              gTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCt         >  2:1108403/1‑150 (MQ=255)
                                                                                                                                                  gcCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTa                                                                                >  1:837770/1‑75 (MQ=255)
                                                                                                                                                  gcCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTa                                                                                <  2:837770/75‑1 (MQ=255)
                                                                                                                                                  gcCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTCCCt     >  2:1866730/1‑150 (MQ=255)
                                                                                                                                                     gggggCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTCCCTGGc  >  1:1872583/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
AGCAACCTGGTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTCCCTGGC  >  NC_002947/4673446‑4673744

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: