Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,609,045 G→C 5.9% R202G (CGC→GGC)  waaA ← 3‑deoxy‑D‑manno‑octulosonic acid transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,609,0450GC 5.9% 366.2 / 12.3 118R202G (CGC→GGC) waaA3‑deoxy‑D‑manno‑octulosonic acid transferase
Reads supporting (aligned to +/- strand):  ref base G (52/59);  new base C (4/3);  total (56/62)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-01

Only 100 of 121 total aligned reads displayed.

GATCAAGGCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTGGCCA  >  NC_002947/5608896‑5609194
                                                                                                                                                       |                                                                                                                                                     
gATCAAGGCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCAcgcg                                                                                                                                                       >  2:488951/1‑150 (MQ=255)
    aaGGCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCa                                                                                                                                                   >  2:2213085/1‑150 (MQ=255)
     aGGCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAg                                                                                                                                                  >  2:893157/1‑150 (MQ=255)
      ggCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTTTCGTAGCCGGTGACTTGCACGCGCTCAgg                                                                                                                                                 <  1:64076/150‑1 (MQ=255)
       gCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAggg                                                                                                                                                >  2:409253/1‑150 (MQ=255)
       gCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTACCCGCGCttgggg                                                                                                                                                >  2:662666/1‑145 (MQ=255)
        cATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGc                                                                                                                                               <  2:777358/150‑1 (MQ=255)
          tCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGt                                                                                                                                             <  2:946933/150‑1 (MQ=255)
          tCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGt                                                                                                                                             <  1:2360547/150‑1 (MQ=255)
           cTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTg                                                                                                                                            >  1:1826626/1‑150 (MQ=255)
                ccGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCCCTGAGCGCGTGCaagt                                                                                                                                       <  2:1385651/150‑4 (MQ=255)
                  gtcgtGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAg                                                                                                                                     <  1:1236981/150‑1 (MQ=255)
                  gtcgtGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAg                                                                                                                                     <  2:424223/150‑1 (MQ=255)
                       gggTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCAc                                                                                                                                <  2:1424144/150‑1 (MQ=255)
                        ggTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACg                                                                                                                               <  2:139660/150‑1 (MQ=255)
                             tGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAgc                                                                                                                          >  1:833049/1‑150 (MQ=255)
                             tGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAgc                                                                                                                          >  1:801334/1‑150 (MQ=255)
                               gcggcAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAgcgc                                                                                                                        >  2:2031485/1‑150 (MQ=255)
                               gcggcAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTAAATGGTCAGGTAAAACCCGATCGAGCCTGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAgcgc                                                                                                                        <  1:1488110/150‑1 (MQ=255)
                                 ggcAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTg                                                                                                                      >  2:910025/1‑150 (MQ=255)
                                     aTCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGcc                                                                                                                  >  1:1063428/1‑150 (MQ=255)
                                     aTCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGGGGCAGCAGCTGTTCGTCGGTGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCCAGGAGACGGTGGGcc                                                                                                                  >  1:1993889/1‑150 (MQ=255)
                                        cacacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGTTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCt                                                                                                                 >  2:1657254/1‑148 (MQ=255)
                                        cacacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCt                                                                                                                 <  1:1657254/148‑1 (MQ=255)
                                         acacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGTTCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAGGCGCCGGGACTggg                                                                                                              >  1:992199/1‑150 (MQ=255)
                                         acacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTg                                                                                                                      <  2:245590/142‑1 (MQ=255)
                                         acacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCgg                                                                                                              >  1:970004/1‑150 (MQ=255)
                                         acacaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCgg                                                                                                              >  2:2307838/1‑150 (MQ=255)
                                           acaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTgg                                                                                                                     >  2:581423/1‑141 (MQ=255)
                                           acaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGtt                                                                                                            >  1:125333/1‑150 (MQ=255)
                                           acaGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGtt                                                                                                            <  2:125333/150‑1 (MQ=255)
                                              gggCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGTGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCg                                                                                                         <  1:1145558/150‑1 (MQ=255)
                                              gggCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCg                                                                                                         <  1:1647005/150‑1 (MQ=255)
                                                gcgcTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGt                                                                                                       <  2:2270660/150‑1 (MQ=255)
                                                  gcTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCt                                                                                                     <  2:1219711/150‑1 (MQ=255)
                                                  gcTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGCCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCt                                                                                                     >  2:1598445/1‑150 (MQ=255)
                                                    tGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGc                                                                                                   >  2:1733827/1‑150 (MQ=255)
                                                    tGGCCGCCCCCCAACTGTTCACGAAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGc                                                                                                   >  1:1674906/1‑150 (MQ=255)
                                                     ggCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCa                                                                                                  >  1:828345/1‑150 (MQ=255)
                                                        cgccCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCCGCTGTTCGCCGATGGTCAGTTCGAACTTGATCTAGCCGGTGACTTGCACGCGCGCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACgg                                                                                               >  2:1073653/1‑150 (MQ=255)
                                                        cgccCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACgg                                                                                               <  1:2213085/150‑1 (MQ=255)
                                                         gccCCCCACCTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGc                                                                                              <  1:184937/150‑1 (MQ=255)
                                                          ccccccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACGTGCACGCCCTGAGCGCGTGCaag                                                                                                                                        >  2:1768638/1‑105 (MQ=255)
                                                           cccccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGa                                                                                            <  2:1826626/150‑1 (MQ=255)
                                                           cccccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGa                                                                                            <  2:208272/150‑1 (MQ=255)
                                                            ccccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGAt                                                                                           <  1:357032/150‑1 (MQ=255)
                                                            ccccAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGAt                                                                                           >  1:1477548/1‑150 (MQ=255)
                                                                aaCTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCTCCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAgg                                                                                       <  1:662666/150‑1 (MQ=255)
                                                                aaCTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAgg                                                                                       <  1:1733827/150‑1 (MQ=255)
                                                                aaCTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAgg                                                                                       <  1:2307838/150‑1 (MQ=255)
                                                                    gTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCa                                                                                   >  1:1158537/1‑150 (MQ=255)
                                                                    gTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCa                                                                                   >  2:1476788/1‑150 (MQ=255)
                                                                     ttCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCAt                                                                                  >  1:717866/1‑150 (MQ=255)
                                                                     ttCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCAt                                                                                  >  1:2267500/1‑150 (MQ=255)
                                                                      tCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCCCTGAGCGCGTGCaa                                                                                                                                         <  2:252786/94‑2 (MQ=255)
                                                                      tCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCCCTGAGCGCGTGCaa                                                                                                                                         >  1:252786/1‑93 (MQ=255)
                                                                       cACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTAGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCAtct                                                                                <  1:1106434/150‑1 (MQ=255)
                                                                       cACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCCCAGATGCCGGGAGGGGAAGTCGGGGGCTCGGGGACGGGCGCGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTGCGGTCTGCACGGCGATCAGGCTCAtct                                                                                >  1:647042/1‑150 (MQ=255)
                                                                        aCGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGaa                                                                                                                            >  1:923894/1‑105 (MQ=255)
                                                                        aCGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTc                                                                                                                >  2:1929957/1‑117 (MQ=255)
                                                                        aCGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTc                                                                                                                <  1:1929957/117‑1 (MQ=255)
                                                                          gCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCgg                                                                             <  2:1776608/150‑1 (MQ=255)
                                                                              tGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGTTTTCg                                                                                                         >  1:1947970/1‑118 (MQ=255)
                                                                              tGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCg                                                                                                         <  2:1947970/118‑1 (MQ=255)
                                                                                 ccTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCCGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCCCATCTCGGCCAGCAt                                                                      <  1:1716177/150‑1 (MQ=255)
                                                                                   tCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAg                                                                                                                                     <  2:1447275/85‑1 (MQ=255)
                                                                                   tCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAg                                                                                                                                     >  1:1447275/1‑85 (MQ=255)
                                                                                   tCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCgg                                                                             <  1:974564/141‑1 (MQ=255)
                                                                                   tCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCgg                                                                             >  2:974564/1‑141 (MQ=255)
                                                                                        cgTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCgg                                                               <  2:2422115/150‑1 (MQ=255)
                                                                                           ggCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCCCATCTCGGCCAGCATCGGCCGGGTc                                                            <  1:1077064/150‑1 (MQ=255)
                                                                                                 agcTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTg                                                      <  1:1476788/150‑1 (MQ=255)
                                                                                                  gcTGTTCGGCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTgg                                                     >  2:90965/1‑150 (MQ=255)
                                                                                                     gTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGcca                                                  >  2:15909/1‑148 (MQ=255)
                                                                                                        cgtcgATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAgc                                               <  2:970004/150‑1 (MQ=255)
                                                                                                         gtcgATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCgg                                                               >  1:395896/1‑133 (MQ=255)
                                                                                                           cgATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGcc                                            <  2:833049/150‑1 (MQ=255)
                                                                                                            gATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGccc                                           >  1:134557/1‑150 (MQ=255)
                                                                                                                gTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCCATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAg                                       <  1:661961/150‑1 (MQ=255)
                                                                                                                      gTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACg                                  >  2:1968564/1‑150 (MQ=255)
                                                                                                                      gTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACg                                  >  1:330265/1‑150 (MQ=255)
                                                                                                                      gTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACg                                  >  1:836314/1‑150 (MQ=255)
                                                                                                                            cTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCg                                                                >  2:472520/1‑114 (MQ=255)
                                                                                                                            cTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCg                                                                <  1:472520/114‑1 (MQ=255)
                                                                                                                            cTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAACGCGCCCATAGCCACGGGCCGa                            <  2:956562/150‑1 (MQ=255)
                                                                                                                            cTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGa                            <  2:2267500/150‑1 (MQ=255)
                                                                                                                             ttGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAg                           <  1:915693/150‑1 (MQ=255)
                                                                                                                               gATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAgcg                         <  2:1063428/150‑1 (MQ=255)
                                                                                                                                aTCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAgcgc                        >  1:1579031/1‑150 (MQ=255)
                                                                                                                                aTCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAgcgc                        <  2:647042/150‑1 (MQ=255)
                                                                                                                                  cGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTc                      <  2:828345/150‑1 (MQ=255)
                                                                                                                                   gAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCACGCTGACCATAGCCACGGGCCGAGCGCTCa                     <  1:1520440/150‑1 (MQ=255)
                                                                                                                                   gAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCa                     <  2:717866/150‑1 (MQ=255)
                                                                                                                                           ctgggtGCACGCCCTGAGCGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAgc              >  2:1153857/6‑150 (MQ=255)
                                                                                                                                               ttGCACGCGCGCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTGCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAgcgcgc          >  2:1675695/1‑150 (MQ=255)
                                                                                                                                                  cACGCGCTCAGGGCGTGCACCCAGGTCTCGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCAtt       >  2:728276/1‑150 (MQ=255)
                                                                                                                                                   aCGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTg      <  2:1993889/150‑1 (MQ=255)
                                                                                                                                                   aCGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTg      >  1:1342150/1‑150 (MQ=255)
                                                                                                                                                     gcgcTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTGGc    >  2:1588797/1‑150 (MQ=255)
                                                                                                                                                       gcTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTGGCCa  >  2:2385958/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
GATCAAGGCATCTTCACCGTCGTGGGTGCTGGCGGCAATCCACACAGGGCGCTGGCCGCCCCCCAACTGTTCACGCAATGCCCTCGCGCGTGGCAGCAGCTGTTCGTCGATGGTCAGGTCGAACTTGATCGAGCCGGTGACTTGCACGCGCTCAGGGCGTGCACCCAGGTCACGGAAGCGCTGGGCCTCGGTTTCGGTCTGCACGGCGATCAGGCTCATCTCGGCCAGCATCGGCCGGGTCAGCCTGGCAAAGCGCCCATAGCCACGGGCCGAGCGCTCAGACAAGCGCGCATTGGCCA  >  NC_002947/5608896‑5609194

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: