Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,457,256 T→G 5.3% intergenic (+98/+58) mpl → / ← yjgA UDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase/ribosome‑associated UPF0307 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,457,2560TG 5.3% 223.1 / 2.7 86intergenic (+98/+58)mpl/yjgAUDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase/ribosome‑associated UPF0307 family protein
Reads supporting (aligned to +/- strand):  ref base T (46/35);  new base G (2/3);  total (48/38)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.99e-01

ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTCTCTTTCTTCGCG  >  NC_000913/4457109‑4457404
                                                                                                                                                   |                                                                                                                                                    
aTCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTTAGCATAGTGCCTGATGCGACCCTTGATGAGTCtt                                                                                                                                                     >  2:1182484/1‑149 (MQ=255)
aTCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCtt                                                                                                                                                     <  1:158819/149‑1 (MQ=255)
aTCAGAAACTGCTGGAAGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGCCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCtt                                                                                                                                                     >  1:233793/1‑149 (MQ=255)
 tCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTa                                                                                                                                                    >  1:286289/1‑149 (MQ=255)
  cAGTAACTGCTGGATGGGCTGGCTAAGAAGGCGGAAGCGGCGCAGTAATTCGGCGTTAGGCAGATATAGTATGCGTAGGCAGGATAAGGCGTTGACGCCGCATCCGGCATTTGAGCATAGTGCCTGTTTCGACTCTTGATGCGTCTTAt                                                                                                                                                   <  1:727281/149‑1 (MQ=37)
  cAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTAt                                                                                                                                                   <  2:1021207/149‑1 (MQ=255)
    gAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCCCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCa                                                                                                                                                 >  1:60867/1‑149 (MQ=255)
     aaaCTGCTGTATGGGCTGGCGACGACGGCGGAAGCCGCGCAGTAATTCGGCATTAGCCAGATATAGTATGCGTAGTCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAg                                                                                                                                                <  1:794705/149‑1 (MQ=38)
      acctactggctggGCTGGCGAAGAAGGCGGAAGCCGCGCAGTATTTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTGCACGCCGGATCCGGCATTAGAGCATAGTGCCTGATAAGACGCATCAAGCGTCGCATCAgg                                                                                                                                               <  1:265605/144‑1 (MQ=255)
       aCTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGc                                                                                                                                              >  2:926556/1‑149 (MQ=255)
        ctgctgGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGca                                                                                                                                             >  1:902759/1‑148 (MQ=255)
        ctgctgGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATAAGACGCATCAAGCGTCGCATCAGGca                                                                                                                                             <  1:736806/149‑2 (MQ=38)
          gctgGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTa                                                                                                                                           <  1:1159539/149‑1 (MQ=255)
           ctgGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCCTAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTAc                                                                                                                                          >  2:218249/1‑149 (MQ=255)
              gATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAgg                                                                                                                                       >  2:906933/1‑149 (MQ=255)
               attggCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGt                                                                                                                                      <  1:450414/146‑1 (MQ=255)
                tGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTg                                                                                                                                     <  2:356644/149‑1 (MQ=255)
                  ggCTGGCGATGAAGGCGGAAGCCGCGCATTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCa                                                                                                                                   <  1:1038183/149‑1 (MQ=255)
                   gCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCCTCCGGCATTTGAGCATAGTGCCTGATGCGCCGCTTGATGCGTCTTATCAGGCCTACAGGTGCAc                                                                                                                                  >  2:876026/1‑149 (MQ=255)
                    cTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGAGGCGTCTTATCAGGCCTACAGGTGCAcc                                                                                                                                 >  2:125423/1‑149 (MQ=255)
                        cGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTAt                                                                                                                             <  1:101895/149‑1 (MQ=255)
                         gaagaaGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATc                                                                                                                            >  1:744611/1‑149 (MQ=255)
                         gaagaaGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTCATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATc                                                                                                                            >  1:808009/1‑149 (MQ=255)
                            gcagGCGGAAGCCGCGCAGTAATGCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCgg                                                                                                                         <  2:182799/147‑1 (MQ=255)
                             aagGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACAGTATCCGGc                                                                                                                        >  1:567067/1‑149 (MQ=255)
                               ggCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCaa                                                                                                                      <  1:904195/149‑1 (MQ=255)
                                  ggAAGCCGCGCAGTAATTAGGCCTGAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAAcc                                                                                                                   >  1:694254/1‑149 (MQ=255)
                                    aaGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCAt                                                                                                                 <  2:1023537/149‑1 (MQ=255)
                                        cgcgCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCCTCCGGCATTTGAGCATAGTGCCTGCTGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTc                                                                                                             >  2:920906/1‑149 (MQ=255)
                                        cgcgCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTc                                                                                                             >  1:553231/1‑149 (MQ=255)
                                             aGTAATTCGGCCTTAGCCAGATATAGTATTCGTAGGCCGTATAAGGCGTGCACGCCGCATCCTGCATATGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGGGCACCGTATCCGGCAAACCATACTCACATc                                                                                                        <  2:740017/149‑1 (MQ=38)
                                             aGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATc                                                                                                        >  2:842350/1‑149 (MQ=255)
                                                   tCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCaacaac                                                                                                  >  1:155277/1‑149 (MQ=255)
                                                     ggCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGa                                                                                                >  2:199662/1‑149 (MQ=255)
                                                        cTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGaaaa                                                                                             >  1:1008460/1‑149 (MQ=255)
                                                         ttAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAAt                                                                                            >  1:227856/1‑149 (MQ=255)
                                                             ccAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAATTTAcc                                                                                        >  2:513776/1‑149 (MQ=255)
                                                             ccAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTAcc                                                                                        >  1:1153013/1‑149 (MQ=255)
                                                             ccAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTAcc                                                                                        >  1:160610/1‑149 (MQ=255)
                                                               aGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATGAGGCCTACAGGTGCACCGTCTCCGGCAAACCATACTCACATCAACAACGAAAATGACCCt                                                                                      >  2:1178225/1‑149 (MQ=255)
                                                                     aGTATGCGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTtc                                                                                <  1:724952/149‑1 (MQ=255)
                                                                       tATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTtctc                                                                              <  2:456938/149‑1 (MQ=255)
                                                                       tATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTtctc                                                                              <  2:567067/149‑1 (MQ=255)
                                                                        aTGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTcc                                                                             <  1:849496/149‑1 (MQ=255)
                                                                         tGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGCAAATTACCCTTCGTTCTCCg                                                                            >  2:861771/1‑149 (MQ=255)
                                                                          tcaaaagCCGGATGCGGCGTGAACGCCTTAACCGGCCGTAGAGCATATTGCCTGATGCGACGCTTGCTGCGTCTTATCATGCCTACAGGTGCACCGTATCCGGCAATCAATACTAACATCAACAACGAAAATTACCCATCGTTATACGc                                                                           >  2:653009/7‑149 (MQ=19)
                                                                            aaatgCCGGATGCGGCGTGAACGCCTTATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTa                                                                         >  1:224904/5‑149 (MQ=38)
                                                                             tAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTaa                                                                        >  2:309330/1‑149 (MQ=255)
                                                                               ggCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACt                                                                      >  1:651228/1‑149 (MQ=255)
                                                                                gCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTc                                                                     <  2:892144/149‑1 (MQ=255)
                                                                                  cGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTcgc                                                                   >  2:579826/1‑149 (MQ=255)
                                                                                    gATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCgcgc                                                                 <  2:744611/149‑1 (MQ=255)
                                                                                    gATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCgcgc                                                                 <  1:161128/149‑1 (MQ=255)
                                                                                    gATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCgcgc                                                                 <  1:260712/149‑1 (MQ=255)
                                                                                    gATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCgcgc                                                                 <  2:808009/149‑1 (MQ=255)
                                                                                    gATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCgcgc                                                                 <  2:845250/149‑1 (MQ=255)
                                                                                       aaGGCGGTCACGCCGCATCCGGCATTTGAGCATAGTGGCTGCTGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGa                                                              <  2:770888/149‑1 (MQ=255)
                                                                                          gCGTTCACGCCGCGTCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTGATCCGGCCTACAGGTGCACCGTATCCGGCCATCCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTCACTCGCGCAGATAc                                                           >  2:778727/1‑149 (MQ=39)
                                                                                          gCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATAc                                                           >  1:347300/1‑149 (MQ=255)
                                                                                          gCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTAGCCGGCAAACCATACTCACATCAACAACGAAATTTACCCTTCGTTCTCCGCTAACTCGCGGAGATAc                                                           >  2:931731/1‑149 (MQ=255)
                                                                                            gTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCAGGCAAACCATACTCACATCAACAACGCAAATTACCCTTCGATCTCCGCGAACTCGCGCAGATACTg                                                         >  1:215586/1‑149 (MQ=255)
                                                                                              tCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGa                                                       <  2:553231/149‑1 (MQ=255)
                                                                                                 cgccgcATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCAAATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGAAGATACTGGaaat                                                    >  2:299631/1‑148 (MQ=255)
                                                                                                    cgcATCCTGCATTTGATCATATTGCCTGATGCGATGCTTTAGGCGTCTAAGCAGGCCTACAGGTGCCCCTTATCCGGCAAACCATACTCAGATCAACAACGAAAATTACCCTTCGTTCTCCGCAAACTCGCGCAGATACTGGAAAATCt                                                 <  2:1013070/149‑1 (MQ=255)
                                                                                                    cgcATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCt                                                 >  2:831255/1‑149 (MQ=255)
                                                                                                     gcATCCTGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCATACAGGTTCACCGTTTCCGGATAACCATACTCACATCAACTACGAAAATTACCCTTCGTTCTCGGCTAACTCGCGCAGATACTGGAAATTCTg                                                <  2:697102/149‑1 (MQ=38)
                                                                                                     gcATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTg                                                <  2:638085/149‑1 (MQ=255)
                                                                                                         ccGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGt                                            >  2:633372/1‑149 (MQ=255)
                                                                                                         ccGGAATTTGAGCATAGTGCCTGATGCGCCGCTTGATGAGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCCAACCATAATAACATCAACAACTAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGt                                            >  2:917221/1‑149 (MQ=255)
                                                                                                            gCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAATTTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCg                                         >  1:799604/1‑149 (MQ=255)
                                                                                                                  gAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTc                                   >  2:29883/1‑149 (MQ=255)
                                                                                                                  gAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTc                                   <  1:659994/149‑1 (MQ=255)
                                                                                                                     cATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAATCCATACTCACCTCAACAACGATAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTCGTAAATCTGGCTTGCGGATTTcggc                                >  1:1022081/1‑149 (MQ=38)
                                                                                                                      aTAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGATCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTcggcg                               <  1:770007/149‑1 (MQ=255)
                                                                                                                        aGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCggcggc                             <  2:1008460/149‑1 (MQ=255)
                                                                                                                            ccTGATAAGACGCATCAAGCGTCGCATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTAt                         >  1:1032119/1‑149 (MQ=255)
                                                                                                                                ataagACGCATCAAGCGTCGCATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCTGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTccc                     >  1:1131912/5‑149 (MQ=255)
                                                                                                                                  gCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCATCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCctt                   <  2:927093/149‑1 (MQ=255)
                                                                                                                                    gACGCTTGATGCGTCTTATCAGGCCTACAGGTGCCCCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCgttct                 >  1:656306/1‑149 (MQ=255)
                                                                                                                                    gACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATCCTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATATCGGCGGCTTATTCCcttct                 >  1:263069/1‑149 (MQ=255)
                                                                                                                                    gACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGTCGGCTTATTCCcttct                 >  1:563643/1‑149 (MQ=255)
                                                                                                                                    gACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCcttct                 >  1:864353/1‑149 (MQ=255)
                                                                                                                                     aCGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCcttctt                <  1:926556/149‑1 (MQ=255)
                                                                                                                                        cTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCttctttct             <  2:1153013/149‑1 (MQ=255)
                                                                                                                                           taTGCGTCTTATCATGCCTATAGGTGCACCGTATACGGCAAACCATACTCACATCAACAACTAAAATTACACTTCGTTCTCCGCTAACTAGCGAAGAAACTGGAAAATCATGCGTGCGGATTTCTTATGCTTAATCACTtatttctcat          >  2:382706/2‑147 (MQ=255)
                                                                                                                                           gATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTACCTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTtctttctctt          <  2:356087/149‑1 (MQ=255)
                                                                                                                                                gTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTGTTTCTCTttcttc     >  2:1018501/1‑149 (MQ=255)
                                                                                                                                                gTCGTATCACGCCTACAGGTGCACCGTCTCCGGCAAACCATACTCACATCAACACCGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTCTCTttcttc     <  2:1095953/149‑1 (MQ=255)
                                                                                                                                                  cTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTCTCTTTCTTcgc   <  1:199662/149‑1 (MQ=255)
                                                                                                                                                   ttATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTCTCTTTCTTcgcg  >  2:508658/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATCAGAAACTGCTGGATGGGCTGGCGAAGAAGGCGGAAGCCGCGCAGTAATTCGGCCTTAGCCAGATATAGTATGCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATTTGAGCATAGTGCCTGATGCGACGCTTGATGCGTCTTATCAGGCCTACAGGTGCACCGTATCCGGCAAACCATACTCACATCAACAACGAAAATTACCCTTCGTTCTCCGCTAACTCGCGCAGATACTGGAAAATCTGGCGTGCGGATTTCGGCGGCTTATTCCCTTCTTTCTCTTTCTTCGCG  >  NC_000913/4457109‑4457404

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