Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,116,874 C→A 27.8% E24* (GAA→TAA)  dinI ← DNA damage‑inducible protein I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,116,8740CA27.8% 652.2 / 207.5 300E24* (GAA→TAA) dinIDNA damage‑inducible protein I
Reads supporting (aligned to +/- strand):  ref base C (85/131);  new base A (30/54);  total (115/185)
Fisher's exact test for biased strand distribution p-value = 5.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 301 total aligned reads displayed.

CTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATT  >  CP009273/1116726‑1117021
                                                                                                                                                    |                                                                                                                                                   
cTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCACTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTc                                                                                                                                                     <  2:2673931/149‑1 (MQ=255)
   tCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCgcc                                                                                                                                                  <  1:1895070/149‑1 (MQ=255)
   tCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTAgcc                                                                                                                                                  <  2:145915/149‑1 (MQ=255)
    cccACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCgcca                                                                                                                                                 >  2:4296332/1‑149 (MQ=255)
    cccACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCgcca                                                                                                                                                 <  1:4659292/149‑1 (MQ=255)
    cccACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTAgcca                                                                                                                                                 >  1:3891948/1‑149 (MQ=255)
       aCGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCgccagcc                                                                                                                                              >  1:1039697/1‑149 (MQ=255)
       aCGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCgccagcc                                                                                                                                              >  1:2017775/1‑149 (MQ=255)
        cGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCgccagcca                                                                                                                                             <  2:124867/149‑1 (MQ=255)
         gTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGccagccag                                                                                                                                            <  1:2976382/149‑1 (MQ=255)
         gTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGccagccag                                                                                                                                            <  1:3778083/149‑1 (MQ=255)
         gTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTAGccagccag                                                                                                                                            >  1:1098606/1‑149 (MQ=255)
            tCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGc                                                                                                                                         <  1:410125/149‑1 (MQ=255)
              cTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTAGCCAGCCAGGGCGt                                                                                                                                       >  1:4697234/1‑149 (MQ=255)
                ggAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCa                                                                                                                                     >  1:1544522/1‑149 (MQ=255)
                ggAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCa                                                                                                                                     >  1:3052940/1‑149 (MQ=255)
                 gagaATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCaa                                                                                                                                    <  2:1590481/149‑1 (MQ=255)
                 gagaATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTAGCCAGCCAGGGCGTCaa                                                                                                                                    <  1:919821/149‑1 (MQ=255)
                  agaATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAAt                                                                                                                                   <  1:2545160/149‑1 (MQ=255)
                  agaATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTAGCCAGCCAGGGCGTCAAt                                                                                                                                   <  1:2211334/149‑1 (MQ=255)
                   gaATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATa                                                                                                                                  >  2:4167459/1‑149 (MQ=255)
                     aTTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTAGCCAGCCAGGGCGTCAATAGc                                                                                                                                <  2:2349207/149‑1 (MQ=255)
                      tttCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGcc                                                                                                                               <  1:1988325/149‑1 (MQ=255)
                        tCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGcccc                                                                                                                             <  2:3655346/149‑1 (MQ=255)
                          aCTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTAGCCAGCCAGGGCGTCAATAGCCCCAg                                                                                                                           >  1:705597/1‑149 (MQ=255)
                           cTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTAGCCAGCCAGGGCGTCAATAGCCCCAGc                                                                                                                          <  1:3669456/149‑1 (MQ=255)
                             aaTGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTg                                                                                                                        >  1:2368085/1‑149 (MQ=255)
                               tGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGc                                                                                                                      >  2:2563066/1‑149 (MQ=255)
                               tGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGc                                                                                                                      >  1:4622179/1‑149 (MQ=255)
                                gcgcTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTAGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCa                                                                                                                     <  2:694628/149‑1 (MQ=255)
                                 cgcTGTTTATCTTCTTTTGTCGCGCCACTAACCGAGAAATTAGTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCaa                                                                                                                    <  1:2633/149‑1 (MQ=255)
                                  gcTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTAGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAAt                                                                                                                   <  1:4431663/149‑1 (MQ=255)
                                   cTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTAGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATg                                                                                                                  <  1:354637/149‑1 (MQ=255)
                                    tGTTTATCTTCTTTTGTCGCACCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATgg                                                                                                                 <  2:4611401/149‑1 (MQ=255)
                                     gTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGGACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGga                                                                                                                <  2:836246/149‑1 (MQ=255)
                                      tttATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGgag                                                                                                               <  1:3381679/149‑1 (MQ=255)
                                        tATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGaa                                                                                                             >  1:1282559/1‑149 (MQ=255)
                                          tcttctTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGt                                                                                                           >  1:1435489/1‑149 (MQ=255)
                                          tcttctTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTAGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGt                                                                                                           <  1:4171394/149‑1 (MQ=255)
                                            ttcttTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTAGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGttt                                                                                                         <  1:4342854/149‑1 (MQ=255)
                                             tcttTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGtttt                                                                                                        >  1:3137733/1‑149 (MQ=255)
                                              cttTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTc                                                                                                       >  2:1936000/1‑149 (MQ=255)
                                                     cgcgCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTAGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATgg                                                                                                <  1:4492868/149‑1 (MQ=255)
                                                       cgcCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTg                                                                                              >  2:1167181/1‑149 (MQ=255)
                                                         ccAATAACCGATAAATGATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGAc                                                                                            <  1:1235327/149‑1 (MQ=255)
                                                            aTAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGACGTTTTCGCTATGGTGACTTc                                                                                         >  1:1695319/1‑149 (MQ=255)
                                                            aTAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTAGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTc                                                                                         >  2:3840124/1‑149 (MQ=255)
                                                              atccGATATATTATTCGATGCGGCTTTACGTAGCGATACGTGGACTTCATTATCTGGAAACTCATTCTGAATTCGGCGGGAAAGTTAGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGAtttata                                                                                       <  2:303568/147‑4 (MQ=255)
                                                                ccGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAAtt                                                                                     <  1:4261041/149‑1 (MQ=255)
                                                                   aTAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGACTACGGCGGGCAAGTTAGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGc                                                                                  <  2:289506/149‑1 (MQ=255)
                                                                    tAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTAGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCa                                                                                 >  2:3580981/1‑149 (MQ=255)
                                                                    tAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTAGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCa                                                                                 <  2:360963/149‑1 (MQ=255)
                                                                     aaaTTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCat                                                                                <  2:1642475/149‑1 (MQ=255)
                                                                     aaaTTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCat                                                                                <  2:3529486/149‑1 (MQ=255)
                                                                               cTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGccccc                                                                      <  1:1764203/149‑1 (MQ=255)
                                                                               cTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGccccc                                                                      <  1:3014108/149‑1 (MQ=255)
                                                                                tGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCTCATAATAGCCCCCt                                                                     <  1:2451074/149‑1 (MQ=255)
                                                                                tGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTAGCAAAATAGCCCCCt                                                                     >  1:4432421/1‑149 (MQ=255)
                                                                                  cGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGt                                                                   <  1:777173/149‑1 (MQ=255)
                                                                                      aTAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGaa                                                               <  1:4747392/149‑1 (MQ=255)
                                                                                        aaCGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAata                                                             <  1:1025228/149‑1 (MQ=255)
                                                                                         aCGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTACGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAatat                                                            <  1:3188371/149‑1 (MQ=255)
                                                                                         aCGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTAGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTAAATTAGCATAATAGCCCCCAGTTGAatat                                                            >  1:1721893/1‑149 (MQ=255)
                                                                                            tACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATATCCCCCTGTTGAATATGCt                                                         >  1:2273734/1‑149 (MQ=255)
                                                                                              ccGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCAGTTGAATATGCTgg                                                       >  1:2562098/1‑149 (MQ=255)
                                                                                               cGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATATCCCCCTGTTGAATATGCTGGt                                                      >  2:4021115/1‑149 (MQ=255)
                                                                                               cGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGt                                                      <  2:1990332/149‑1 (MQ=255)
                                                                                               cGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTAGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGt                                                      <  2:1669811/149‑1 (MQ=255)
                                                                                                 aTACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTAGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGtta                                                    <  1:205789/149‑1 (MQ=255)
                                                                                                     gTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGttattta                                                >  2:340906/1‑149 (MQ=255)
                                                                                                        gCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATAc                                             <  2:2824340/149‑1 (MQ=255)
                                                                                                          cTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAg                                           <  2:2501147/149‑1 (MQ=255)
                                                                                                          cTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAg                                           <  2:2258657/149‑1 (MQ=255)
                                                                                                          cTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTAGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAg                                           <  2:3891948/149‑1 (MQ=255)
                                                                                                             cATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTa                                        <  2:593790/149‑1 (MQ=255)
                                                                                                                  tCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATa                                   <  2:795303/149‑1 (MQ=255)
                                                                                                                    aGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAAc                                 <  2:3688671/149‑1 (MQ=255)
                                                                                                                      gAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCt                               <  2:1853087/149‑1 (MQ=255)
                                                                                                                         aCGCATACTGAATACGGCGGGAAAGTTAGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAAt                            <  2:3144965/149‑1 (MQ=255)
                                                                                                                          cGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTAGCTATTGTGACTTCAATTCGCATAATAACCCCCAGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATg                           >  1:3995005/1‑149 (MQ=255)
                                                                                                                           gCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTGTTCGCTATGGTGCCTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGGAAAATAACCTAATGa                          >  2:554081/1‑149 (MQ=255)
                                                                                                                           gCATACTGAATACGGCGGGAAAGTTAGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTGTTCGCTATGGTGACTTCAATTAGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAAACAAATgt                          >  1:4052751/1‑148 (MQ=255)
                                                                                                                             aTACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGaca                        >  2:2496934/1‑149 (MQ=255)
                                                                                                                              tACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGacaa                       >  2:4549684/1‑149 (MQ=255)
                                                                                                                              tACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTAGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGacaa                       >  2:1416428/1‑149 (MQ=255)
                                                                                                                                 tGAATACGGCGGGAAAGTTAGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAg                    >  1:1802966/1‑149 (MQ=255)
                                                                                                                                  gAATACGGCGGGAAAGTTAGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAgg                   <  1:2756038/149‑1 (MQ=255)
                                                                                                                                   caTACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGa                  <  2:519666/148‑1 (MQ=255)
                                                                                                                                   aagaCGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGa                  <  2:3114331/146‑1 (MQ=255)
                                                                                                                                    agaCGGCGGGAAAGTTAGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGaa                 <  2:4240667/147‑1 (MQ=255)
                                                                                                                                     tACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAg                <  2:2207642/149‑1 (MQ=255)
                                                                                                                                        ggcggGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTa             >  2:4350251/1‑149 (MQ=255)
                                                                                                                                           gggAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGa          <  2:4622643/149‑1 (MQ=255)
                                                                                                                                            ggAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGAt         >  1:4306556/1‑149 (MQ=255)
                                                                                                                                             gAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGAtt        <  1:2869183/149‑1 (MQ=255)
                                                                                                                                                aGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGAttttt     <  1:528907/149‑1 (MQ=255)
                                                                                                                                                  ttCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTAt   <  1:4199577/149‑1 (MQ=255)
                                                                                                                                                  gtCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTAt   <  2:3210952/148‑1 (MQ=255)
                                                                                                                                                   tCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTAtt  <  2:2017775/149‑1 (MQ=255)
                                                                                                                                                   tCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATGCTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTAtt  <  2:2841726/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CTTTCCCACGTTTCCTGGAGAATTTCACTAATGCGCTGTTTATCTTCTTTTGTCGCGCCAATAACCGATAAATTATTCGCTGCGGCATAACGTACCGATACGTGGCCTTCATTATCAGGAAACGCATACTGAATACGGCGGGAAAGTTCGCCAGCCAGGGCGTCAATAGCCCCAGCTGGCAATGGAGAAGTTTTCGCTATGGTGACTTCAATTCGCATAATAGCCCCCTGTTGAATATACTGGTTATTTATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATT  >  CP009273/1116726‑1117021

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: