Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 515,004 G→T 35.2% S73* (TCG→TAG)  tesA ← multifunctional acyl‑CoA thioesterase 1; protease I; lysophospholipase L1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273515,0040GT35.2% 427.5 / 239.1 257S73* (TCG→TAG) tesAmultifunctional acyl‑CoA thioesterase 1; protease I; lysophospholipase L1
Reads supporting (aligned to +/- strand):  ref base G (74/92);  new base T (43/48);  total (117/140)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 259 total aligned reads displayed.

TTGACATCCTGCAAAATCTGGCGCAGCGTTTGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGGCGCTCAGGCTATCACCCAGAATCAATAACGTGTCCGCTGCGGCGGCA  >  CP009273/514857‑515153
                                                                                                                                                   |                                                                                                                                                     
ttGACATCCTGCAAAATCTGGCGCAGCGTTTGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAg                                                                                                                                                     <  2:3526988/150‑1 (MQ=255)
ttGACATCCTGCAAAATCTGGCGCAGCGTTTGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAg                                                                                                                                                     >  2:1518306/1‑150 (MQ=255)
 tGACATCCTGCAAAATCTGGCGCAGCGTTTGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAgg                                                                                                                                                    >  1:2648582/1‑150 (MQ=255)
   aCATCCTGCAAAATCTGGCGCAGCGTTTGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGgtg                                                                                                                                                  >  2:1412266/1‑150 (MQ=255)
     aTCCTGCAAAATCTGGCGCAGCGTTTGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTc                                                                                                                                                >  2:1157018/1‑150 (MQ=255)
     aTCCTGCAAAATCTGGCGCAGCGTTTGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTc                                                                                                                                                >  2:787265/1‑150 (MQ=255)
     aTCCTGCAAAATCTGGCGCAGCGTTTGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTc                                                                                                                                                >  2:3737943/1‑150 (MQ=255)
      tCCTGCAAAATCTGGCGCAGCGTTTGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTcg                                                                                                                                               <  2:2453066/150‑1 (MQ=255)
      tCCTGCAAAATCTGGCGCAGCGTTTGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTcg                                                                                                                                               <  2:619959/150‑1 (MQ=255)
        cTGCAAAATCTGGCGCAGCGTTTGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTcgcc                                                                                                                                             <  2:446567/150‑1 (MQ=255)
          gCAAAATCTGGCGCAGCGTTTGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTcgccgc                                                                                                                                           >  2:3894104/1‑150 (MQ=255)
          gCAAAATCTGGCGCAGCGTTTGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTcgccgc                                                                                                                                           <  1:3560287/150‑1 (MQ=255)
            aaaaTCTGGCGCAGCGTTTGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTg                                                                                                                                         >  2:2316150/1‑150 (MQ=255)
              aaTCTGGCGCAGCGTTTGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGAt                                                                                                                                       <  2:2200443/150‑1 (MQ=255)
               aTCTGGCGCAGCGTTTGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATg                                                                                                                                      <  2:396776/150‑1 (MQ=255)
                 cTGGCGCAGCGTTTGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCt                                                                                                                                    <  2:2942885/150‑1 (MQ=255)
                   ggCGCAGCGTTTGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTgg                                                                                                                                  <  2:2269050/150‑1 (MQ=255)
                        aGCGTTTGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTa                                                                                                                             <  1:1885997/150‑1 (MQ=255)
                          cGTTTGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAAc                                                                                                                           <  1:2960147/150‑1 (MQ=255)
                             ttGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTAc                                                                                                                        >  1:3423882/1‑150 (MQ=255)
                              tGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTAcc                                                                                                                       <  2:3422740/150‑1 (MQ=255)
                              tGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTAcc                                                                                                                       <  2:1627981/150‑1 (MQ=255)
                              tGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTAcc                                                                                                                       >  1:483235/1‑150 (MQ=255)
                               gCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCg                                                                                                                      >  1:67805/1‑150 (MQ=255)
                                cTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGa                                                                                                                     <  1:1016042/150‑1 (MQ=255)
                                cTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGa                                                                                                                     <  1:4609971/150‑1 (MQ=255)
                                  cGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACg                                                                                                                   <  1:2438236/150‑1 (MQ=255)
                                  cGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACg                                                                                                                   <  2:1163145/150‑1 (MQ=255)
                                   ggTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGt                                                                                                                  >  1:2065144/1‑150 (MQ=255)
                                    gTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGtt                                                                                                                 >  2:3605917/1‑150 (MQ=255)
                                    gTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGtt                                                                                                                 <  2:1902754/150‑1 (MQ=255)
                                    gTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGtt                                                                                                                 >  2:1674845/1‑150 (MQ=255)
                                         cTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGCAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGGTGCTGGCAGGAACTACCGACGGGTTAct                                                                                                            >  2:1648496/1‑150 (MQ=255)
                                           gtgGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTActct                                                                                                          <  1:888697/150‑1 (MQ=255)
                                           gtgGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTActct                                                                                                          <  1:1342676/150‑1 (MQ=255)
                                            tgGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTg                                                                                                         <  2:1030595/150‑1 (MQ=255)
                                                tGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCAc                                                                                                     >  1:3972374/1‑150 (MQ=255)
                                                  aaaaCCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACtt                                                                                                   >  1:2797201/1‑150 (MQ=255)
                                                  aaaaCCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACtt                                                                                                   <  1:1850289/150‑1 (MQ=255)
                                                   aaaCCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTa                                                                                                  >  1:2505649/1‑150 (MQ=255)
                                                      ccACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCa                                                                                               <  1:2673574/150‑1 (MQ=255)
                                                      ccACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCa                                                                                               >  2:2927945/1‑150 (MQ=255)
                                                       cACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCAt                                                                                              <  2:4225929/150‑1 (MQ=255)
                                                        aCGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCAtt                                                                                             <  1:1217475/150‑1 (MQ=255)
                                                        aCGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCAtt                                                                                             >  2:2713107/1‑150 (MQ=255)
                                                          gCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTca                                                                                           <  2:677642/150‑1 (MQ=255)
                                                           cAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTcaa                                                                                          <  2:1176370/150‑1 (MQ=255)
                                                           cAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTcaa                                                                                          >  2:3151067/1‑150 (MQ=255)
                                                             aaCCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTcaaca                                                                                        >  1:3230563/1‑150 (MQ=255)
                                                              aCCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTcaacaa                                                                                       <  2:2009790/150‑1 (MQ=255)
                                                              aCCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTcaacaa                                                                                       >  1:3303732/1‑150 (MQ=255)
                                                                cGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAagg                                                                                     >  1:2289481/1‑150 (MQ=255)
                                                                cGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAagg                                                                                     <  2:4007143/150‑1 (MQ=255)
                                                                cGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAagg                                                                                     <  2:1548014/150‑1 (MQ=255)
                                                                cGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAagg                                                                                     <  1:2342414/150‑1 (MQ=255)
                                                                 gTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAaggc                                                                                    <  2:4150670/150‑1 (MQ=255)
                                                                   cATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAaggcag                                                                                  <  1:997226/150‑1 (MQ=255)
                                                                    aTTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAaggcagg                                                                                 >  2:4598213/1‑150 (MQ=255)
                                                                       gccgccCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGcaggcca                                                                              <  2:1056547/150‑1 (MQ=255)
                                                                         cgccCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGcaggccagg                                                                            <  2:4496566/150‑1 (MQ=255)
                                                                          gccCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGcaggccaggc                                                                           >  1:1718512/1‑150 (MQ=255)
                                                                            ccAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCcg                                                                         <  1:3051804/150‑1 (MQ=255)
                                                                            ccAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCcg                                                                         <  1:1555350/150‑1 (MQ=255)
                                                                             cAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCcgc                                                                        >  1:4523992/1‑150 (MQ=255)
                                                                              aGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCcgcg                                                                       >  2:1201233/1‑150 (MQ=255)
                                                                                      cAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGa                                                               >  1:2070959/1‑150 (MQ=255)
                                                                                        gCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACa                                                             >  2:3683206/1‑150 (MQ=255)
                                                                                           cccAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTc                                                          <  2:3938879/150‑1 (MQ=255)
                                                                                           cccAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTc                                                          <  2:1828601/150‑1 (MQ=255)
                                                                                            ccAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCg                                                         <  1:1627729/150‑1 (MQ=255)
                                                                                            ccAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCg                                                         >  2:3858103/1‑150 (MQ=255)
                                                                                            ccAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCg                                                         <  1:41715/150‑1 (MQ=255)
                                                                                             cAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGa                                                        >  1:4346117/1‑150 (MQ=255)
                                                                                              aaCGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTAGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGAt                                                       >  2:3061992/1‑150 (MQ=255)
                                                                                                cgcgGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATAc                                                     <  1:4638392/150‑1 (MQ=255)
                                                                                                  cgGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATAccc                                                   <  2:2093804/150‑1 (MQ=255)
                                                                                                  cgGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATAccc                                                   >  1:2342511/1‑150 (MQ=255)
                                                                                                     cTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGgc                                                <  1:3894104/150‑1 (MQ=255)
                                                                                                     cTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGgc                                                >  2:1095969/1‑150 (MQ=255)
                                                                                                      tGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGgcg                                               >  2:3726642/1‑150 (MQ=255)
                                                                                                      tGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGgcg                                               >  1:2471296/1‑150 (MQ=255)
                                                                                                         tGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGGCGCTc                                            <  2:2648310/150‑1 (MQ=255)
                                                                                                         tGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGGCGCTc                                            >  2:3779878/1‑150 (MQ=255)
                                                                                                          gCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGGCGCTCa                                           <  2:4066524/150‑1 (MQ=255)
                                                                                                                    cagAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGGCGCTCAGGCTATCAcc                                 <  1:4132095/150‑1 (MQ=255)
                                                                                                                          cGGAAGGCGCGCCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGGCGCTCAGGCTATCACCCAGAAt                           <  1:1441086/150‑1 (MQ=255)
                                                                                                                          cGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGGCGCTCAGGCTATCACCCAGAAt                           <  1:3958460/150‑1 (MQ=255)
                                                                                                                            gAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGGCGCTCAGGCTATCACCCAGAATCa                         <  2:4319682/150‑1 (MQ=255)
                                                                                                                                    gcCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGGCGCTCAGGCTATCACCCAGAATCAATAACgtg                 >  1:255255/1‑150 (MQ=255)
                                                                                                                                    gcCAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGGCGCTCAGGCTATCACCCAGAATCAATAACgtg                 >  1:2572108/1‑150 (MQ=255)
                                                                                                                                     ccAGTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGGCGCTCAGGCTATCACCCAGAATCAATAACgtgt                <  1:729902/150‑1 (MQ=255)
                                                                                                                                       aGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGGCGCTCAGGCTATCACCCAGAATCAATAACGTGTcc              <  2:1881850/150‑1 (MQ=255)
                                                                                                                                       aGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGGCGCTCAGGCTATCACCCAGAATCAATAACGTGTcc              >  1:4603216/1‑150 (MQ=255)
                                                                                                                                        gTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGGCGCTCAGGCTATCACCCAGAATCAATAACGTGTCCg             >  1:38534/1‑150 (MQ=255)
                                                                                                                                        gTCCTTGTTGCTAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGGCGCTCAGGCTATCACCCAGAATCAATAACGTGTCCg             <  2:1337950/150‑1 (MQ=255)
                                                                                                                                        gTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGGCGCTCAGGCTATCACCCAGAATCAATAACGTGTCCg             <  1:4291908/150‑1 (MQ=255)
                                                                                                                                        gTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGACAGGCCTCGCTGGCAGACATTCGATACCCGGCGCTCAGGCTATCACCCAGAATCAATAACGTGTCCg             >  2:1109495/1‑150 (MQ=255)
                                                                                                                                               ttgCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGGCGCTCAGGCTATCACCCAGAATCAATAACGTGTCCGCTgcggc      >  1:3408337/1‑150 (MQ=255)
                                                                                                                                                  cGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGGCGCTCAGGCTATCACCCAGAATCAATAACGTGTCCGCTGCggcggc   >  2:2140981/1‑150 (MQ=255)
                                                                                                                                                   gAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGGCGCTCAGGCTATCACCCAGAATCAATAACGTGTCCGCTGCGGCGGCa  <  1:4137464/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTGACATCCTGCAAAATCTGGCGCAGCGTTTGCTCGGTTTGCTGTGGCTGAAAACCACGCAAACCGTCATTGCCGCCCAGTTCAACCAGCACCCAACGCGGCTGATGCTGTTTCAGCAGAGCCGGAAGGCGCGCCAGTCCTTGTTGCGAGGTGTCGCCGCTGATGCTGGCATTAACTACCGACGTTTTACTCTGCCACTTATCATTCAACAAGGCAGGCCAGGCCGCGCTGGCAGACATTCGATACCCGGCGCTCAGGCTATCACCCAGAATCAATAACGTGTCCGCTGCGGCGGCA  >  CP009273/514857‑515153

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: