Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,445,719 A→T 90.0% T40S (ACC→TCC)  mpl → UDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,445,7190AT90.0% 681.9 / 54.1 228T40S (ACC→TCC) mplUDP‑N‑acetylmuramate:L‑alanyl‑gamma‑D‑glutamyl‑ meso‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base A (8/15);  new base T (104/101);  total (112/116)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 229 total aligned reads displayed.

GAGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTT  >  CP009273/4445572‑4445863
                                                                                                                                                   |                                                                                                                                                
  gCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTcc                                                                                                                                                <  1:3107637/148‑1 (MQ=255)
  gCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTcc                                                                                                                                                <  2:2982085/148‑1 (MQ=255)
   cGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCt                                                                                                                                               >  2:1306067/1‑148 (MQ=255)
    ggCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCtt                                                                                                                                              <  2:1615912/148‑1 (MQ=255)
      cAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTAc                                                                                                                                            >  1:779512/1‑148 (MQ=255)
        gATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACtt                                                                                                                                          >  2:1428040/1‑148 (MQ=255)
         aTTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTg                                                                                                                                         >  1:2307589/1‑148 (MQ=255)
          ttAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTga                                                                                                                                        >  2:1112696/1‑148 (MQ=255)
           tAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTgag                                                                                                                                       >  2:82814/1‑148 (MQ=255)
           tAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTgag                                                                                                                                       <  1:2461460/148‑1 (MQ=255)
           tAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTgag                                                                                                                                       >  2:1914303/1‑148 (MQ=255)
             aaaaaGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGaa                                                                                                                                     >  2:3100151/1‑148 (MQ=255)
              aaaaGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAg                                                                                                                                    >  2:1785175/1‑148 (MQ=255)
               aaaGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGc                                                                                                                                   <  2:3224270/148‑1 (MQ=255)
               aaaGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGc                                                                                                                                   <  2:3224287/148‑1 (MQ=255)
                aaGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCa                                                                                                                                  <  2:1848978/148‑1 (MQ=255)
                aaGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCa                                                                                                                                  <  1:2182862/148‑1 (MQ=255)
                 aGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCaa                                                                                                                                 >  2:2513487/1‑148 (MQ=255)
                 aGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCaa                                                                                                                                 <  2:2783964/148‑1 (MQ=255)
                  ggATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAg                                                                                                                                >  2:1293774/1‑148 (MQ=255)
                  ggATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAg                                                                                                                                >  2:2999590/1‑148 (MQ=255)
                   taTAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAgg                                                                                                                               >  1:2444934/2‑148 (MQ=255)
                          aCGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAg                                                                                                                        <  2:2216971/148‑1 (MQ=255)
                              aTGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGa                                                                                                                    >  1:1158008/1‑148 (MQ=255)
                              aTGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGa                                                                                                                    >  1:1158022/1‑148 (MQ=255)
                              aTGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGa                                                                                                                    >  1:2031283/1‑148 (MQ=255)
                              aTGCGCAGTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGa                                                                                                                    <  1:173384/148‑1 (MQ=255)
                                  gcATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCa                                                                                                                <  2:1164534/148‑1 (MQ=255)
                                     ttCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAggg                                                                                                             >  1:865975/1‑148 (MQ=255)
                                      tCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAgggg                                                                                                            <  2:2443732/148‑1 (MQ=255)
                                      tCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAgggg                                                                                                            >  1:1563878/1‑148 (MQ=255)
                                      tCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAgggg                                                                                                            >  1:326146/1‑148 (MQ=255)
                                        atatTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTa                                                                                                          <  1:1989415/148‑1 (MQ=255)
                                        atatTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTa                                                                                                          >  1:842414/1‑148 (MQ=255)
                                        atatTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTa                                                                                                          <  1:81862/148‑1 (MQ=255)
                                         tatTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTAc                                                                                                         >  2:1292822/1‑148 (MQ=255)
                                          atTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACg                                                                                                        <  2:2710020/148‑1 (MQ=255)
                                           ttttAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGa                                                                                                       >  2:2197471/1‑148 (MQ=255)
                                           ttttAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGa                                                                                                       <  2:479126/148‑1 (MQ=255)
                                                    tttGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATgccagcca                                                                                              >  2:1534357/1‑148 (MQ=255)
                                                          gCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCTGTGGTACGATGCCAGCCAGCTCGa                                                                                        >  1:1607254/1‑148 (MQ=255)
                                                          gCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGa                                                                                        >  1:1607870/1‑148 (MQ=255)
                                                          gCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGa                                                                                        <  2:706891/148‑1 (MQ=255)
                                                           cacaTTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAg                                                                                       >  2:1459355/1‑148 (MQ=255)
                                                                   tGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGcc                                                                               >  1:604419/1‑148 (MQ=255)
                                                                     ggCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCgg                                                                             <  1:1070844/148‑1 (MQ=255)
                                                                     ggCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCgg                                                                             >  1:1788104/1‑148 (MQ=255)
                                                                      gCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGa                                                                            <  1:2878013/148‑1 (MQ=255)
                                                                       cGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGAt                                                                           >  1:2990433/1‑148 (MQ=255)
                                                                       cGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTCACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGAt                                                                           <  2:2761235/148‑1 (MQ=255)
                                                                       cGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCCCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGAt                                                                           >  1:2990748/1‑148 (MQ=255)
                                                                          tCTGGCGATGCTGGCGCGCCAGTTCGGCCATGAAGTAACGGGGTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGAGCTg                                                                        >  1:1455828/1‑148 (MQ=255)
                                                                          tCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTg                                                                        >  1:1957322/1‑148 (MQ=255)
                                                                          tCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTg                                                                        >  2:2829210/1‑148 (MQ=255)
                                                                          tCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTg                                                                        <  1:364995/148‑1 (MQ=255)
                                                                          tCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTg                                                                        >  1:1455864/1‑148 (MQ=255)
                                                                           cTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTgg                                                                       <  2:1922313/148‑1 (MQ=255)
                                                                                 aTGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGatta                                                                 <  1:2202493/148‑1 (MQ=255)
                                                                                 aTGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGatta                                                                 <  2:1250552/148‑1 (MQ=255)
                                                                                  tGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGattat                                                                <  1:2435032/148‑1 (MQ=255)
                                                                                    cTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTg                                                              <  2:2131395/148‑1 (MQ=255)
                                                                                       gcgcgcCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCa                                                           <  2:1588653/148‑1 (MQ=255)
                                                                                        cgcgcCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCaa                                                          <  2:3255051/148‑1 (MQ=255)
                                                                                                ttAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCACGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGa                                                  >  1:1650736/1‑148 (MQ=255)
                                                                                                  aGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGAcc                                                <  1:1427993/148‑1 (MQ=255)
                                                                                                        tGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGa                                          >  2:630034/1‑148 (MQ=255)
                                                                                                         gAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGaa                                         <  1:2553663/148‑1 (MQ=255)
                                                                                                         gAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGaa                                         <  1:2613848/148‑1 (MQ=255)
                                                                                                         gAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGaa                                         <  1:420788/148‑1 (MQ=255)
                                                                                                          aaGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGaaa                                        >  2:1687921/1‑148 (MQ=255)
                                                                                                          aaGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGaaa                                        >  1:543100/1‑148 (MQ=255)
                                                                                                            gTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATc                                      <  1:1469809/148‑1 (MQ=255)
                                                                                                             tAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCTACGCCATGACCCGTGGAAATcc                                     <  1:2174285/148‑1 (MQ=255)
                                                                                                              aaCGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCg                                    <  1:2605433/148‑1 (MQ=255)
                                                                                                               aCGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCgt                                   >  2:3271615/1‑148 (MQ=255)
                                                                                                                cGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCgtg                                  <  2:1282775/148‑1 (MQ=255)
                                                                                                                 ggggTCGGACGCCACTGTGTATCCGCCGATGAGCTCCTTCCTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCgtgt                                 <  2:1716597/148‑1 (MQ=255)
                                                                                                                 gggTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCgtgt                                 <  2:2321589/148‑1 (MQ=255)
                                                                                                                 gggTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCgtgt                                 <  2:2321949/148‑1 (MQ=255)
                                                                                                                   gTTCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCgtgtgt                               >  2:2686216/1‑148 (MQ=255)
                                                                                                                     tCGGACGCCAATGTGTATCCGCCGATGAGCTCCTTCCTGGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTgg                             <  2:1727415/148‑1 (MQ=255)
                                                                                                                      cGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGa                            <  2:50774/148‑1 (MQ=255)
                                                                                                                      cGGACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGa                            <  1:1872764/148‑1 (MQ=255)
                                                                                                                        gACGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAg                          >  1:1749060/1‑148 (MQ=255)
                                                                                                                         aCGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGc                         >  1:2135616/1‑148 (MQ=255)
                                                                                                                         aCGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGc                         <  2:1358398/148‑1 (MQ=255)
                                                                                                                         aCGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGc                         >  1:2061638/1‑148 (MQ=255)
                                                                                                                          cGCCAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCg                        <  2:207397/148‑1 (MQ=255)
                                                                                                                           gCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCgg                       <  2:289033/148‑1 (MQ=255)
                                                                                                                             cAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTa                     <  1:1243432/148‑1 (MQ=255)
                                                                                                                             cAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTa                     <  1:2970173/148‑1 (MQ=255)
                                                                                                                             cAATGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTa                     >  2:1842738/1‑148 (MQ=255)
                                                                                                                              aaTGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTAc                    <  1:2311835/148‑1 (MQ=255)
                                                                                                                               aTGTGTATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTACt                   >  2:1398622/1‑148 (MQ=255)
                                                                                                                                    tATCCGCCGATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGaaa              >  2:2901633/1‑148 (MQ=255)
                                                                                                                                       ccgccgATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGaaaaaa           <  1:1306067/148‑1 (MQ=255)
                                                                                                                                          ccgATGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGaaaaaaaat        <  1:1297473/148‑3 (MQ=255)
                                                                                                                                             aTGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATcc     <  2:1957322/148‑1 (MQ=255)
                                                                                                                                             aTGAGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATcc     >  1:603233/1‑148 (MQ=255)
                                                                                                                                                aGCTCCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCtt  <  2:604419/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GAGCGGCAGATTAAAAAAGGATAGTGACGTATGCGCATTCATATTTTAGGAATTTGTGGCACATTTATGGGCGGTCTGGCGATGCTGGCGCGCCAGTTAGGCCATGAAGTAACGGGTTCGGACGCCAATGTGTATCCGCCGATGAGCACCTTACTTGAGAAGCAAGGCATTGAGCTGATTCAGGGGTACGATGCCAGCCAGCTCGAGCCGCAGCCGGATCTGGTGATTATTGGCAACGCCATGACCCGTGGAAATCCGTGTGTGGAAGCGGTACTGGAAAAAAATATCCCTT  >  CP009273/4445572‑4445863

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