Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,678,544 G→C 7.2% 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,678,5440GC 7.2% 875.8 / 44.1 262P107P (CCG→CCCPP_4137outer membrane siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (141/102);  new base C (7/12);  total (148/114)
Fisher's exact test for biased strand distribution p-value = 9.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 283 total aligned reads displayed.

AGCAACCTGGTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTCCCTGG  >  NC_002947/4678395‑4678692
                                                                                                                                                     |                                                                                                                                                    
aGCAACCTGGTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCg                                                                                                                                                      >  2:1080250/1‑150 (MQ=255)
aGCAACCTGGTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCg                                                                                                                                                      >  1:3430972/1‑150 (MQ=255)
  cAACCTGGTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCggg                                                                                                                                                    <  2:5517934/150‑1 (MQ=255)
  cAACCTGGTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCggg                                                                                                                                                    <  2:4593156/150‑1 (MQ=255)
       tggtggACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGgcgc                                                                                                                                               >  2:4845451/1‑150 (MQ=255)
       tggtggACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGgcgc                                                                                                                                               >  1:1256577/1‑150 (MQ=255)
         gtggACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCAc                                                                                                                                             >  1:736751/1‑150 (MQ=255)
          tggACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACg                                                                                                                                            <  2:2170497/150‑1 (MQ=255)
           ggACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACgg                                                                                                                                           <  2:4159139/150‑1 (MQ=255)
           ggACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACgg                                                                                                                                           <  1:1028903/150‑1 (MQ=255)
            gACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACg                                                                                                                                            <  2:5314812/148‑1 (MQ=255)
            gACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACGGt                                                                                                                                          <  2:3088626/150‑1 (MQ=255)
               aTGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTcc                                                                                                                                       <  2:822541/150‑1 (MQ=255)
                ggCAGCGCGGGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGAGGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACGGttgag                                                                                                                                      <  2:5622144/149‑4 (MQ=255)
                 gCAGCGCGTGGGGCAAGGCCGCGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCTGGGGTTTTCGCCCAATCGGCAGGTAAGGTTGGCATCAAGCTGTCGATACGCGGATCGGGCCTAAACCGTGCGCCGGGGGCGCACGGTTCCgg                                                                                                                                     <  1:4157832/150‑1 (MQ=255)
                 gCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCgg                                                                                                                                     <  2:3611520/150‑1 (MQ=255)
                 gCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCgg                                                                                                                                     >  1:2044629/1‑150 (MQ=255)
                 gCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACGGttgagg                                                                                                                                     <  2:3430600/150‑5 (MQ=255)
                   aGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACGGttgaggcc                                                                                                                                   <  2:1773592/150‑7 (MQ=255)
                    gcgcGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCGCGGTTCCGgtgt                                                                                                                                  <  1:4886167/150‑1 (MQ=255)
                    gcgcGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGgtgt                                                                                                                                  >  1:1221674/1‑150 (MQ=255)
                      gcgTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCt                                                                                                                                <  2:2714509/150‑1 (MQ=255)
                        gTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTac                                                                                                                              <  2:3875162/150‑1 (MQ=255)
                        gTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTac                                                                                                                              <  2:1046716/150‑1 (MQ=255)
                        gTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTac                                                                                                                              >  1:887445/1‑150 (MQ=255)
                         tGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTaca                                                                                                                             <  2:4671574/150‑1 (MQ=255)
                          ggggCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTacac                                                                                                                            <  2:4509827/150‑1 (MQ=255)
                             gCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACGGtt                                                                                                                                         <  1:288511/134‑1 (MQ=255)
                                aGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGtt                                                                                                                      <  1:4985589/150‑1 (MQ=255)
                                  gCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGGGTCTACACCATGTTCg                                                                                                                    >  2:5120267/1‑150 (MQ=255)
                                    cGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGAc                                                                                                                  >  1:821046/1‑150 (MQ=255)
                                          gCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTg                                                                                                            >  2:4473129/1‑150 (MQ=255)
                                          gCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTg                                                                                                            >  1:169426/1‑150 (MQ=255)
                                          gCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTg                                                                                                            >  2:4305180/1‑150 (MQ=255)
                                           ccAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGc                                                                                                           >  2:5156961/1‑150 (MQ=255)
                                               cAACCAGGACGTGCTGGCCTACCCGCCTGGGGTGTTCGCCCAATCGGCAGTTAACGATGGCATCATGCTGTCGATACGCGGCTCGGGCCTCAACCTTGCGCCGGGGGCGCACGGTTCCGGTTTCTACACCATGTTCGAGGGCCTGCCGtt                                                                                                       <  2:4883274/150‑1 (MQ=255)
                                               cAACCAGGACGTGCTGGCCTACCAGCCTGGGGTGTTCGCCAAATCGGCTGGTAACTATGGCATCACGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGtt                                                                                                       <  1:3560980/150‑1 (MQ=255)
                                                   cAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGAcc                                                                                                   >  1:4516301/1‑150 (MQ=255)
                                                            cTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGggc                                                                                          >  2:3986707/1‑150 (MQ=255)
                                                               gCCTACCAGCCTGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCGGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGggcggc                                                                                       <  1:4280785/150‑1 (MQ=255)
                                                                 cTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCAc                                                                                     >  1:4847012/1‑150 (MQ=255)
                                                                  tACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACg                                                                                    >  1:3339508/1‑150 (MQ=255)
                                                                        ccagGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTAc                                                                              >  1:1014777/1‑150 (MQ=255)
                                                                           ggggTTTTCGCCAAATCGGCATTTAACTATGGAATCAAGCTGTCTATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGAACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGaa                                                                           <  1:2579205/150‑1 (MQ=255)
                                                                           ggggTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGaa                                                                           >  2:2596158/1‑150 (MQ=255)
                                                                           ggggTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGaa                                                                           <  2:2212859/150‑1 (MQ=255)
                                                                           ggggTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGaa                                                                           >  1:3518076/1‑150 (MQ=255)
                                                                            gggTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAAc                                                                          >  1:69612/1‑150 (MQ=255)
                                                                            gggTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAAc                                                                          >  2:2822673/1‑150 (MQ=255)
                                                                             ggTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACt                                                                         <  1:1344246/150‑1 (MQ=255)
                                                                              gtgtTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTg                                                                        >  2:3928864/1‑150 (MQ=255)
                                                                                gtTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGtt                                                                      <  1:3327942/150‑1 (MQ=255)
                                                                                 ttCGCCCAATCGGCAGGTAACGATAGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTc                                                                     <  2:1206379/150‑1 (MQ=255)
                                                                                    gCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGaa                                                                  >  1:2457612/1‑150 (MQ=255)
                                                                                     cccAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCggg                                                                                                                                                    <  1:4360665/67‑1 (MQ=255)
                                                                                     cccAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAAc                                                                 <  1:2449579/150‑1 (MQ=255)
                                                                                     cccAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAAc                                                                 >  2:3779753/1‑150 (MQ=255)
                                                                                     cccAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAAc                                                                 <  2:1474109/150‑1 (MQ=255)
                                                                                     cccAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAAc                                                                 <  2:3012587/150‑1 (MQ=255)
                                                                                        aaTCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGc                                                              >  2:43974/1‑150 (MQ=255)
                                                                                           cGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCtgt                                                           <  2:3740897/150‑1 (MQ=255)
                                                                                               aGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCt                                                       >  1:5117433/1‑150 (MQ=255)
                                                                                                    aCGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGcc                                                  <  1:4627406/150‑1 (MQ=255)
                                                                                                      gATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGt                                                >  1:2952782/1‑150 (MQ=255)
                                                                                                      gATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGt                                                >  2:3252096/1‑150 (MQ=255)
                                                                                                        tGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGc                                              <  2:959687/150‑1 (MQ=255)
                                                                                                         ggCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGGGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGcc                                             <  1:3741663/150‑1 (MQ=255)
                                                                                                           cATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGa                                           >  1:233118/1‑150 (MQ=255)
                                                                                                            aTCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGaa                                          >  1:2344896/1‑150 (MQ=255)
                                                                                                             tCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAg                                         <  2:407752/150‑1 (MQ=255)
                                                                                                             tCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAg                                         >  2:1533212/1‑150 (MQ=255)
                                                                                                                aGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGtgc                                      >  1:1819992/1‑150 (MQ=255)
                                                                                                                aGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGGTGCCCGGCCCGGGCGGCACGCCGTACGAACTGGTCGAACCGCTGTGGCTGAGCCGTGCCGAAGtgc                                      >  1:1934872/1‑150 (MQ=255)
                                                                                                                  cTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGtgctg                                    >  1:109429/1‑150 (MQ=255)
                                                                                                                  cTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGtgctg                                    >  1:5136455/1‑150 (MQ=255)
                                                                                                                   tGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGtgctgc                                   <  2:2924485/150‑1 (MQ=255)
                                                                                                                   tGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGtgctgc                                   <  2:2781552/150‑1 (MQ=255)
                                                                                                                   tGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGtgctgc                                   >  2:1846269/1‑150 (MQ=255)
                                                                                                                     tCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGt                                 >  2:633272/1‑150 (MQ=255)
                                                                                                                     tCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGt                                 <  2:4344630/150‑1 (MQ=255)
                                                                                                                       gATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTgg                               >  1:2672572/1‑150 (MQ=255)
                                                                                                                       gATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTgg                               >  2:1336828/1‑150 (MQ=255)
                                                                                                                       gATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTgg                               >  1:5178353/1‑150 (MQ=255)
                                                                                                                            gcgGCTCGGGCCTCATCCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCa                          >  2:3500324/1‑150 (MQ=255)
                                                                                                                            gcgGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCa                          >  1:4052675/1‑150 (MQ=255)
                                                                                                                            gcgGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCa                          >  1:1444729/1‑150 (MQ=255)
                                                                                                                            gcgGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCATGCCGAAGTGCTGCGTGGCGCCa                          >  1:5563279/1‑150 (MQ=255)
                                                                                                                                cTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATgg                      <  1:2822673/150‑1 (MQ=255)
                                                                                                                                    ggCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCttt                  <  2:2272818/150‑1 (MQ=255)
                                                                                                                                    ggCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCttt                  <  2:3543031/150‑1 (MQ=255)
                                                                                                                                        tCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGAcc              >  1:3647758/1‑150 (MQ=255)
                                                                                                                                        tCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGAcc              >  2:238525/1‑150 (MQ=255)
                                                                                                                                          aaCCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCaa            >  2:770452/1‑150 (MQ=255)
                                                                                                                                           aCCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAg           >  2:1352008/1‑150 (MQ=255)
                                                                                                                                           aCCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAg           >  1:3073245/1‑150 (MQ=255)
                                                                                                                                            ccGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAgg          >  1:1581924/1‑150 (MQ=255)
                                                                                                                                            ccGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAgg          >  1:3889587/1‑150 (MQ=255)
                                                                                                                                            ccGTGCGCCCCCGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAgg          >  1:4000412/1‑150 (MQ=255)
                                                                                                                                                  gcCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTCCCt    >  1:1414870/1‑150 (MQ=255)
                                                                                                                                                   ccGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCATGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTCCCTg   <  2:5563279/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGCAACCTGGTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTCCCTGG  >  NC_002947/4678395‑4678692

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