Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 91,765 C→A 6.6% P210Q (CCG→CAG)  murF → UDP‑N‑acetylmuramoyl‑tripeptide:D‑alanyl‑D‑ alanine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927391,7650CA 6.6% 974.2 / 26.3 281P210Q (CCG→CAG) murFUDP‑N‑acetylmuramoyl‑tripeptide:D‑alanyl‑D‑ alanine ligase
Reads supporting (aligned to +/- strand):  ref base C (132/127);  new base A (6/14);  total (138/143)
Fisher's exact test for biased strand distribution p-value = 1.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 287 total aligned reads displayed.

GCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCTCGCACGG  >  CP009273/91617‑91912
                                                                                                                                                    |                                                                                                                                                   
gcgAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGcc                                                                                                                                                     <  1:266389/149‑1 (MQ=255)
gcgAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGcc                                                                                                                                                     >  2:2324493/1‑149 (MQ=255)
  gAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCgg                                                                                                                                                   <  1:4164980/149‑1 (MQ=255)
   aaCCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGa                                                                                                                                                  >  1:2497322/1‑149 (MQ=255)
    aCCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGaa                                                                                                                                                 <  2:3886148/149‑1 (MQ=255)
       aTCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAAc                                                                                                                                              >  1:4370135/1‑149 (MQ=255)
       aTCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAAc                                                                                                                                              >  2:3706604/1‑149 (MQ=255)
          aGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGt                                                                                                                                           >  1:3801259/1‑149 (MQ=255)
            ggCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTAt                                                                                                                                         <  1:1035137/149‑1 (MQ=255)
               gAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGc                                                                                                                                      <  2:668988/149‑1 (MQ=255)
                aaaTAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGcc                                                                                                                                     <  2:3078271/149‑1 (MQ=255)
                  aTAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCAt                                                                                                                                   >  1:1810622/1‑149 (MQ=255)
                  aTAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCAt                                                                                                                                   >  2:2781213/1‑149 (MQ=255)
                  aTAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCAGGAAAACGGTATCGCCAt                                                                                                                                   <  2:1914128/149‑1 (MQ=255)
                     gCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTAt                                                                                                                                <  2:3131935/149‑1 (MQ=255)
                      ccTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATg                                                                                                                               >  2:4329036/1‑149 (MQ=255)
                       cTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGa                                                                                                                              <  1:1860774/149‑1 (MQ=255)
                       cTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGa                                                                                                                              >  2:4547141/1‑149 (MQ=255)
                       cTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCAGGCAAACGGTATCGCCATTATGa                                                                                                                              >  2:1028225/1‑149 (MQ=255)
                         ggACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAAc                                                                                                                            >  2:1568771/1‑149 (MQ=255)
                         ggACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAAc                                                                                                                            >  1:3466727/1‑149 (MQ=255)
                          gACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCAGGAAAACGGTATCGCCATTATGAACg                                                                                                                           <  2:704294/149‑1 (MQ=255)
                                gAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGaca                                                                                                                     >  1:2862900/1‑149 (MQ=255)
                                 aGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGacaa                                                                                                                    >  1:715011/1‑149 (MQ=255)
                                  gTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGacaac                                                                                                                   <  1:2474968/149‑1 (MQ=255)
                                       aCTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGa                                                                                                              <  1:1627178/149‑1 (MQ=255)
                                         tCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACt                                                                                                            >  2:2534487/1‑149 (MQ=255)
                                         tCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACt                                                                                                            <  2:2294317/149‑1 (MQ=255)
                                           gCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGATGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTgg                                                                                                          <  2:2433947/149‑1 (MQ=255)
                                            cccGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCTCATCTGGATGTTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGc                                                                                                         <  2:2871133/149‑1 (MQ=255)
                                             ccGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCt                                                                                                        >  1:2742084/1‑149 (MQ=255)
                                               ggAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGa                                                                                                      >  1:2936227/1‑149 (MQ=255)
                                                 atgCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGATGGTTTTGGCTCGCTTGGGGGTGTCGCGAAAGCTAAAGTTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCAATTATGAACGCCGACAACAACGACTGGCTGAAc                                                                                                    <  1:2648563/147‑1 (MQ=255)
                                                 aaGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAAc                                                                                                    >  1:3420526/1‑149 (MQ=255)
                                                 aaGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAAc                                                                                                    >  1:3964960/1‑149 (MQ=255)
                                                   gCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTg                                                                                                  >  1:1252542/1‑149 (MQ=255)
                                                    cTGCGCTGGTAAACAACCTGGCAGGGGCGCATGTGGATGGTTTTGGCTCGCTTGCGGGGGTCGCGAAAGCGAAAGGTTAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTgg                                                                                                 <  1:153031/149‑1 (MQ=255)
                                                            gTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGt                                                                                         <  2:3374134/149‑1 (MQ=255)
                                                               aacaacCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAAt                                                                                      >  2:3744038/1‑149 (MQ=255)
                                                                acaacCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAAtt                                                                                     >  2:1033459/1‑149 (MQ=255)
                                                                  aacCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTgg                                                                                   <  2:3542898/149‑1 (MQ=255)
                                                                  aacCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTgg                                                                                   >  1:1730865/1‑149 (MQ=255)
                                                                   acCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGc                                                                                  >  2:4383968/1‑149 (MQ=255)
                                                                    ccTGGCCGCGGCGATTCTGGATGGTTTTGGCTCGCTTGCGGGTGTCGCTAAAGCTAATGGTGAAATATTTTGCGGCCTGCCGGAAAACGGTATCGCAATTATCAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCt                                                                                 <  1:4114486/149‑1 (MQ=255)
                                                                      tGGCAGCGGCTCATCTGGTTGGTTTTGGCTCGGTGGGGGGTGTCGCGAAAGCTAATGGTGAAATCTTTAGCGGCCTGCGGGAAAACGGTAGCGCAATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCa                                                                               <  1:4453398/149‑1 (MQ=255)
                                                                      tGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAATGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCa                                                                               <  1:1471588/149‑1 (MQ=255)
                                                                       ggCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCAc                                                                              >  1:4004282/1‑149 (MQ=255)
                                                                        gCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACg                                                                             >  2:1434307/1‑149 (MQ=255)
                                                                            cggcgCATCGGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCaaa                                                                         <  2:4243397/149‑1 (MQ=255)
                                                                             ggcgCATCTGGAAGGTTTGGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGCCTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAg                                                                        <  2:390916/149‑1 (MQ=255)
                                                                                  aTCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTgg                                                                   <  1:3984721/149‑1 (MQ=255)
                                                                                  aTCGGGATGTTTTTGGCTCGCTTGGGGGTGTCGCTAACGCTAATGGTTAAATATTTAGCGGCCTGCCTCAAAACGGTATCGCAATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTgg                                                                   <  1:4059862/149‑1 (MQ=255)
                                                                                   tCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGc                                                                  >  1:2601044/1‑149 (MQ=255)
                                                                                   tCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGc                                                                  >  2:1278922/1‑149 (MQ=255)
                                                                                    cTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCAGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCACAGTGTGGCg                                                                 >  1:308262/1‑149 (MQ=255)
                                                                                       taaGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCAGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGttt                                                              <  1:2357419/148‑1 (MQ=255)
                                                                                        aaGTTTTTGGCTCGCTGGGGGGTGTCGCGAACGCGAATGGTGAAATGTTTAGCGGCCTGCCGGAAACCGGTATCGCCTTTTTGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTtc                                                             <  2:839640/149‑1 (MQ=255)
                                                                                          gtttttGGCTCGTGTGGGGGGGTCGCGAACGCTAATGTTGAAATCTTTAGCGGCCTGCCGTAAAACGGTATCGCAATTTTGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTtctc                                                           <  1:1284653/147‑1 (MQ=255)
                                                                                            ttttGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCAc                                                         >  1:3385498/1‑149 (MQ=255)
                                                                                            tttgGGCTCGCTTGGGGGTGTCGCGAAAGCGAATGGTTAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCAATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCAc                                                         <  1:3692942/149‑1 (MQ=255)
                                                                                                ggCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCaa                                                     >  2:1143036/1‑149 (MQ=255)
                                                                                                ggCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCaa                                                     <  2:3258360/149‑1 (MQ=255)
                                                                                                       ttGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCCTTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATgccgcc                                              >  2:4463342/1‑149 (MQ=255)
                                                                                                          cGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAAc                                           <  2:2936227/149‑1 (MQ=255)
                                                                                                          cGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAAc                                           <  1:3905512/149‑1 (MQ=255)
                                                                                                          cGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAAc                                           <  1:4272397/149‑1 (MQ=255)
                                                                                                            ggTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAg                                         <  1:1750366/149‑1 (MQ=255)
                                                                                                            ggTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAg                                         <  2:172886/149‑1 (MQ=255)
                                                                                                              tgtCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCg                                       <  2:522905/149‑1 (MQ=255)
                                                                                                                tCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGAt                                     >  1:2037730/1‑149 (MQ=255)
                                                                                                                  gcgaaagcgaaaGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGAttt                                   <  2:2169231/149‑1 (MQ=255)
                                                                                                                    gaaagcgaaagGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCa                                 >  2:894816/1‑149 (MQ=255)
                                                                                                                      aagcgaaagGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGACCTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCAcc                               >  2:3766629/1‑149 (MQ=255)
                                                                                                                      aagcgaaagGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCAcc                               >  1:4549775/1‑149 (MQ=255)
                                                                                                                        gcgaaagGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGc                             >  1:2208381/1‑149 (MQ=255)
                                                                                                                              ggTGAAATCTTTAGCGGCCTGCAGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCaa                       <  2:536255/149‑1 (MQ=255)
                                                                                                                               gTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAat                      >  2:159697/1‑149 (MQ=255)
                                                                                                                               gTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAat                      >  1:2361579/1‑149 (MQ=255)
                                                                                                                               gTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAat                      <  2:3613858/149‑1 (MQ=255)
                                                                                                                                tGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAata                     <  1:324319/149‑1 (MQ=255)
                                                                                                                                tGAAATATTTAGCGGCCTGCCGGAAAACGGTATCGCCTTTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAata                     <  1:521750/149‑1 (MQ=255)
                                                                                                                                 gAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAatat                    >  1:1672659/1‑149 (MQ=255)
                                                                                                                                  aaaTCTTTAGCGGCCTGCAGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATc                   <  2:4733446/149‑1 (MQ=255)
                                                                                                                                  aaaTCTTTAGCGGCCGGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATc                   <  1:2883833/149‑1 (MQ=255)
                                                                                                                                   aaTCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATcc                  <  2:3672461/149‑1 (MQ=255)
                                                                                                                                     tCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCAt                >  2:2044466/1‑149 (MQ=255)
                                                                                                                                     tCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCAt                >  2:630179/1‑149 (MQ=255)
                                                                                                                                       tttAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCAtgt              >  2:1863359/1‑149 (MQ=255)
                                                                                                                                        ttAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCAtgtg             <  2:273249/149‑1 (MQ=255)
                                                                                                                                          aGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGAc           >  1:4626191/1‑149 (MQ=255)
                                                                                                                                          aGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGAc           >  1:1250472/1‑149 (MQ=255)
                                                                                                                                           gCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGAcc          >  2:3771085/1‑149 (MQ=255)
                                                                                                                                             ggCCTGCCGTAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCTc        <  1:2006807/149‑1 (MQ=255)
                                                                                                                                             ggCCTGCCGGAAAACGGTATCGCCATTATTAACGCCGACAACAACGACTGGCTGAACTGGGAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCTc        <  1:2827407/149‑1 (MQ=255)
                                                                                                                                             ggCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCTc        >  2:303818/1‑149 (MQ=255)
                                                                                                                                             ggCCTGCCGGAAAACGGTATCGCAATTATTAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCTc        <  1:1056853/149‑1 (MQ=255)
                                                                                                                                              gCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCTCg       >  1:922601/1‑149 (MQ=255)
                                                                                                                                                cTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCTCGCa     <  1:4669830/149‑1 (MQ=255)
                                                                                                                                                cTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCTCGCa     <  1:2272567/149‑1 (MQ=255)
                                                                                                                                                   ccGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCTCGCACgg  <  1:649579/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCGAACCATCAGGGCGAAATAGCCTGGACTGTGAGTCTGACTCGCCCGGAAGCTGCGCTGGTCAACAACCTGGCAGCGGCGCATCTGGAAGGTTTTGGCTCGCTTGCGGGTGTCGCGAAAGCGAAAGGTGAAATCTTTAGCGGCCTGCCGGAAAACGGTATCGCCATTATGAACGCCGACAACAACGACTGGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCTCGCACGG  >  CP009273/91617‑91912

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