Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 627,794 G→A 5.6% 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.6% 1230.6 / 39.0 337D53N (GAC→AAC) ybdDDUF466 family protein
Reads supporting (aligned to +/- strand):  ref base G (139/179);  new base A (10/9);  total (149/188)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 339 total aligned reads displayed.

TTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCGGACGCAG  >  CP009273/627646‑627940
                                                                                                                                                    |                                                                                                                                                  
  cACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGAc                                                                                                                                                  <  2:2746192/149‑1 (MQ=255)
      ggCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACACCTATGTCGAACATATGCGGGTTAACCATCCCGATCACACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGAcgcgc                                                                                                                                              <  2:4915740/149‑1 (MQ=255)
        cAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACgcgcgc                                                                                                                                            <  1:934507/149‑1 (MQ=255)
        cAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACgcgcgc                                                                                                                                            <  2:2332568/149‑1 (MQ=255)
           aaGCCGGAAAATATTTAGGTCAGGCGGCGTAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTAc                                                                                                                                         >  2:1518406/1‑149 (MQ=255)
           aaGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGAACGCGCGCTAc                                                                                                                                         <  2:2725771/149‑1 (MQ=255)
           aaGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGATTTTCAGTGAGCGGCAGGACGCGCGCTAc                                                                                                                                         >  1:907817/1‑149 (MQ=255)
            aGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACg                                                                                                                                        >  1:2106456/1‑149 (MQ=255)
             gCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACgg                                                                                                                                       >  1:765311/1‑149 (MQ=255)
               cGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTg                                                                                                                                     <  1:3671178/149‑1 (MQ=255)
               cGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTg                                                                                                                                     <  1:1922176/149‑1 (MQ=255)
                 gAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGa                                                                                                                                   >  1:3068977/1‑149 (MQ=255)
                    aaTATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGTGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGaaaa                                                                                                                                >  1:88600/1‑149 (MQ=255)
                     atatTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAg                                                                                                                               <  2:1141579/149‑1 (MQ=255)
                      tatTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAgg                                                                                                                              >  2:4082335/1‑149 (MQ=255)
                       atTGAGGTCCGGCGGCGACGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAggc                                                                                                                             <  1:3801791/149‑1 (MQ=255)
                          tAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAggcggc                                                                                                                          <  2:4994379/149‑1 (MQ=255)
                              tCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAACAGGCGgcgcgc                                                                                                                      >  2:3056900/1‑149 (MQ=255)
                                 ggcggcGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCgct                                                                                                                   >  2:1201457/1‑149 (MQ=255)
                                    ggcgAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCgctgct                                                                                                                >  2:773594/1‑149 (MQ=255)
                                      cgAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCgctgctgc                                                                                                              <  2:535796/149‑1 (MQ=255)
                                      cgAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCgctgctgc                                                                                                              <  1:1096357/149‑1 (MQ=255)
                                        aaGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTa                                                                                                            >  2:2273138/1‑149 (MQ=255)
                                            tgatgaTTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTCAAtt                                                                                                        >  2:1834886/1‑149 (MQ=255)
                                            tgatgaTTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAAtt                                                                                                        >  1:284903/1‑149 (MQ=255)
                                            tgatgaTTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAAtt                                                                                                        >  2:4366338/1‑149 (MQ=255)
                                               tgatTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGt                                                                                                     <  1:3271584/149‑1 (MQ=255)
                                               tgatTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGt                                                                                                     <  2:467178/149‑1 (MQ=255)
                                                gatTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTa                                                                                                    <  2:3699153/149‑1 (MQ=255)
                                                 atTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTAc                                                                                                   >  1:2963336/1‑149 (MQ=255)
                                                 atTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTAc                                                                                                   <  2:2917812/149‑1 (MQ=255)
                                                  ttGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACg                                                                                                  <  1:1248782/149‑1 (MQ=255)
                                                  ttGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACg                                                                                                  >  2:4753211/1‑149 (MQ=255)
                                                   tGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGt                                                                                                 <  2:1295034/149‑1 (MQ=255)
                                                    ggTATGCCTGATTACGACACCTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTc                                                                                                <  2:1417296/149‑1 (MQ=255)
                                                     gTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCa                                                                                               >  1:3041404/1‑149 (MQ=255)
                                                        tGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGc                                                                                            <  1:202954/149‑1 (MQ=255)
                                                        tGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGAACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGc                                                                                            <  1:4406488/149‑1 (MQ=255)
                                                         gCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCt                                                                                           <  1:2702247/149‑1 (MQ=255)
                                                          ccTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCtt                                                                                          <  1:4652931/149‑1 (MQ=255)
                                                          ccTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCtt                                                                                          <  1:2683158/149‑1 (MQ=255)
                                                            tGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTa                                                                                        <  2:3222148/149‑1 (MQ=255)
                                                             gATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTaa                                                                                       >  2:222420/1‑149 (MQ=255)
                                                             gATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTaa                                                                                       >  1:3900650/1‑149 (MQ=255)
                                                              aTTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTaaa                                                                                      <  2:3056622/149‑1 (MQ=255)
                                                              aTTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTaaa                                                                                      <  2:4672242/149‑1 (MQ=255)
                                                               ttACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAAc                                                                                     <  1:3983514/149‑1 (MQ=255)
                                                                 aCGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGa                                                                                   >  2:177791/1‑149 (MQ=255)
                                                                   gACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGAACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGAtt                                                                                 >  2:853634/1‑149 (MQ=255)
                                                                    aCAACTTTGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTc                                                                                <  1:5113557/149‑1 (MQ=255)
                                                                    aCAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTc                                                                                >  2:5127292/1‑149 (MQ=255)
                                                                     cAACTATGTCTAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGAACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGATAAATTTGTACGTCAGGCTTTAAACGATTcc                                                                               <  1:701747/149‑1 (MQ=255)
                                                                       aCTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCAc                                                                             >  1:462106/1‑149 (MQ=255)
                                                                       aCTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCAc                                                                             >  2:4047327/1‑149 (MQ=255)
                                                                       aCTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCAc                                                                             >  1:4932785/1‑149 (MQ=255)
                                                                       aCTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCAc                                                                             >  2:1910503/1‑149 (MQ=255)
                                                                         tATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACtt                                                                           >  2:459988/1‑149 (MQ=255)
                                                                         tATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACtt                                                                           >  2:2081661/1‑149 (MQ=255)
                                                                          ctGTCGAACATTGGGGGTTTAACCAGCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAATCGATTCCACttt                                                                          <  1:1777514/148‑1 (MQ=255)
                                                                           tGTCTAAAATAGGGGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGATTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACtttt                                                                         <  1:3115492/149‑1 (MQ=255)
                                                                            gTCTAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACttttt                                                                        <  2:1397735/149‑1 (MQ=255)
                                                                            gTCGAACATATGCGGGTTAACCATCCCTATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACttttt                                                                        <  1:2913789/149‑1 (MQ=255)
                                                                             tCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACtttttt                                                                       <  2:3024292/149‑1 (MQ=255)
                                                                                aaCATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGaa                                                                    >  1:4763654/1‑149 (MQ=255)
                                                                                  cATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACg                                                                  >  2:4441363/1‑149 (MQ=255)
                                                                                     atGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTg                                                               >  2:703806/1‑149 (MQ=255)
                                                                                         gggTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACg                                                           <  1:1460062/149‑1 (MQ=255)
                                                                                             tAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAAc                                                       >  2:3879260/1‑149 (MQ=255)
                                                                                             tAACCAGCCCCATCAACCGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGTAAAAGGCGGCGCGCGCTGCTGCTAAATTGGTACGTCTGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAAc                                                       <  1:3515314/149‑1 (MQ=255)
                                                                                               aCCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTGTTTTGAACGCTGCACGCACCGt                                                     >  2:2459891/1‑149 (MQ=255)
                                                                                                  aTCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATc                                                  <  2:3411546/149‑1 (MQ=255)
                                                                                                         tCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGATTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTc                                           <  1:1011416/149‑1 (MQ=255)
                                                                                                          cAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCa                                          >  1:840598/1‑149 (MQ=255)
                                                                                                           aaaCGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAg                                         >  1:3147268/1‑149 (MQ=255)
                                                                                                           aaaCGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAg                                         >  1:1946473/1‑149 (MQ=255)
                                                                                                             aCGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCg                                       >  1:3768798/1‑149 (MQ=255)
                                                                                                               gCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTACACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTg                                     >  2:3962549/1‑149 (MQ=255)
                                                                                                                  gATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACt                                  >  1:3414598/1‑149 (MQ=255)
                                                                                                                   aTGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTg                                 >  2:4844733/1‑149 (MQ=255)
                                                                                                                    tGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTgg                                <  2:1652574/149‑1 (MQ=255)
                                                                                                                     gACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGc                               <  1:1637549/149‑1 (MQ=255)
                                                                                                                       ccTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAATGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAg                             <  1:2355449/149‑1 (MQ=255)
                                                                                                                       ccTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAg                             <  1:524269/149‑1 (MQ=255)
                                                                                                                         tACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGt                           <  1:4101410/149‑1 (MQ=255)
                                                                                                                         tACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGt                           <  1:1504965/149‑1 (MQ=255)
                                                                                                                           cGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTaa                         <  1:828454/149‑1 (MQ=255)
                                                                                                                            gAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTaat                        <  1:4244281/149‑1 (MQ=255)
                                                                                                                               gagTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTaatgaa                     >  2:2163211/1‑149 (MQ=255)
                                                                                                                               gagTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTGGTACGTCCGGCTTTACACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTaatgaa                     <  1:41525/149‑1 (MQ=255)
                                                                                                                                agTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTaatgaat                    <  2:2670345/149‑1 (MQ=255)
                                                                                                                                 gTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTaatgaatg                   <  1:1818149/149‑1 (MQ=255)
                                                                                                                                    tttCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGAc                <  2:4552627/149‑1 (MQ=255)
                                                                                                                                    tttCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGAc                >  1:4119948/1‑149 (MQ=255)
                                                                                                                                        cGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCa            <  2:253612/149‑1 (MQ=255)
                                                                                                                                            aGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCgg        <  1:2625136/149‑1 (MQ=255)
                                                                                                                                                 cAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCATCGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCGGACGCa   <  2:4193877/149‑1 (MQ=255)
                                                                                                                                                  aGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCGGACGCAg  >  2:4975277/1‑149 (MQ=255)
                                                                                                                                                  aGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCGGACGCAg  <  2:2068061/149‑1 (MQ=255)
                                                                                                                                                  aGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCGGACGCAg  <  2:2268022/149‑1 (MQ=255)
                                                                                                                                                  aGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCGGACGCAg  >  2:4039597/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTCACTGGCAAAAGCCGGAAAATATTTAGGTCAGGCGGCGAAGCTGATGATTGGTATGCCTGATTACGACAACTATGTCGAACATATGCGGGTTAACCATCCCGATCAAACGCCGATGACCTACGAAGAGTTTTTCCGTGAGCGGCAGGACGCGCGCTACGGTGGAAAAGGCGGCGCGCGCTGCTGCTAAATTTGTACGTCAGGCTTTAAACGATTCCACTTTTTTGAACGCTGCACGCAACGTATCTGGTGTCAGCGTGACTGGCAGGTAATGAATGGACTCCACCGGACGCAG  >  CP009273/627646‑627940

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: