Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,521,032 T→A 5.7% L83M (TTG→ATG) ‡ murD → UDP‑N‑acetylmuramoylalanine‑‑D‑glutamate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,521,0320TA 5.7% 332.9 / 8.5 106L83M (TTG→ATG) ‡murDUDP‑N‑acetylmuramoylalanine‑‑D‑glutamate ligase
Reads supporting (aligned to +/- strand):  ref base T (56/44);  new base A (4/2);  total (60/46)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

Only 100 of 118 total aligned reads displayed.

TTGCCTTTGCGGTCGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATG  >  NC_002947/1520883‑1521181
                                                                                                                                                     |                                                                                                                                                     
ttGCCTTTGCGGTCGCCGATACCCGCGAGCAGCCGCCGGAAATGGCTGTACGCCGCCTTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCggcctggcct                                                                                                                                                       <  2:664127/150‑1 (MQ=255)
 tGCCTTTGCGGTCGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCtt                                                                                                                                                      <  2:2171025/150‑1 (MQ=255)
  gCCTTTGCGGTCGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGCAGTTCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTg                                                                                                                                                     <  1:1345397/150‑1 (MQ=255)
  gCCTTTGCGGTCGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTg                                                                                                                                                     >  2:2190178/1‑150 (MQ=255)
    cTTTGCGGTCGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTg                                                                                                                                                     <  1:35445/148‑1 (MQ=255)
    cTTTGCGGTCGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTg                                                                                                                                                     >  2:35445/1‑148 (MQ=255)
    cTTTGCGGTCGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGc                                                                                                                                                   >  1:1467433/1‑150 (MQ=255)
          ggTCGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGcc                                                                                                                                             >  1:421991/1‑150 (MQ=255)
             cGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGc                                                                                                                                          <  2:1930353/150‑1 (MQ=255)
              gCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGcc                                                                                                                                         >  1:1490564/1‑150 (MQ=255)
               ccGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCt                                                                                                                                        >  2:831835/1‑150 (MQ=255)
                cGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCAGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCtt                                                                                                                                       >  2:831439/1‑150 (MQ=255)
                  aTACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTgc                                                                                                                                     <  2:2362904/150‑1 (MQ=255)
                      ccGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGTACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTgcagca                                                                                                                                 >  2:1381653/1‑150 (MQ=255)
                      ccGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTg                                                                                                                                                     >  1:1042323/1‑130 (MQ=255)
                      ccGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTg                                                                                                                                                     <  2:1042323/130‑1 (MQ=255)
                      ccGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTgcagca                                                                                                                                 <  2:1904221/150‑1 (MQ=255)
                      ccGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTgcagca                                                                                                                                 <  1:1131818/150‑1 (MQ=255)
                      ccGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTgcagca                                                                                                                                 >  1:2305785/1‑150 (MQ=255)
                      ccGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTgcagca                                                                                                                                 >  1:481452/1‑150 (MQ=255)
                            gcagcCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGc                                                                                                                           >  1:1603342/1‑150 (MQ=255)
                                ccgccgGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCtt                                                                                                                                                      <  2:548547/119‑1 (MQ=255)
                                ccgccgGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCtt                                                                                                                                                      >  1:548547/1‑119 (MQ=255)
                                  gccgGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGCGGTGGGGTGGTGGGAACTGGACGTGGACTTCCATTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCcgcgcg                                                                                                                     >  1:415064/1‑150 (MQ=255)
                                  gccgGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTgg                                                                                                                                                    <  1:1406173/119‑1 (MQ=255)
                                  gccgGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCGGCAGCCTGctgcaag                                                                                                                     <  2:2020477/150‑7 (MQ=255)
                                  gccgGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCAGGCCTTgg                                                                                                                                                    >  2:1406173/1‑119 (MQ=255)
                                    cgGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCgcgtg                                                                                                                   >  1:1467171/1‑150 (MQ=255)
                                     ggAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTTCAGCAGGCTGCCGCgcgtgg                                                                                                                  >  2:2259466/1‑150 (MQ=255)
                                      gAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCgcgtggc                                                                                                                 >  1:1459053/1‑150 (MQ=255)
                                          tGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGa                                                                                                             >  1:1136022/1‑150 (MQ=255)
                                           ggAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGaa                                                                                                            <  2:1467171/150‑1 (MQ=255)
                                            gAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCgcgtgg                                                                                                                  >  2:2447011/1‑143 (MQ=255)
                                            gAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCgcgtgg                                                                                                                  <  1:2447011/143‑1 (MQ=255)
                                              aaCCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACCCGACGTGCGCTCCCGGTGCCGCGCCAACGAGCGGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCt                                                                                                         <  2:1103988/150‑1 (MQ=255)
                                                  cTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAg                                                                                                           >  2:2065861/1‑144 (MQ=255)
                                                  cTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAg                                                                                                           <  1:2065861/144‑1 (MQ=255)
                                                  cTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCt                                                                                                     >  2:1304864/1‑150 (MQ=255)
                                                  cTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCt                                                                                                     >  2:397400/1‑150 (MQ=255)
                                                      gcCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCg                                                                                                 >  2:1126337/1‑150 (MQ=255)
                                                       ccGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGa                                                                                                >  1:740248/1‑150 (MQ=255)
                                                       ccGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGa                                                                                                >  1:154255/1‑150 (MQ=255)
                                                            aTTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATag                                                                                           >  2:1571544/1‑150 (MQ=255)
                                                             ttACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGaa                                                                                                            >  2:84103/1‑132 (MQ=255)
                                                             ttACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATaga                                                                                          <  1:319067/150‑1 (MQ=255)
                                                             ttACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATaga                                                                                          >  2:1393437/1‑150 (MQ=255)
                                                             ttACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGACGGCTTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATaga                                                                                          >  1:744032/1‑150 (MQ=255)
                                                                 ccGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGTACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGa                                                                                                <  1:882722/140‑1 (MQ=255)
                                                                 ccGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGTACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGa                                                                                                >  2:882722/1‑140 (MQ=255)
                                                                 ccGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTg                                                                                      >  1:1794012/1‑150 (MQ=255)
                                                                   gCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGtt                                                                                    <  1:1791698/150‑1 (MQ=255)
                                                                       gTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGt                                                                                     >  2:408200/1‑145 (MQ=255)
                                                                       gTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGt                                                                                     <  1:408200/145‑1 (MQ=255)
                                                                        tGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTc                                                                               <  1:852070/150‑1 (MQ=255)
                                                                         ggAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCg                                                                              >  1:2444201/1‑150 (MQ=255)
                                                                         ggAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCg                                                                              <  2:2304129/150‑1 (MQ=255)
                                                                          gAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGc                                                                             <  1:2162786/150‑1 (MQ=255)
                                                                                   tgtgGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAg                                                                    <  2:2451588/150‑1 (MQ=255)
                                                                                       gtgAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCAAGGCCAGGCCgggg                                                                                                                                         >  1:1154436/1‑75 (MQ=38)
                                                                                       gtgAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCAAGGCCAGGCCgggg                                                                                                                                         <  2:1154436/77‑3 (MQ=38)
                                                                                           aCTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGAt                                                            <  2:421991/150‑1 (MQ=255)
                                                                                             tGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCgg                                                                                                                                           <  2:1071280/69‑1 (MQ=255)
                                                                                             tGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCgg                                                                                                                                           >  1:1071280/1‑69 (MQ=255)
                                                                                               gACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGcagcag                                                                                                                                >  2:460146/1‑78 (MQ=255)
                                                                                               gACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGcagcag                                                                                                                                <  1:460146/78‑1 (MQ=255)
                                                                                               gACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCgcgtgg                                                                                                                  >  1:1575650/1‑92 (MQ=255)
                                                                                               gACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCgcgtgg                                                                                                                  <  2:1575650/92‑1 (MQ=255)
                                                                                               gACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGtt                                                        <  1:831439/150‑1 (MQ=255)
                                                                                               gACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCCACATAGCGCTGTTTGCCCGCCATGCCAAGGCCCCGCTCGtt                                                        >  2:281151/1‑150 (MQ=255)
                                                                                                     gACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATc                                                  <  2:1136022/150‑1 (MQ=255)
                                                                                                      aCTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATaga                                                                                          >  2:1930370/1‑109 (MQ=255)
                                                                                                      aCTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATaga                                                                                          <  1:1930370/109‑1 (MQ=255)
                                                                                                       cTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGGGGCGTGGAGCTGTCGGGCTACGTAGGGCCGTTTTCGCGCCATGCCAAGGCGCCGATCCTTGccgccag                                                >  2:500794/1‑145 (MQ=255)
                                                                                                          ccTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTgg                                             <  1:1571544/150‑1 (MQ=255)
                                                                                                            tgtgCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGtt                                           >  2:637105/1‑150 (MQ=255)
                                                                                                                ccGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCaa                                       >  2:1789024/1‑150 (MQ=255)
                                                                                                                ccGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCaa                                       >  2:40300/1‑150 (MQ=255)
                                                                                                                  gcgcCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGAGGTGTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACg                                     >  1:1329995/1‑150 (MQ=255)
                                                                                                                     ccAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATc                                                           >  2:2418192/1‑125 (MQ=255)
                                                                                                                      cAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTc                                          <  2:1377640/141‑1 (MQ=255)
                                                                                                                      cAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTc                                          >  1:1377640/1‑141 (MQ=255)
                                                                                                                             cTGTACGTGAGCCCCGGCCTGTCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCAcc                          <  1:1304864/150‑1 (MQ=255)
                                                                                                                             cTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCAcc                          >  2:2372719/1‑150 (MQ=255)
                                                                                                                              tGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCg                         >  2:1011535/1‑150 (MQ=255)
                                                                                                                                tACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCTATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTc                       >  2:1458514/1‑150 (MQ=255)
                                                                                                                                  cGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGc                                                                                                  <  1:1050397/73‑1 (MQ=255)
                                                                                                                                   tgcaaggccGGCGTGGCCAAGGCCAGGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGAGATATAGCGGTGTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCacg                    >  1:290007/8‑150 (MQ=255)
                                                                                                                                    tGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCacga                   <  1:831835/150‑1 (MQ=255)
                                                                                                                                        ggccGGCGTGGCCAAGGCCAGGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTc                                                                               >  1:136254/3‑86 (MQ=37)
                                                                                                                                        ggccGGCGTGGCCAAGGCCAGGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTc                                                                               <  2:136254/84‑1 (MQ=37)
                                                                                                                                        ggccGGCGTGGCCAAGGCCAGGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCt               >  2:1036930/3‑150 (MQ=255)
                                                                                                                                        ccccGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCt               >  2:2184572/1‑150 (MQ=255)
                                                                                                                                        ccccGGCATGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGAGACACGGTGTCTGGCGACGATCTTCTGTTTGCTCTCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACACCAAGAGCACCGTCACGACGCt               >  1:290983/1‑150 (MQ=255)
                                                                                                                                         cccGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTg              <  2:1603342/150‑1 (MQ=255)
                                                                                                                                               ctggcctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGt        <  2:1987441/150‑1 (MQ=255)
                                                                                                                                                tggcctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTAATACGGCCTATGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTg       >  1:746429/1‑150 (MQ=255)
                                                                                                                                                   cctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGaaa    <  2:690659/150‑1 (MQ=255)
                                                                                                                                                    ctTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAAt   <  1:1315547/150‑1 (MQ=255)
                                                                                                                                                    attGGCCACGACGAACTTGCAGCAAGCTACCGCACGTGGCGTTAATCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGACCGTTGCGAAGAGCGTCGCGAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAAt   >  2:2366342/2‑150 (MQ=255)
                                                                                                                                                     ttGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATg  <  1:2372719/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTGCCTTTGCGGTCGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATG  >  NC_002947/1520883‑1521181

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: