Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,056,803 A→T 8.0% M1K (ATG→AAG) † ldtA ← L,D‑transpeptidase linking Lpp to murein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,056,8030AT 8.0% 997.2 / 59.7 292M1K (ATG→AAG) †ldtAL,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base A (134/134);  new base T (7/16);  total (142/150)
Fisher's exact test for biased strand distribution p-value = 8.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-01

Only 100 of 312 total aligned reads displayed.

AGCGGCTGGGTATTGTGATCAGGAACAGTTACAGTAAACGACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCAAGTTCGCCCATGTTGTTGTATATCGCAC  >  CP009273/2056656‑2056951
                                                                                                                                                   |                                                                                                                                                    
          tagtgtgATCAGGAACAGTTACTGTAAACGACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCtt                                                                                                                                           <  2:3095252/146‑1 (MQ=255)
           aTTGTGATCAGTAACATTTACATAAAACGACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCttt                                                                                                                                          <  1:3751731/149‑1 (MQ=255)
            ttttgATCAGGAACATTTACTGTAATCGACTGCCCCAATAAACGGCTACCCTCTGGAGGTAATGGATATGTTACCTCAAGTCTAGTATGGCTGGCATAATGCAGAGCAAATTAGCAAAGAATATTTACACGACGCTTCATGTCCCTTTc                                                                                                                                         <  2:1269594/146‑1 (MQ=255)
            ttGTGATCAGGAACAGTTACAGTAAACGACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTc                                                                                                                                         >  1:1286374/1‑149 (MQ=255)
               tgATCCGGAACAGTTACTGTACACGACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTa                                                                                                                                      <  1:5205060/149‑1 (MQ=255)
                gATCAGGAACAGTTACAGTAAACGACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTAt                                                                                                                                     >  1:3708763/1‑149 (MQ=255)
                gATCAGGAACAGTTACAGTAAACGACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTAt                                                                                                                                     >  2:268803/1‑149 (MQ=255)
                 aTCAGGAACAGTTACAGTAAACGACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATg                                                                                                                                    <  1:408684/149‑1 (MQ=255)
                    aGGAACAGTTACAGTAAACGACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTcg                                                                                                                                 >  2:3524433/1‑149 (MQ=255)
                     ggAACAGTTACAGTAAACGACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTcgc                                                                                                                                >  1:616928/1‑149 (MQ=255)
                      gAACAGTTACAGTAAACGACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTGTCCTATGTcgcg                                                                                                                               >  2:4012098/1‑149 (MQ=255)
                       aaCAGTTACAGTAAACGACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGa                                                                                                                              <  2:4538118/149‑1 (MQ=255)
                       aaCAGTTACAGTAAACGACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGa                                                                                                                              <  2:334117/149‑1 (MQ=255)
                            ttACAGTAAACGACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCTAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCt                                                                                                                         <  2:4286779/149‑1 (MQ=255)
                            ttACAGTAAACGACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCt                                                                                                                         >  1:4488034/1‑149 (MQ=255)
                             aactgaaaaCGACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCTTCATGTCCCTTTCCTATGTCGCGAAAGCTa                                                                                                                        <  1:4792901/143‑1 (MQ=255)
                                aGTAAACGACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATcc                                                                                                                     <  2:731079/149‑1 (MQ=255)
                                  tAAACGACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCTTCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGt                                                                                                                   <  1:3248367/149‑1 (MQ=255)
                                    aaCGACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTa                                                                                                                 >  1:3453404/1‑149 (MQ=255)
                                     aCGACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCTTCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTaa                                                                                                                <  2:1320208/149‑1 (MQ=255)
                                      cGACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAg                                                                                                               <  2:4194536/149‑1 (MQ=255)
                                       gACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGATGCGTCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGt                                                                                                              >  1:760483/1‑149 (MQ=255)
                                        aCTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGTATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGta                                                                                                             <  2:4123878/149‑1 (MQ=255)
                                        aCTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGta                                                                                                             >  1:4623021/1‑149 (MQ=255)
                                        aCTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGta                                                                                                             <  2:4053685/149‑1 (MQ=255)
                                          tGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGtata                                                                                                           <  1:2632967/149‑1 (MQ=255)
                                          tGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGtata                                                                                                           <  1:4422997/149‑1 (MQ=255)
                                            ccccACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGc                                                                                                         <  2:3480230/149‑1 (MQ=255)
                                              ccACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCtt                                                                                                       <  1:5154964/149‑1 (MQ=255)
                                               cACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCttt                                                                                                      >  1:2881330/1‑149 (MQ=255)
                                                aCTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCtttt                                                                                                     <  2:2498590/149‑1 (MQ=255)
                                                 cTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTa                                                                                                    <  2:1237728/149‑1 (MQ=255)
                                                 cTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTa                                                                                                    <  1:2313640/149‑1 (MQ=255)
                                                   aaaCGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATc                                                                                                  <  1:2075198/149‑1 (MQ=255)
                                                   aaaCGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATc                                                                                                  <  1:3979702/149‑1 (MQ=255)
                                                    aaCGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCa                                                                                                 >  1:2018322/1‑149 (MQ=255)
                                                     aCGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAg                                                                                                >  1:5013885/1‑149 (MQ=255)
                                                           aCCCTCTGGAGGTAATGTATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACtttt                                                                                          <  1:1450712/149‑1 (MQ=255)
                                                           aCCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACtttt                                                                                          >  2:3747749/1‑149 (MQ=255)
                                                           aCCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACtttt                                                                                          >  1:1360313/1‑149 (MQ=255)
                                                            cccTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTc                                                                                         <  2:889615/149‑1 (MQ=255)
                                                            cccTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTc                                                                                         >  1:254103/1‑149 (MQ=255)
                                                            cccTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTc                                                                                         <  1:1866460/149‑1 (MQ=255)
                                                             ccTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCg                                                                                        <  2:2702474/149‑1 (MQ=255)
                                                               tctGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGtt                                                                                      <  1:2493810/149‑1 (MQ=255)
                                                                  ggAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGttttt                                                                                   >  2:1845217/1‑149 (MQ=255)
                                                                  ggAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGttttt                                                                                   >  2:4699131/1‑149 (MQ=255)
                                                                     ggTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAAc                                                                                <  1:4606912/149‑1 (MQ=255)
                                                                       tAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTg                                                                              <  1:3781273/149‑1 (MQ=255)
                                                                         aTGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGtt                                                                            >  2:213708/1‑149 (MQ=255)
                                                                          tGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTc                                                                           <  1:2400521/149‑1 (MQ=255)
                                                                              tAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCTTCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAAt                                                                       >  2:561435/1‑149 (MQ=255)
                                                                               aaGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATc                                                                      <  2:3404963/149‑1 (MQ=255)
                                                                                 gTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAg                                                                    >  2:4501476/1‑149 (MQ=255)
                                                                                   tACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGaa                                                                  >  1:3577476/1‑149 (MQ=255)
                                                                                    aCCGCCAGGCTAGTATGGCGGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAg                                                                 <  1:3222276/149‑1 (MQ=255)
                                                                                     ccGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCACGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGt                                                                <  2:1133225/149‑1 (MQ=255)
                                                                                      cGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCCAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTc                                                               >  2:1872450/1‑149 (MQ=255)
                                                                                           ggCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTAtt                                                          <  2:614905/149‑1 (MQ=255)
                                                                                           ggCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTAtt                                                          <  1:1518559/149‑1 (MQ=255)
                                                                                           ggCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTAtt                                                          <  2:367626/149‑1 (MQ=255)
                                                                                             cTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAATAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTcc                                                        <  1:4880113/149‑1 (MQ=255)
                                                                                             cTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTcc                                                        <  1:2704891/149‑1 (MQ=255)
                                                                                             cTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTcc                                                        <  1:101994/149‑1 (MQ=255)
                                                                                              tAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTccc                                                       >  2:4257304/1‑149 (MQ=255)
                                                                                              tAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAAGCAGAAGTCGTATTccc                                                       >  2:4257929/1‑149 (MQ=255)
                                                                                                     gctggcAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGa                                                <  2:3551210/149‑1 (MQ=255)
                                                                                                      ctggcAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGaa                                               <  2:2911992/149‑1 (MQ=255)
                                                                                                       tggcAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAAc                                              <  2:1656757/149‑1 (MQ=255)
                                                                                                       tggcAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAAc                                              <  2:1374699/149‑1 (MQ=255)
                                                                                                       tggcAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAAc                                              >  1:2819373/1‑149 (MQ=255)
                                                                                                            aaaaGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAatat                                         <  1:2772785/149‑1 (MQ=255)
                                                                                                            aaaaGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAatat                                         >  1:3136935/1‑149 (MQ=255)
                                                                                                               aGCAGAGCAAATGAGCAAAGAATATTTACACGACGCTTCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTAc                                      <  1:975733/149‑1 (MQ=255)
                                                                                                                   gagCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGgc                                  <  1:4582872/149‑1 (MQ=255)
                                                                                                                   gagCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGgc                                  >  1:4239746/1‑149 (MQ=255)
                                                                                                                     gCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGgcag                                >  2:4341669/1‑149 (MQ=255)
                                                                                                                        aaTGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCaa                             <  1:1734081/149‑1 (MQ=255)
                                                                                                                        aaTGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCACGGTAGAACAATATTACTGGCAGCaa                             >  2:4093998/1‑149 (MQ=255)
                                                                                                                           gAGCAAAGAATATTTACACGACTCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCAAGtt                          >  1:1916787/1‑149 (MQ=255)
                                                                                                                           gAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCAAGtt                          >  2:2255254/1‑149 (MQ=255)
                                                                                                                            aGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCAAGTTc                         >  1:3040739/1‑149 (MQ=255)
                                                                                                                            aGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCAAGTTc                         >  1:3657248/1‑149 (MQ=255)
                                                                                                                                   aaTATTTACACGACGCTTCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCAAGTTCGCCCAtg                  >  1:823938/1‑149 (MQ=255)
                                                                                                                                     tatTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCAAGTTCGCCCAtgtt                >  2:1650311/1‑149 (MQ=255)
                                                                                                                                     tatTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCAAGTTCGCCCAtgtt                >  2:1089410/1‑149 (MQ=255)
                                                                                                                                       tttACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCAAGTTCGCCCAtgttgt              >  1:3161903/1‑149 (MQ=255)
                                                                                                                                       tttACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCAAGTTCGCCCAtgttgt              <  2:3708763/149‑1 (MQ=255)
                                                                                                                                           cacGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCAAGTTCGCCCATGTTGTTGta          <  1:1581155/149‑1 (MQ=255)
                                                                                                                                            acgacgCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCAAGTTCGCCCCTGTTGTTGtat         >  2:4867445/1‑149 (MQ=255)
                                                                                                                                             cgacgCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCAAGTTCGCCCATGTTGTTGtata        <  1:2142307/149‑1 (MQ=255)
                                                                                                                                              gacgCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCAAGTTCGCCCATGTTGTTGTatat       >  2:4916226/1‑149 (MQ=255)
                                                                                                                                               acgCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCAAGTTCGCCCATGTTGTTGTATATc      >  1:4163270/1‑149 (MQ=255)
                                                                                                                                                cgCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCAAGTTCGCCCATGTTGTTGTATATCg     >  2:4484446/1‑149 (MQ=255)
                                                                                                                                                  catcatGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTTTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCAAGTTCGCCCATGTTGTTGTATATCGca   >  2:2336172/1‑149 (MQ=255)
                                                                                                                                                   atcatGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCATAAGTCGTATTCCCCGGTAGAACAATATTACTGTCAGCAAGTTCGCCCATGTTGTTGTATATCGcac  >  2:3879466/1‑149 (MQ=255)
                                                                                                                                                   atcatGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCAAGTTCGCCCATGTTGTTGTATATCGcac  >  2:3831939/1‑149 (MQ=255)
                                                                                                                                                   atcatGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCAAGTTCGCCCATGTTGTTGTATATCGcac  >  2:5084567/1‑149 (MQ=255)
                                                                                                                                                   atcatGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCAAGTTCGCCCATGTTGTTGTATATCGcac  >  1:309889/1‑149 (MQ=255)
                                                                                                                                                   atcatGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCAAGTTCGCCCATGTTGTTGTATATCGcac  >  2:2981625/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGCGGCTGGGTATTGTGATCAGGAACAGTTACAGTAAACGACTGCCCCACTAAACGGCTACCCTCTGGAGGTAATGGATAAGTTACCGCCAGGCTAGTATGGCTGGCAAAAAGCAGAGCAAATGAGCAAAGAATATTTACACGACGCATCATGTCCCTTTCCTATGTCGCGAAAGCTATCCGTTAAGTATAGCTTTTATCAGACTTTTCGTTTTTAACTGTTCAAATCAGAAGTCGTATTCCCCGGTAGAACAATATTACTGGCAGCAAGTTCGCCCATGTTGTTGTATATCGCAC  >  CP009273/2056656‑2056951

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: