Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 627,794 G→A 5.3% D53N (GAC→AAC)  ybdD → DUF466 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273627,7940GA 5.3% 1020.0 / 28.4 284D53N (GAC→AAC) ybdDDUF466 family protein
Reads supporting (aligned to +/- strand):  ref base G (132/137);  new base A (6/9);  total (138/146)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

Only 100 of 298 total aligned reads displayed.

TTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCGGACGCAGGG  >  CP009273/627646‑627942
                                                                                                                                                    |                                                                                                                                                    
ttCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATCTGCGGGTTAACCCTCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAgg                                                                                                                                                      >  2:4413801/1‑149 (MQ=255)
ttCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAgg                                                                                                                                                      >  1:2845604/1‑149 (MQ=255)
 tCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGa                                                                                                                                                     <  1:845475/149‑1 (MQ=255)
 tCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGa                                                                                                                                                     <  1:707370/149‑1 (MQ=255)
    cTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGAcgc                                                                                                                                                  >  2:2236988/1‑149 (MQ=255)
         aaaaGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCt                                                                                                                                             >  1:987586/1‑149 (MQ=255)
         aaaaGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCt                                                                                                                                             >  1:3879291/1‑149 (MQ=255)
            aGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACg                                                                                                                                          >  2:657588/1‑149 (MQ=255)
              ccGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGAACGCGCGCTACGGt                                                                                                                                        <  2:3821598/149‑1 (MQ=255)
               cGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTg                                                                                                                                       >  2:2942189/1‑149 (MQ=255)
                 gAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGa                                                                                                                                     >  1:1101655/1‑149 (MQ=255)
                 gAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGa                                                                                                                                     >  2:3672285/1‑149 (MQ=255)
                     atatTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAg                                                                                                                                 <  2:838994/149‑1 (MQ=255)
                      tatTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAgg                                                                                                                                >  2:878828/1‑149 (MQ=255)
                      tatTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAgg                                                                                                                                >  1:4512916/1‑149 (MQ=255)
                      tatTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAgg                                                                                                                                <  1:3536746/149‑1 (MQ=255)
                        tttAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAggcg                                                                                                                              <  1:1890630/149‑1 (MQ=255)
                           aGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGAGGCAGGACGCGCGCTACGGTGGAAAAGgcggcg                                                                                                                           >  2:1942302/1‑149 (MQ=255)
                            ggTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGgcgc                                                                                                                          >  2:4582715/1‑149 (MQ=255)
                            ggTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGgcgc                                                                                                                          >  1:339428/1‑149 (MQ=255)
                              tCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGgcgcgc                                                                                                                        <  2:547637/149‑1 (MQ=255)
                              tCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGgcgcgc                                                                                                                        <  2:4396490/149‑1 (MQ=255)
                                   cggcgAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCgctgc                                                                                                                   <  2:3660757/149‑1 (MQ=255)
                                   cggcgAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGAACGCGCGCTACGGTGGAAAAGGCGGCGCGCgctgc                                                                                                                   <  2:1837827/149‑1 (MQ=255)
                                         aGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTaa                                                                                                             >  2:4038929/1‑149 (MQ=255)
                                          gCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTaaa                                                                                                            <  2:1404020/149‑1 (MQ=255)
                                          gCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTaaa                                                                                                            <  2:2769485/149‑1 (MQ=255)
                                            tgatgaTTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAAtt                                                                                                          >  2:3406889/1‑149 (MQ=255)
                                            tgatgaTTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAAtt                                                                                                          <  2:3049298/149‑1 (MQ=255)
                                             gatgatTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAAttt                                                                                                         <  1:2159120/149‑1 (MQ=255)
                                               tgatTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGt                                                                                                       >  1:3160773/1‑149 (MQ=255)
                                                 atTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTAc                                                                                                     <  2:3107478/149‑1 (MQ=255)
                                                  ttGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACg                                                                                                    >  1:424418/1‑149 (MQ=255)
                                                       aTGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAgg                                                                                               <  2:2035057/149‑1 (MQ=255)
                                                            tGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTa                                                                                          <  2:1724780/149‑1 (MQ=255)
                                                            tGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTa                                                                                          >  2:1017156/1‑149 (MQ=255)
                                                            tGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTa                                                                                          <  2:2077332/149‑1 (MQ=255)
                                                             gATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGAACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTaa                                                                                         >  1:3954613/1‑149 (MQ=255)
                                                             gATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACAAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTaa                                                                                         >  2:2469795/1‑149 (MQ=255)
                                                              aTTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTaaa                                                                                        <  1:1346568/149‑1 (MQ=255)
                                                              aTTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTaaa                                                                                        <  2:1486154/149‑1 (MQ=255)
                                                                 aCGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGa                                                                                     <  1:2489147/149‑1 (MQ=255)
                                                                 aCGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGa                                                                                     <  2:2572563/149‑1 (MQ=255)
                                                                   gACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGAtt                                                                                   >  1:555773/1‑149 (MQ=255)
                                                                       aCTATGTCGAACATATGCGTGTTAACCATCCCGATCAAACGCCGATGACCTAGGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCAc                                                                               <  1:3751404/149‑1 (MQ=255)
                                                                          aTGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGTGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACttt                                                                            >  1:434034/1‑149 (MQ=255)
                                                                          aTGTCGAACATATGCGGGTTAACCAGCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACttt                                                                            <  1:183716/149‑1 (MQ=255)
                                                                           tGTCGAACATTTGCGGGTTAACCAGCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTAAGGCTTTAAACGATTCCACtttt                                                                           <  1:338193/149‑1 (MQ=255)
                                                                             tCGAACATATGCGGGTTAACAATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACtttttt                                                                         >  1:1937952/1‑149 (MQ=255)
                                                                              cGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTCAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTg                                                                        >  2:2132498/1‑149 (MQ=255)
                                                                              cGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTg                                                                        >  2:3313193/1‑149 (MQ=255)
                                                                                aaCATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGaa                                                                      >  2:155833/1‑149 (MQ=255)
                                                                                   atatGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGc                                                                   >  1:3746095/1‑149 (MQ=255)
                                                                                   atatGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGc                                                                   <  2:4203078/149‑1 (MQ=255)
                                                                                   atatGAGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTTAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTTTACGTCAGGCTTTAAACGATTCAACTTTTTTGAACGc                                                                   >  2:3466602/1‑149 (MQ=255)
                                                                                    tatGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCt                                                                  >  1:1239097/1‑149 (MQ=255)
                                                                                      tGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGc                                                                >  2:1224939/1‑149 (MQ=255)
                                                                                      tGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGc                                                                >  1:2109508/1‑149 (MQ=255)
                                                                                         gggTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACg                                                             >  1:1019824/1‑149 (MQ=255)
                                                                                           gTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCa                                                           >  1:289574/1‑149 (MQ=255)
                                                                                           gTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCa                                                           >  1:2630007/1‑149 (MQ=255)
                                                                                             tAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGAACGCGCGCTACGGTGGTAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAAc                                                         <  2:2887649/149‑1 (MQ=255)
                                                                                               aCCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAATCGATTCCACTTGTTTGAACGCTGCTCGCAACGt                                                       >  2:2830845/1‑149 (MQ=255)
                                                                                               aCCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGt                                                       >  1:255624/1‑149 (MQ=255)
                                                                                                ccATCCCGATCAAACGCCGATGACCTACGACGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGATGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTa                                                      <  2:3548718/149‑1 (MQ=255)
                                                                                                ccATCCCGATCAAACGCCGATGACATACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCCGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTa                                                      <  1:2009179/149‑1 (MQ=255)
                                                                                                 cATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTAt                                                     >  2:1169072/1‑149 (MQ=255)
                                                                                                  aTCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATc                                                    >  2:4243943/1‑149 (MQ=255)
                                                                                                   tCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCt                                                   >  1:1718214/1‑149 (MQ=255)
                                                                                                   tCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCt                                                   >  1:501494/1‑149 (MQ=255)
                                                                                                   tCCCGATCAAACGCCGATGACCTACGAAGAGTTGTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGATAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCt                                                   <  2:3259731/149‑1 (MQ=255)
                                                                                                      cGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGAACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGgt                                                <  2:4626041/149‑1 (MQ=255)
                                                                                                       gATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGgtg                                               >  2:3841954/1‑149 (MQ=255)
                                                                                                        aTCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGgtgt                                              >  1:3324136/1‑149 (MQ=255)
                                                                                                        aTCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGgtgt                                              <  2:3790835/149‑1 (MQ=255)
                                                                                                         tCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTc                                             <  2:638848/149‑1 (MQ=255)
                                                                                                          cAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCa                                            >  2:434481/1‑149 (MQ=255)
                                                                                                           acaCGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAg                                           <  1:1601857/147‑1 (MQ=255)
                                                                                                                    tGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTgg                                  <  1:3101643/149‑1 (MQ=255)
                                                                                                                    tGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTgg                                  >  2:3209641/1‑149 (MQ=255)
                                                                                                                     gACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGc                                 >  1:116456/1‑149 (MQ=255)
                                                                                                                      aCCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCa                                >  2:450211/1‑149 (MQ=255)
                                                                                                                       ccTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAg                               >  2:2370706/1‑149 (MQ=255)
                                                                                                                       ccTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAg                               >  1:4672204/1‑149 (MQ=255)
                                                                                                                       ccTACGAAGAGTTTTTCCGTGAGCGGCAGAACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAg                               <  2:1599736/149‑1 (MQ=255)
                                                                                                                        cTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAgg                              <  2:1314854/149‑1 (MQ=255)
                                                                                                                        cTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAgg                              <  2:3754866/149‑1 (MQ=255)
                                                                                                                          aCGAAGAGTTTTTCCGTGAGCGGCAGAACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTa                            >  1:2258657/1‑149 (MQ=255)
                                                                                                                          aCGAAGAGTTTTTCCGTGAGCGGCAGAACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTa                            >  1:3202033/1‑149 (MQ=255)
                                                                                                                            gAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTaat                          >  1:2674788/1‑149 (MQ=255)
                                                                                                                                  tttttCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATgg                    >  2:1364515/1‑149 (MQ=255)
                                                                                                                                   ttttCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCATGACTGGCAGGTAATGAATGGa                   <  2:1488319/149‑1 (MQ=255)
                                                                                                                                    tttCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGAc                  >  2:606022/1‑149 (MQ=255)
                                                                                                                                    tttCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGAc                  >  2:732843/1‑149 (MQ=255)
                                                                                                                                         gTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCAc             <  1:4241990/149‑1 (MQ=255)
                                                                                                                                              cGGCAGGACGCGCGCTACGGTGTAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCGGAc        <  1:1456891/149‑1 (MQ=255)
                                                                                                                                                 cAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCGGACGCa     <  2:1142706/149‑1 (MQ=255)
                                                                                                                                                  aGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAACTTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCGGACGCAg    <  2:3010265/149‑1 (MQ=255)
                                                                                                                                                    gACGCGCGCTACGGTGGATAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCGGACGCAggg  <  2:3369837/149‑1 (MQ=255)
                                                                                                                                                    gACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCGGACGCAggg  >  1:4049500/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCGGACGCAGGG  >  CP009273/627646‑627942

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: