Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,521,039 C→G 5.4% T85R (ACG→AGG)  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,0390CG 5.4% 352.2 / 9.3 111T85R (ACG→AGG) murDUDP‑N‑acetylmuramoylalanine‑‑D‑glutamate ligase
Reads supporting (aligned to +/- strand):  ref base C (56/49);  new base G (4/2);  total (60/51)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

Only 100 of 111 total aligned reads displayed.

GGTCGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGC  >  NC_002947/1520893‑1521186
                                                                                                                                                  |                                                                                                                                                   
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)
            ccGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTgcagca                                                                                                                                      <  1:1131818/150‑1 (MQ=255)
            ccGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTgcagca                                                                                                                                      >  1:481452/1‑150 (MQ=255)
            ccGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTgcagca                                                                                                                                      >  1:2305785/1‑150 (MQ=255)
            ccGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTgcagca                                                                                                                                      <  2:1904221/150‑1 (MQ=255)
                        gccgGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGCGGTGGGGTGGTGGGAACTGGACGTGGACTTCCATTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCcgcgcg                                                                                                                          >  1:415064/1‑150 (MQ=255)
                        gccgGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCGGCAGCCTGctgcaag                                                                                                                          <  2:2020477/150‑7 (MQ=255)
                          cgGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCgcgtg                                                                                                                        >  1:1467171/1‑150 (MQ=255)
                           ggAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAAGGAGCAGTAGGTGGGTCCCGGACTGGCTTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCgcgtgg                                                                                                                       <  2:683933/150‑1 (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                                                                                                                       <  1:2447011/143‑1 (MQ=255)
                                    aaCCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACCCGACGTGCGCTCCCGGTGCCGCGCCAACGAGCGGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCt                                                                                                              <  2:1103988/150‑1 (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)
                                            gcCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCg                                                                                                      <  1:1188250/150‑1 (MQ=255)
                                             ccGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGa                                                                                                     >  1:740248/1‑150 (MQ=255)
                                             ccGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGa                                                                                                     >  1:154255/1‑150 (MQ=255)
                                                   ttACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGaa                                                                                                                 >  2:84103/1‑132 (MQ=255)
                                                   ttACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATaga                                                                                               >  2:1393437/1‑150 (MQ=255)
                                                   ttACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATaga                                                                                               <  1:319067/150‑1 (MQ=255)
                                                   ttACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGACGGCTTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATaga                                                                                               >  1:744032/1‑150 (MQ=255)
                                                   ttACCCGCAGGGGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGaa                                                                                                                 <  1:84103/132‑1 (MQ=255)
                                                    tACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGAAATagag                                                                                              >  2:257666/1‑150 (MQ=255)
                                                       ccGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGTACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGa                                                                                                     >  2:882722/1‑140 (MQ=255)
                                                       ccGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTg                                                                                           >  1:1794012/1‑150 (MQ=255)
                                                        cGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGt                                                                                          <  2:1794012/150‑1 (MQ=255)
                                                         gCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGtt                                                                                         <  1:1791698/150‑1 (MQ=255)
                                                             gTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGt                                                                                          <  1:408200/145‑1 (MQ=255)
                                                             gTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGt                                                                                          >  2:408200/1‑145 (MQ=255)
                                                              tGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTc                                                                                    <  2:2305785/150‑1 (MQ=255)
                                                              tGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTc                                                                                    <  1:852070/150‑1 (MQ=255)
                                                               ggAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCg                                                                                   <  2:2304129/150‑1 (MQ=255)
                                                                gAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGc                                                                                  <  1:2162786/150‑1 (MQ=255)
                                                                         tgtgGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAg                                                                         <  2:2451588/150‑1 (MQ=255)
                                                                             gtgAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCAAGGCCAGGCCgggg                                                                                                                                              <  2:1154436/77‑3 (MQ=38)
                                                                             gtgAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCAAGGCCAGGCCgggg                                                                                                                                              >  1:1154436/1‑75 (MQ=38)
                                                                                 aCTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGAt                                                                 <  2:421991/150‑1 (MQ=255)
                                                                                   tGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCgg                                                                                                                                                <  2:1071280/69‑1 (MQ=255)
                                                                                   tGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCgg                                                                                                                                                >  1:1071280/1‑69 (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:1529084/1‑150 (MQ=255)
                                                                                     gACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGtt                                                             <  1:831439/150‑1 (MQ=255)
                                                                                     gACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCCACATAGCGCTGTTTGCCCGCCATGCCAAGGCCCCGCTCGtt                                                             >  2:281151/1‑150 (MQ=255)
                                                                                        gTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGcc                                                          >  2:2120947/1‑150 (MQ=255)
                                                                                           gACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATc                                                       <  2:1136022/150‑1 (MQ=255)
                                                                                            aCTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATaga                                                                                               <  1:1930370/109‑1 (MQ=255)
                                                                                            aCTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATaga                                                                                               >  2:1930370/1‑109 (MQ=255)
                                                                                             cTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGGGGCGTGGAGCTGTCGGGCTACGTAGGGCCGTTTTCGCGCCATGCCAAGGCGCCGATCCTTGccgccag                                                     >  2:500794/1‑145 (MQ=255)
                                                                                                ccTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTgg                                                  <  1:1571544/150‑1 (MQ=255)
                                                                                                 cTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGt                                                 <  1:397400/150‑1 (MQ=255)
                                                                                                  tgtgCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGtt                                                >  2:637105/1‑150 (MQ=255)
                                                                                                      ccGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCaa                                            >  2:40300/1‑150 (MQ=255)
                                                                                                      ccGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCaa                                            >  2:1789024/1‑150 (MQ=255)
                                                                                                        gcgcCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGAGGTGTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACg                                          >  1:1329995/1‑150 (MQ=255)
                                                                                                           ccAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATc                                                                <  1:2418192/125‑1 (MQ=255)
                                                                                                           ccAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATc                                                                >  2:2418192/1‑125 (MQ=255)
                                                                                                            cAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTc                                               <  2:1377640/141‑1 (MQ=255)
                                                                                                            cAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCaa                                      <  1:1839188/150‑1 (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)
                                                                                                                      tACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTc                            >  1:1398869/1‑150 (MQ=255)
                                                                                                                        cGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGc                                                                                                       <  1:1050397/73‑1 (MQ=255)
                                                                                                                        cGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGc                                                                                                       >  2:1050397/1‑73 (MQ=255)
                                                                                                                         tgcaaggccGGCGTGGCCAAGGCCAGGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGAGATATAGCGGTGTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCacg                         >  1:290007/8‑150 (MQ=255)
                                                                                                                          tGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCacga                        <  1:831835/150‑1 (MQ=255)
                                                                                                                              ggccGGCGTGGCCAAGGCCAGGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTc                                                                                    <  2:136254/84‑1 (MQ=37)
                                                                                                                              ggccGGCGTGGCCAAGGCCAGGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTc                                                                                    >  1:136254/3‑86 (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:1459053/150‑1 (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)
                                                                                                                                           ttGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATg       <  2:415064/150‑1 (MQ=255)
                                                                                                                                             ggCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCCTCAGTGGTTCCACCGCCACGGGCACCGGCTCGACGCTGGTCGGTGAAATGgc     <  2:1467433/150‑1 (MQ=255)
                                                                                                                                              gCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGgcc    <  2:154255/150‑1 (MQ=255)
                                                                                                                                              gCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGgcc    >  2:678920/1‑150 (MQ=255)
                                                                                                                                                cACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCa                           >  1:403341/1‑125 (MQ=255)
                                                                                                                                                cACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCa                           <  2:403341/125‑1 (MQ=255)
                                                                                                                                                cACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGgccgc  <  1:1381653/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GGTCGCCGATACCCGCGAGCAGCCGCCGGAACTGGAAACCCTGCGCCGTGATTACCCGCAGGTGGAAGTGCGCTGTGGTGAACTGGACGTGGACTTCCTGTGCCGCGCCAACGAGCTGTACGTGAGCCCCGGCCTGGCCTTGGCCACGCCGGCCTTGCAGCAGGCTGCCGCGCGTGGCGTGAAGCTGTCTGGCGACATAGAGCTGTTTGCTCGCCATGCCAAGGCGCCGATCGTTGCCATCAGTGGTTCCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGC  >  NC_002947/1520893‑1521186

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