Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,673,595 G→C 7.0% 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.0% 447.1 / 21.4 142P107P (CCG→CCCPP_4137outer membrane siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (78/54);  new base C (2/8);  total (80/62)
Fisher's exact test for biased strand distribution p-value = 2.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.13e-01

Only 100 of 149 total aligned reads displayed.

GCAACCTGGTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTCCCTG  >  NC_002947/4673447‑4673742
                                                                                                                                                    |                                                                                                                                                   
gCAACCTGGTGGACATGGAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGTTCGGGCCTCAACCGTGCGCCgg                                                                                                                                                    <  2:733105/150‑1 (MQ=255)
  aaCCTGGTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCgggg                                                                                                                                                  <  1:2569540/150‑1 (MQ=255)
   aCCTGGTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCgtcag                                                                                                                                                 <  2:699450/150‑5 (MQ=255)
          ggACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACgg                                                                                                                                          <  1:1245099/150‑1 (MQ=255)
            aCATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACGGtt                                                                                                                                        <  1:2659990/150‑1 (MQ=255)
              aTGCAGCGCGTGGGGCAAGGCCGAGTGGTCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATAGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACGGttga                                                                                                                                      <  2:319392/150‑3 (MQ=255)
              aTGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTcc                                                                                                                                      >  2:2924717/1‑150 (MQ=255)
               tGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGGTGGCCTCAAAGCTAGGGGGGTCGACTTAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCg                                                                                                                                     <  1:1898175/150‑1 (MQ=255)
                  aGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTACCGgtg                                                                                                                                  >  2:553033/1‑150 (MQ=255)
                      cgTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTa                                                                                                                              <  2:2380018/150‑1 (MQ=255)
                          gggCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACTGTGCGCCCCCGGCGCACGGtt                                                                                                                                        >  2:1354651/1‑136 (MQ=255)
                            gCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCAt                                                                                                                        >  1:1984623/1‑150 (MQ=255)
                            gCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGATGTc                                                                                                                                >  1:877784/1‑142 (MQ=255)
                              aaGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGt                                                                                                                      <  1:833089/150‑1 (MQ=255)
                               aGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGtt                                                                                                                     <  2:2822008/150‑1 (MQ=255)
                                 gCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCgg                                                                                                                                                    <  1:638659/117‑1 (MQ=255)
                                   cGAGTTGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGAc                                                                                                                 <  2:1509111/150‑1 (MQ=255)
                                    gAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTc                                                                                                                                <  2:2502495/134‑1 (MQ=255)
                                    gAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTc                                                                                                                                >  1:2502495/1‑134 (MQ=255)
                                         gCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTg                                                                                                           >  1:1039817/1‑150 (MQ=255)
                                          ccAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGc                                                                                                          >  2:655735/1‑150 (MQ=255)
                                            aGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCg                                                                                                        >  1:1651899/1‑150 (MQ=255)
                                               aaCCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTg                                                                                                     >  1:1771504/1‑150 (MQ=255)
                                               aaCCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACGGt                                                                                                                                         <  2:976965/114‑1 (MQ=255)
                                                  cAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGgtgt                                                                                                                                 <  1:774432/119‑1 (MQ=255)
                                                  cAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCCCCGGCGCACg                                                                                                                                           <  1:2856460/109‑1 (MQ=255)
                                                  cAGGACGTGCTGGCCTACCAGCCAGGGGTGTACGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGAcc                                                                                                  >  2:1890558/1‑150 (MQ=255)
                                                         tGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCgg                                                                                           >  1:596450/1‑150 (MQ=255)
                                                          gCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCggg                                                                                          <  1:2387570/150‑1 (MQ=255)
                                                            tGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGACTTGACCGGGCGCCGGTGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGggcg                                                                                        >  1:370915/1‑150 (MQ=255)
                                                                cTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGa                                                                                                                  >  2:1286481/1‑120 (MQ=255)
                                                                cTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGa                                                                                                                  <  1:1286481/120‑1 (MQ=255)
                                                                cTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGcc                                                                                                             >  2:1470609/1‑125 (MQ=255)
                                                                cTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGcc                                                                                                             <  1:1470609/125‑1 (MQ=255)
                                                                cTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGggcgg                                                                                       <  2:490967/147‑1 (MQ=255)
                                                                  aCCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCTGCCCGGGCGGCACGc                                                                                  >  1:722816/1‑150 (MQ=255)
                                                                   ccagccagGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGgt                                                                                                                                   >  2:2921369/1‑100 (MQ=255)
                                                                   ccagccagGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGcc                                                                                 >  1:2906404/1‑150 (MQ=255)
                                                                         agGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGa                                                                           >  1:2225746/1‑150 (MQ=255)
                                                                             gtgtTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACt                                                                        >  2:2354991/1‑149 (MQ=255)
                                                                             gtgtTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACt                                                                        <  1:2354991/149‑1 (MQ=255)
                                                                              tgtTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGt                                                                      <  2:2225746/150‑1 (MQ=255)
                                                                               gtTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGtt                                                                     <  2:2980649/150‑1 (MQ=255)
                                                                                     ccAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGAc                                                                                                   <  1:740299/114‑1 (MQ=255)
                                                                                     ccAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTAc                                                                             >  2:552794/1‑136 (MQ=255)
                                                                                     ccAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTAc                                                                             <  1:552794/136‑1 (MQ=255)
                                                                                      cAATCGGCAGGTAACGATGGCATGAAGCTGACGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCACACGGTTCCGGTGTCTACACCATCTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTAAGCACTGTTCGAAccc                                                              >  2:301018/1‑149 (MQ=255)
                                                                                       aaTCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGc                                                             <  1:655735/150‑1 (MQ=255)
                                                                                            gCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTAc                                                                             <  1:1024438/129‑1 (MQ=255)
                                                                                            gCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTAc                                                                             >  2:1024438/1‑129 (MQ=255)
                                                                                            gCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTgg                                                        >  1:1830352/1‑150 (MQ=255)
                                                                                             cAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGTCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGGCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGc                                                       <  1:2924717/150‑1 (MQ=255)
                                                                                              aGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCt                                                      >  1:259545/1‑150 (MQ=255)
                                                                                              aGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCt                                                      >  2:1938390/1‑150 (MQ=255)
                                                                                              aGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCt                                                      >  1:1749544/1‑150 (MQ=255)
                                                                                                gTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGa                                                    <  2:596450/150‑1 (MQ=255)
                                                                                                   aCGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGcc                                                 >  2:165830/1‑150 (MQ=255)
                                                                                                     gATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGCGCCGt                                               >  2:1432242/1‑150 (MQ=255)
                                                                                                     gATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGt                                               >  2:1131700/1‑150 (MQ=255)
                                                                                                     gATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGCTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGt                                               <  2:1984623/150‑1 (MQ=255)
                                                                                                       tGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATg                                                                                                                       >  2:1790722/1‑76 (MQ=255)
                                                                                                       tGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGc                                             <  2:2220322/150‑1 (MQ=255)
                                                                                                        ggCAACAAGCTGTCGATCGGCGGCGCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGcc                                            <  2:2475723/150‑1 (MQ=255)
                                                                                                          cATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTc                                                                    <  2:2815808/124‑1 (MQ=255)
                                                                                                          cATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTc                                                                    >  1:2815808/1‑124 (MQ=255)
                                                                                                           aTCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGaa                                         >  1:2444730/1‑150 (MQ=255)
                                                                                                            tCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTGCGCACCGCTGTGCCTTAGCCCTGCCGAAg                                        >  2:592109/1‑150 (MQ=255)
                                                                                                              aaGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGAc                                                                                                                 >  2:1358025/1‑75 (MQ=255)
                                                                                                              aaGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGAc                                                                                                                 <  1:1358025/75‑1 (MQ=255)
                                                                                                              aaGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGtg                                      >  1:505154/1‑150 (MQ=255)
                                                                                                                   gTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGcc                                                 >  2:833134/1‑134 (MQ=255)
                                                                                                                      gATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTacac                                                                                                                           <  2:460353/57‑1 (MQ=255)
                                                                                                                      gATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTacac                                                                                                                           >  1:460353/1‑57 (MQ=255)
                                                                                                                      gATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTTGCTGAGCCGTGCCGAAGTGCTGCGTgg                              >  2:1864505/1‑150 (MQ=255)
                                                                                                                      gATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTgg                              >  1:2962269/1‑150 (MQ=255)
                                                                                                                       aTACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGCACTGTTCGAACCGCTGTCGCTGAGCCGTGCCGAAGTGCTGCGTGgc                             <  1:2008305/150‑1 (MQ=255)
                                                                                                                       aTACGCGGCTCGGGCCTCAACCGTGCGCCGGGAGACCACCTTTCCGGTGTCTAGACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCCAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGgc                             <  2:1780175/150‑1 (MQ=255)
                                                                                                                           gcgGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCa                         >  2:324936/1‑150 (MQ=255)
                                                                                                                             ggCTCGGGCCTCAACCGTGCGCCGGGGGCGCACTGTTCCGGTGTCTACACCATGTTCGACGGCCTTCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAAt                       <  2:2045860/150‑1 (MQ=255)
                                                                                                                             ggCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCtgt                                                          >  2:953966/1‑115 (MQ=255)
                                                                                                                             ggCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCtgt                                                          <  1:953966/115‑1 (MQ=255)
                                                                                                                             ggCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTATAGAAAGAATTATTCTCCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAAt                       >  2:755500/1‑150 (MQ=255)
                                                                                                                                tCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTTCGTGGCGCCAATGGc                    >  2:711213/1‑150 (MQ=255)
                                                                                                                                 cGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCt                   <  2:744605/150‑1 (MQ=255)
                                                                                                                                   ggCCTCATCCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGGACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCACTGGCttt                 <  2:72771/150‑1 (MQ=255)
                                                                                                                                    accggaaCCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGa                                                    >  1:1318648/6‑114 (MQ=255)
                                                                                                                                      ctaaaCCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTTTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGAc              <  2:850784/147‑1 (MQ=255)
                                                                                                                                      cTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGa                                                    <  1:2272442/112‑1 (MQ=255)
                                                                                                                                      cTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGa                                                    >  2:2272442/1‑112 (MQ=255)
                                                                                                                                         aaCCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGc                                             >  2:1261412/1‑116 (MQ=255)
                                                                                                                                         aaCCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGc                                             <  1:1261412/116‑1 (MQ=255)
                                                                                                                                         aaCCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCaa           <  1:165830/150‑1 (MQ=255)
                                                                                                                                         aaCCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCaa           >  2:67914/1‑150 (MQ=255)
                                                                                                                                          aCCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAg          <  2:1225772/150‑1 (MQ=255)
                                                                                                                                              tGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACt                                                                        <  1:21863/84‑1 (MQ=255)
                                                                                                                                              tGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACt                                                                        >  2:21863/1‑84 (MQ=255)
                                                                                                                                               gcgcCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTTCCGAATTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTcc     <  1:592109/150‑1 (MQ=255)
                                                                                                                                                  ccGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTCCCTg  >  1:1155572/1‑150 (MQ=255)
                                                                                                                                                   cGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGa                                                    >  1:2116705/1‑99 (MQ=255)
                                                                                                                                                   cGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGa                                                    <  2:2116705/99‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCAACCTGGTGGACATGCAGCGCGTGGGGCAAGGCCGAGTGGCCAGCAACCAGGACGTGCTGGCCTACCAGCCAGGGGTGTTCGCCCAATCGGCAGGTAACGATGGCATCAAGCTGTCGATACGCGGCTCGGGCCTCAACCGTGCGCCGGGGGCGCACGGTTCCGGTGTCTACACCATGTTCGACGGCCTGCCGTTGACCGGCCCGGGCGGCACGCCTTACGAACTGTTCGAACCGCTGTGGCTGAGCCGTGCCGAAGTGCTGCGTGGCGCCAATGGCTTTGACCAAGGCTCCCTG  >  NC_002947/4673447‑4673742

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