Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,924,898 T→A 13.9% D36V (GAT→GTT)  yebG ← DNA damage‑inducible protein regulated by LexA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,924,8980TA13.9% 933.9 / 121.0 298D36V (GAT→GTT) yebGDNA damage‑inducible protein regulated by LexA
Reads supporting (aligned to +/- strand):  ref base T (88/168);  new base A (21/21);  total (109/189)
Fisher's exact test for biased strand distribution p-value = 5.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 299 total aligned reads displayed.

AGCACCGACGACCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGTATATT  >  CP009273/1924750‑1925045
                                                                                                                                                    |                                                                                                                                                   
     cgacgaCCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATccgcc                                                                                                                                                <  1:275492/149‑1 (MQ=255)
     cgacgaCCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATccgcc                                                                                                                                                <  2:3101345/149‑1 (MQ=255)
       acgacCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGt                                                                                                                                              >  2:432921/1‑149 (MQ=255)
       acgacCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGt                                                                                                                                              <  1:1685022/149‑1 (MQ=255)
        cgacCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTa                                                                                                                                             >  2:3214153/1‑149 (MQ=255)
          acCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATc                                                                                                                                           <  2:188650/149‑1 (MQ=255)
          acCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATc                                                                                                                                           <  2:4153957/149‑1 (MQ=255)
           ccTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCa                                                                                                                                          <  2:3952992/149‑1 (MQ=255)
            cTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCaa                                                                                                                                         <  1:4028309/149‑1 (MQ=255)
            cTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCaa                                                                                                                                         <  1:4228874/149‑1 (MQ=255)
             tgtgGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAg                                                                                                                                        <  1:1526881/149‑1 (MQ=255)
             tgtgGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAg                                                                                                                                        <  2:2605784/149‑1 (MQ=255)
               tgGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCa                                                                                                                                      <  1:3802679/149‑1 (MQ=255)
               tgGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCa                                                                                                                                      >  1:2495540/1‑149 (MQ=255)
                 gggACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCAtt                                                                                                                                    >  1:368100/1‑149 (MQ=255)
                 gggACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCAtt                                                                                                                                    <  1:5173563/149‑1 (MQ=255)
                    aCGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTg                                                                                                                                 <  2:4650874/149‑1 (MQ=255)
                     cGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGt                                                                                                                                <  2:4032223/149‑1 (MQ=255)
                      ggTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTc                                                                                                                               >  1:1148716/1‑149 (MQ=255)
                       gTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCa                                                                                                                              <  1:3482586/149‑1 (MQ=255)
                       gTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTCTCAAGCAGATCCGCCGTATCAAGCATTTTGTCa                                                                                                                              <  1:804309/149‑1 (MQ=255)
                         aaTTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATa                                                                                                                            <  1:4164618/149‑1 (MQ=255)
                             tgacGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCAGAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCa                                                                                                                        <  1:4420188/146‑1 (MQ=255)
                             tGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATATTGTCATAAGCa                                                                                                                        >  2:5056944/1‑149 (MQ=255)
                              gCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCAt                                                                                                                       <  2:4405714/149‑1 (MQ=255)
                               ccGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATc                                                                                                                      <  1:4085674/149‑1 (MQ=255)
                                  gTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGc                                                                                                                   <  2:2322927/149‑1 (MQ=255)
                                   ttttCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCt                                                                                                                  <  2:2455126/149‑1 (MQ=255)
                                    tttCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCtt                                                                                                                 <  2:5000524/149‑1 (MQ=255)
                                     ttCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTc                                                                                                                <  2:1494717/149‑1 (MQ=255)
                                      tCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTcc                                                                                                               <  1:3021838/149‑1 (MQ=255)
                                      tCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTcc                                                                                                               >  1:4287522/1‑149 (MQ=255)
                                       cAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCt                                                                                                              <  2:3656079/149‑1 (MQ=255)
                                        aGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCtt                                                                                                             <  1:4953312/149‑1 (MQ=255)
                                         gAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCttt                                                                                                            <  1:1353771/149‑1 (MQ=255)
                                          aaTGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCtttt                                                                                                           <  2:567258/149‑1 (MQ=255)
                                           aTGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCttttt                                                                                                          <  2:4757017/149‑1 (MQ=255)
                                               tGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTg                                                                                                      <  1:895647/149‑1 (MQ=255)
                                               tGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTg                                                                                                      <  1:2145144/149‑1 (MQ=255)
                                               tGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTg                                                                                                      >  1:1835552/1‑149 (MQ=255)
                                                       acaTCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTGTTTGCTGGTAAACGa                                                                                              >  2:4175524/1‑149 (MQ=255)
                                                        caTCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGAc                                                                                             <  1:3980697/149‑1 (MQ=255)
                                                          tatttttGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGATCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGTTTTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACAt                                                                                           <  2:2868626/147‑1 (MQ=255)
                                                            tttttGTTCTGCCAGCCATAGCGAAAGGTCTTCACGTTGCTCGTCTTCCATTTGCAATGGATAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACAttt                                                                                         <  2:2756921/149‑1 (MQ=255)
                                                             ttttGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACAtttt                                                                                        <  1:5060828/149‑1 (MQ=255)
                                                             ttttGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTATTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACAtttt                                                                                        <  2:2901713/149‑1 (MQ=255)
                                                                 gTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTctct                                                                                    >  1:4505924/1‑149 (MQ=255)
                                                                       ccagccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGccc                                                                              >  2:2404007/1‑149 (MQ=255)
                                                                        cagccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCt                                                                             <  2:1790074/149‑1 (MQ=255)
                                                                         agccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTc                                                                            >  1:5216420/1‑149 (MQ=255)
                                                                         agccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTc                                                                            >  1:2623654/1‑149 (MQ=255)
                                                                          gccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCa                                                                           >  2:811951/1‑149 (MQ=255)
                                                                           ccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCAc                                                                          >  1:1545939/1‑149 (MQ=255)
                                                                           ccaTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCAc                                                                          <  2:2153783/149‑1 (MQ=255)
                                                                            caTAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACg                                                                         >  2:3001469/1‑149 (MQ=255)
                                                                              tAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGaa                                                                       <  1:5156184/149‑1 (MQ=255)
                                                                               aGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAAt                                                                      <  1:1507583/149‑1 (MQ=255)
                                                                               aGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAAt                                                                      <  1:3754377/149‑1 (MQ=255)
                                                                               aGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAAt                                                                      <  1:3199633/149‑1 (MQ=255)
                                                                                gCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATg                                                                     <  2:885358/149‑1 (MQ=255)
                                                                                 cGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGa                                                                    <  2:3067393/149‑1 (MQ=255)
                                                                                     aGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTAc                                                                >  2:625942/1‑149 (MQ=255)
                                                                                         cTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTAt                                                            <  1:1377303/149‑1 (MQ=255)
                                                                                           tCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTAttt                                                          <  1:4910982/149‑1 (MQ=255)
                                                                                            cACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTa                                                         >  1:1209240/1‑149 (MQ=255)
                                                                                             aCGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTaa                                                        <  2:1150058/149‑1 (MQ=255)
                                                                                              cGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAAc                                                       <  2:125138/149‑1 (MQ=255)
                                                                                              cGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAAc                                                       <  2:2959880/149‑1 (MQ=255)
                                                                                              cGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAAc                                                       <  1:1640612/149‑1 (MQ=255)
                                                                                               gTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACt                                                      <  1:216020/149‑1 (MQ=255)
                                                                                                    tCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAAc                                                 <  2:752600/149‑1 (MQ=255)
                                                                                                    tCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAAc                                                 >  2:1322031/1‑149 (MQ=255)
                                                                                                    tCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAAc                                                 <  2:3426159/149‑1 (MQ=255)
                                                                                                     cGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAAcc                                                >  1:4821061/1‑149 (MQ=255)
                                                                                                       tCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGAAACAAGCATTAATTAATAAGCATCGGCTTCCTTTTTGCTTGTAAACTACATTATCTATACGCCATCACGAATCACTAAGTATATAAATAAAACCGc                                              >  1:4010012/1‑149 (MQ=255)
                                                                                                        cTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGcc                                             >  1:2651407/1‑149 (MQ=255)
                                                                                                        cTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGcc                                             <  1:2570982/149‑1 (MQ=255)
                                                                                                        cTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGcc                                             <  1:4728600/149‑1 (MQ=255)
                                                                                                         ttCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCa                                            <  2:696649/149‑1 (MQ=255)
                                                                                                          tCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCAt                                           <  2:5186229/149‑1 (MQ=255)
                                                                                                             aTTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTg                                        <  1:997648/149‑1 (MQ=255)
                                                                                                                  cACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGcc                                   <  1:3361322/149‑1 (MQ=255)
                                                                                                                   aCTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCct                                  >  2:5245305/1‑149 (MQ=255)
                                                                                                                   aCTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCct                                  <  2:399071/149‑1 (MQ=255)
                                                                                                                    cTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCctc                                 <  1:382485/149‑1 (MQ=255)
                                                                                                                     tGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCctct                                <  1:3577206/149‑1 (MQ=255)
                                                                                                                      ggAGAATTTGTCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCtctc                               >  2:453275/1‑149 (MQ=255)
                                                                                                                       gagaATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCa                              >  1:131912/1‑149 (MQ=255)
                                                                                                                            tttGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCataata                         <  2:1977918/149‑1 (MQ=255)
                                                                                                                            tttGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCataata                         >  2:2533831/1‑149 (MQ=255)
                                                                                                                               gTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACt                      <  2:4617697/149‑1 (MQ=255)
                                                                                                                                tCAGCCAGGTGTCAAGCAGAACCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACtg                     <  2:155122/149‑1 (MQ=255)
                                                                                                                                       ggTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGAtttt              <  1:3127584/149‑1 (MQ=255)
                                                                                                                                       ggTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGAtttt              >  2:5079735/1‑149 (MQ=255)
                                                                                                                                           tCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATAc          <  2:779747/149‑1 (MQ=255)
                                                                                                                                                 aGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGtata    <  1:1343888/149‑1 (MQ=255)
                                                                                                                                                  gATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGTatat   <  1:1535620/149‑1 (MQ=255)
                                                                                                                                                   aTCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGTATAtt  <  1:2243424/149‑1 (MQ=255)
                                                                                                                                                   aTCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGTATAtt  <  1:4446054/149‑1 (MQ=255)
                                                                                                                                                   aTCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGTATAtt  <  2:1607376/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGCACCGACGACCTGTGGGGACGGTAATTTGCCGGTTTTCAGAATGGTGCTCAGCACATCTTTTTGTTCTGCCAGCCATAGCGAAAGGGCTTCACGTTGCTCGTCTTCCATTTGCACTGGAGAATTTGTCAGCCAGGTGTCAAGCAGATCCGCCGTATCAAGCATTTTGTCATAAGCATCGGCTTCCTTTTTGCTGGTAAACGACATTTTCTCTTCGCCCTCACGAATGACTACGTATTTAACTTCAACCGCCATTTGCAGCCTCTCATAATAACTGTGATTTTATACAGTATATT  >  CP009273/1924750‑1925045

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: