Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 495,504 A→T 59.7% intergenic (‑75/‑77) aes ← / → gsk acetyl esterase/inosine/guanosine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273495,5040AT59.7% 251.0 / inf 298intergenic (‑75/‑77)aes/gskacetyl esterase/inosine/guanosine kinase
Reads supporting (aligned to +/- strand):  ref base A (62/58);  new base T (97/81);  total (159/139)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 302 total aligned reads displayed.

GAAGAGTATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTACTTCCCCGTAAACGCACGCGATCCGCTGCTTCAGCA  >  CP009273/495357‑495651
                                                                                                                                                   |                                                                                                                                                   
 aaGAGTATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGaaa                                                                                                                                                   >  1:751843/1‑149 (MQ=255)
 aaGAGTATTCACAACGGTCCTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGAGATTGaaa                                                                                                                                                   >  1:3994071/1‑149 (MQ=255)
   gagTATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACt                                                                                                                                                 <  1:952552/149‑1 (MQ=255)
    agTATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTg                                                                                                                                                <  2:4569659/149‑1 (MQ=255)
    agTATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTg                                                                                                                                                <  2:3944267/149‑1 (MQ=255)
    agTATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTg                                                                                                                                                <  2:1976303/149‑1 (MQ=255)
      tATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCAGCCTGCTGAAATTGAAACTGGc                                                                                                                                              >  1:3252281/1‑149 (MQ=255)
       aTTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCt                                                                                                                                             >  1:3324368/1‑149 (MQ=255)
       aTTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCt                                                                                                                                             >  1:2660296/1‑149 (MQ=255)
          cacaACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCg                                                                                                                                          >  1:836434/1‑149 (MQ=255)
           acaACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGt                                                                                                                                         >  2:2291560/1‑149 (MQ=255)
             aaCGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGtct                                                                                                                                       >  1:4246438/1‑149 (MQ=255)
              aCGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGtctc                                                                                                                                      >  1:5224952/1‑149 (MQ=255)
               cGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTcg                                                                                                                                     <  1:495343/149‑1 (MQ=255)
                    ttCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCgcgcgc                                                                                                                                <  1:2930716/149‑1 (MQ=255)
                     tCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGCAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCt                                                                                                                               >  1:250228/1‑149 (MQ=255)
                        tttCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTccc                                                                                                                            <  1:2179174/149‑1 (MQ=255)
                         ttCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCg                                                                                                                           <  1:1745570/149‑1 (MQ=255)
                             gcagAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAg                                                                                                                       >  1:2390410/1‑149 (MQ=255)
                               agAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGAt                                                                                                                     <  1:1263213/149‑1 (MQ=255)
                                gAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGAtt                                                                                                                    <  1:5098112/149‑1 (MQ=255)
                                     aaGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGtt                                                                                                               <  2:1278937/149‑1 (MQ=255)
                                         tCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCTAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATTTTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACa                                                                                                           <  2:1761527/149‑1 (MQ=255)
                                         tCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTGCGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACa                                                                                                           <  2:3284474/149‑1 (MQ=255)
                                          ccAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACAt                                                                                                          >  2:580729/1‑149 (MQ=255)
                                          ccAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACAt                                                                                                          <  1:3416713/149‑1 (MQ=255)
                                           cAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACAtt                                                                                                         >  1:1005283/1‑149 (MQ=255)
                                            aGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTc                                                                                                        <  2:1597939/149‑1 (MQ=255)
                                              aaCAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTcgc                                                                                                      <  1:2941952/149‑1 (MQ=255)
                                               aCAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTcgcc                                                                                                     <  2:751843/149‑1 (MQ=255)
                                                cAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTcgccg                                                                                                    >  1:3855433/1‑149 (MQ=255)
                                                 aGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTcgccgc                                                                                                   >  2:2211213/1‑149 (MQ=255)
                                                 aGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTcgccgc                                                                                                   <  1:2032140/149‑1 (MQ=255)
                                                  ggTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCt                                                                                                  >  2:1031892/1‑149 (MQ=255)
                                                        ttgtttTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATACTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTa                                                                                            >  1:1481952/1‑149 (MQ=255)
                                                          gtttTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAac                                                                                          >  1:4549525/1‑149 (MQ=255)
                                                              tCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAaccacc                                                                                      <  1:4286241/149‑1 (MQ=255)
                                                                    ttCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTaaa                                                                                >  1:4887824/1‑149 (MQ=255)
                                                                        tACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAaaca                                                                            <  1:3574330/149‑1 (MQ=255)
                                                                        tACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAaaca                                                                            <  1:2124003/149‑1 (MQ=255)
                                                                         aCAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAaacaa                                                                           <  1:1307617/149‑1 (MQ=255)
                                                                          cAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAaacaac                                                                          <  2:4331681/149‑1 (MQ=255)
                                                                          cAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAaacaac                                                                          <  2:1406978/149‑1 (MQ=255)
                                                                              aaCTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATg                                                                      >  1:19090/1‑149 (MQ=255)
                                                                                 tCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGaaa                                                                   >  1:282322/1‑149 (MQ=255)
                                                                                   cTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAAtt                                                                 >  1:4417299/1‑149 (MQ=255)
                                                                                    tttCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAAttt                                                                >  2:5229808/1‑149 (MQ=255)
                                                                                    tttCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAAttt                                                                >  2:2608561/1‑149 (MQ=255)
                                                                                     ttCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTc                                                               <  1:668985/149‑1 (MQ=255)
                                                                                     ttCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTc                                                               <  2:3562913/149‑1 (MQ=255)
                                                                                     ttCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTc                                                               <  2:3089946/149‑1 (MQ=255)
                                                                                     ttCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTc                                                               >  1:141128/1‑149 (MQ=255)
                                                                                        aaaTTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCg                                                            >  2:644469/1‑149 (MQ=255)
                                                                                         aaTTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCgg                                                           >  1:3887141/1‑149 (MQ=255)
                                                                                         aaTTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCgg                                                           <  2:2673554/149‑1 (MQ=255)
                                                                                          aTTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGt                                                          <  1:2744813/149‑1 (MQ=255)
                                                                                          aTTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGt                                                          >  1:1696796/1‑149 (MQ=255)
                                                                                           ttACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTa                                                         >  1:4053491/1‑149 (MQ=255)
                                                                                           ttACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTa                                                         <  2:214145/149‑1 (MQ=255)
                                                                                            tACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTaa                                                        >  2:4819480/1‑149 (MQ=255)
                                                                                             aCGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTaaa                                                       <  2:3324368/149‑1 (MQ=255)
                                                                                             aCGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTaaa                                                       <  2:4835791/149‑1 (MQ=255)
                                                                                              cGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAAc                                                      <  2:494333/149‑1 (MQ=255)
                                                                                              cGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAAc                                                      <  1:1146130/149‑1 (MQ=255)
                                                                                                   ttGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTaaa                                                 >  1:4741989/1‑149 (MQ=255)
                                                                                                   ttGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTaaa                                                 <  1:2986982/149‑1 (MQ=255)
                                                                                                   ttGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTaaa                                                 >  2:3028566/1‑149 (MQ=255)
                                                                                                    tGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAAt                                                <  2:3228126/149‑1 (MQ=255)
                                                                                                    tGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAAt                                                <  2:5030673/149‑1 (MQ=255)
                                                                                                     gTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACAGTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATTAAATTTCCCGGTAAACGTAAATc                                               <  2:2970723/149‑1 (MQ=255)
                                                                                                        aGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCaa                                            >  2:4525181/1‑149 (MQ=255)
                                                                                                         ggAACCACTGCCATGAAAATGCGATCCCGCCAGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCaaa                                           >  1:2867317/1‑149 (MQ=255)
                                                                                                            aCCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACAt                                        <  1:398164/149‑1 (MQ=255)
                                                                                                            aCCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACAt                                        <  2:2825430/149‑1 (MQ=255)
                                                                                                            aCCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACAt                                        >  2:3529925/1‑149 (MQ=255)
                                                                                                             ccACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACAtt                                       >  2:3425529/1‑149 (MQ=255)
                                                                                                             ccACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACAtt                                       <  1:462555/149‑1 (MQ=255)
                                                                                                             ccACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACAtt                                       >  2:327220/1‑149 (MQ=255)
                                                                                                              cACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTa                                      <  2:2385768/149‑1 (MQ=255)
                                                                                                              cACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTa                                      <  1:2229310/149‑1 (MQ=255)
                                                                                                              cACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTa                                      >  1:2087031/1‑149 (MQ=255)
                                                                                                               aCTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTAc                                     <  1:3759717/149‑1 (MQ=255)
                                                                                                                cTGCCATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTACt                                    >  2:4725299/1‑149 (MQ=255)
                                                                                                                cTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTACt                                    <  2:2938626/149‑1 (MQ=255)
                                                                                                                cTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTACt                                    >  2:2724871/1‑149 (MQ=255)
                                                                                                                cTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTACt                                    <  1:1883543/149‑1 (MQ=255)
                                                                                                                    cATGAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTACTTccc                                <  2:408297/149‑1 (MQ=255)
                                                                                                                       gAAAATGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTACTTCCCCGt                             >  2:3989182/1‑149 (MQ=255)
                                                                                                                        aaaaTGCGATCCCGCCTGCTGATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTACTTCCCCGTa                            >  2:4294150/1‑149 (MQ=255)
                                                                                                                          aaTGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTACTTCCCCGTaaa                          <  1:3822366/149‑1 (MQ=255)
                                                                                                                           aTGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTACTTCCCCGTAAac                         <  2:3639593/149‑1 (MQ=255)
                                                                                                                           aTGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTACTTCCCCGTAAac                         >  2:1917762/1‑149 (MQ=255)
                                                                                                                            tGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTACTTCCCCGTAAacg                        >  2:380336/1‑149 (MQ=255)
                                                                                                                              cGATCCCGTCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTACTTCCCCGTAAacgca                      <  1:1551259/149‑1 (MQ=255)
                                                                                                                                 tCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTACTTCCCCGTAAacgcacgc                   <  2:1187514/149‑1 (MQ=255)
                                                                                                                                  cccGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTACTTCCCCGTAAACGCAcgcg                  >  2:87581/1‑149 (MQ=255)
                                                                                                                                        tgctgATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTACTTCCCCGTAAACGCACGCGATCCgc            >  2:1009072/1‑149 (MQ=255)
                                                                                                                                         gctgATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTACTTCCCCGTAAACGCACGCGATCCgct           >  2:1670619/1‑149 (MQ=255)
                                                                                                                                            gATATTGTAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTACTTCCCCGTAAACGCACGCGATCCgctgct        <  1:4160349/149‑1 (MQ=255)
                                                                                                                                                  gAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTACTTCCCCGTAAACGCACGCGATCCGCTGCTTCAGCa  >  2:1124860/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GAAGAGTATTCACAACGGTCTTCATTTCAGCAGAAATAAGGTCCAGAACAGGTAGTTTGTTTTCCGGCTTCATACAAAAACTCCTTTCAAATTACGTCATTGTAAGGAACCACTGCCATGAAAATGCGATCCCGCCTGCTGATATTGAAACTGGCTGCGTCTCGCGCGCTCCCGTCAGATTGTGTTAACATTCGCCGCTCAGTTAACCACCCGTAAAAACAACCATGAAATTTCCCGGTAAACGTAAATCCAAACATTACTTCCCCGTAAACGCACGCGATCCGCTGCTTCAGCA  >  CP009273/495357‑495651

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: