Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,896,944 C→A 35.3% V40L (GTA→TTA)  bglB ← cryptic phospho‑beta‑glucosidase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,896,9440CA35.3% 437.6 / 271.5 275V40L (GTA→TTA) bglBcryptic phospho‑beta‑glucosidase B
Reads supporting (aligned to +/- strand):  ref base C (66/112);  new base A (51/46);  total (117/158)
Fisher's exact test for biased strand distribution p-value = 1.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-01

Only 100 of 277 total aligned reads displayed.

AATTCGCGCCCAGGCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATG  >  CP009273/3896795‑3897092
                                                                                                                                                     |                                                                                                                                                    
aaTTCGCGCCCAGGCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAc                                                                                                                                                      <  2:2938638/150‑1 (MQ=255)
 aTTCGCGCCCAGGCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACg                                                                                                                                                     <  2:756980/150‑1 (MQ=255)
 aTTCGCGCCCAGGCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACg                                                                                                                                                     >  1:4537558/1‑150 (MQ=255)
 aTTCGCGCCCAGGCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACg                                                                                                                                                     >  1:792825/1‑150 (MQ=255)
    cgcgCCCAGGCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCa                                                                                                                                                  <  1:3258093/150‑1 (MQ=255)
    cgcgCCCAGGCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCa                                                                                                                                                  <  1:2346704/150‑1 (MQ=255)
     gcgcCCAGGCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCAt                                                                                                                                                 <  1:4230749/150‑1 (MQ=255)
      cgcCCAGGCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATg                                                                                                                                                <  1:3871045/150‑1 (MQ=255)
        cacAGGAAATTGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACTCGATATCTTCAGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGat                                                                                                                                              <  2:2348375/148‑2 (MQ=255)
         ccAGGCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAgg                                                                                                                                             >  1:2249613/1‑150 (MQ=255)
         ccAGGCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGAAGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAgg                                                                                                                                             >  1:4578109/1‑150 (MQ=255)
          cAGGCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGc                                                                                                                                            >  2:1564581/1‑150 (MQ=255)
                aTGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTaaa                                                                                                                                      <  1:4303800/150‑1 (MQ=255)
                 tGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAAt                                                                                                                                     <  1:331818/150‑1 (MQ=255)
                  ggAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATc                                                                                                                                    <  1:49221/150‑1 (MQ=255)
                   gAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCt                                                                                                                                   <  1:1272256/150‑1 (MQ=255)
                   gAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCt                                                                                                                                   <  2:3965674/150‑1 (MQ=255)
                   gAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCt                                                                                                                                   <  2:2816415/150‑1 (MQ=255)
                       tACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAgg                                                                                                                               <  1:2690221/150‑1 (MQ=255)
                           cAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCga                                                                                                                           <  1:4570571/150‑1 (MQ=255)
                           cAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCga                                                                                                                           <  1:3211703/150‑1 (MQ=255)
                            aGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCgag                                                                                                                          >  1:3236261/1‑150 (MQ=255)
                               cAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATc                                                                                                                       <  1:1710421/150‑1 (MQ=255)
                                aGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATcc                                                                                                                      >  2:3453953/1‑150 (MQ=255)
                                 ggTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATccc                                                                                                                     >  1:4144484/1‑150 (MQ=255)
                                  gTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCt                                                                                                                    >  2:1095556/1‑150 (MQ=255)
                                       gCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGcc                                                                                                               <  2:2463047/150‑1 (MQ=255)
                                       gCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGcc                                                                                                               <  2:2610871/150‑1 (MQ=255)
                                        cccATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCa                                                                                                              >  1:4268422/1‑150 (MQ=255)
                                         ccATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCAt                                                                                                             <  1:2789366/150‑1 (MQ=255)
                                          cATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATc                                                                                                            <  2:4429953/150‑1 (MQ=255)
                                          cATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATc                                                                                                            >  2:3372344/1‑150 (MQ=255)
                                           aTCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCt                                                                                                           >  2:651291/1‑150 (MQ=255)
                                                  cAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGcc                                                                                                    >  2:842510/1‑150 (MQ=255)
                                                   aaaTAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCa                                                                                                   <  1:867247/150‑1 (MQ=255)
                                                   aaaTAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCa                                                                                                   >  2:911576/1‑150 (MQ=255)
                                                   aaaTAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCa                                                                                                   >  2:3118447/1‑150 (MQ=255)
                                                    aaTAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAg                                                                                                  <  1:422581/150‑1 (MQ=255)
                                                     aTAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAgg                                                                                                 >  1:2960762/1‑150 (MQ=255)
                                                     aTAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAgg                                                                                                 >  2:3876058/1‑150 (MQ=255)
                                                                  ttCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTg                                                                                    <  2:2737045/150‑1 (MQ=255)
                                                                   tCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGa                                                                                   <  2:2989079/150‑1 (MQ=255)
                                                                     cGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGAtt                                                                                 >  1:2052408/1‑150 (MQ=255)
                                                                      gggTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTg                                                                                <  1:637467/150‑1 (MQ=255)
                                                                      gggTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTg                                                                                <  2:360422/150‑1 (MQ=255)
                                                                       ggTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTgg                                                                               >  1:3406261/1‑150 (MQ=255)
                                                                        gTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGc                                                                              >  1:2592173/1‑150 (MQ=255)
                                                                         tAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCa                                                                             <  2:1961312/150‑1 (MQ=255)
                                                                         tAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCa                                                                             >  1:1240217/1‑150 (MQ=255)
                                                                           aCGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGc                                                                           >  2:4487757/1‑150 (MQ=255)
                                                                            cGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCt                                                                          <  1:1193485/150‑1 (MQ=255)
                                                                             ggTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTg                                                                         >  1:1350827/1‑150 (MQ=255)
                                                                               tGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGtt                                                                       >  1:3813186/1‑150 (MQ=255)
                                                                                gATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTg                                                                      <  1:2153316/150‑1 (MQ=255)
                                                                                 aTAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTgc                                                                     <  1:334650/150‑1 (MQ=255)
                                                                                     aaaTCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCa                                                                 >  1:3425449/1‑150 (MQ=255)
                                                                                             ggCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCTTGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCTTCTTCCTGGCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAg                                                         <  1:1652447/150‑1 (MQ=255)
                                                                                             ggCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAg                                                         <  2:1339235/150‑1 (MQ=255)
                                                                                              gCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGa                                                        <  1:3538114/150‑1 (MQ=255)
                                                                                                 aCATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAAt                                                     <  1:1411220/150‑1 (MQ=255)
                                                                                                 aCATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAAt                                                     <  2:3587858/150‑1 (MQ=255)
                                                                                                 aCATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAAt                                                     <  2:1545533/150‑1 (MQ=255)
                                                                                                  cATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATg                                                    >  1:90537/1‑150 (MQ=255)
                                                                                                  cATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATg                                                    >  2:1285189/1‑150 (MQ=255)
                                                                                                   aTCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGt                                                   <  1:544107/150‑1 (MQ=255)
                                                                                                    tCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGtt                                                  <  2:24091/150‑1 (MQ=255)
                                                                                                    tCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGtt                                                  <  1:2419056/150‑1 (MQ=255)
                                                                                                     cTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGttt                                                 >  1:3119997/1‑150 (MQ=255)
                                                                                                       ttGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCt                                               <  2:1851189/150‑1 (MQ=255)
                                                                                                       ttGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCt                                               <  1:2268212/150‑1 (MQ=255)
                                                                                                         gATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTgg                                             <  1:2536051/150‑1 (MQ=255)
                                                                                                         gATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTgg                                             >  1:1745044/1‑150 (MQ=255)
                                                                                                            atTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCCTGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGaaa                                          >  2:2119786/1‑150 (MQ=255)
                                                                                                             ttCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAAt                                         <  1:946440/150‑1 (MQ=255)
                                                                                                              tctcTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATg                                        <  2:4138215/150‑1 (MQ=255)
                                                                                                               ctctTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGc                                       <  1:4396475/150‑1 (MQ=255)
                                                                                                                 ctTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCtt                                     >  2:2459144/1‑150 (MQ=255)
                                                                                                                 ctTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCtt                                     <  2:1841793/150‑1 (MQ=255)
                                                                                                                   ttCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTc                                   <  2:4265456/150‑1 (MQ=255)
                                                                                                                    tCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCa                                  <  1:2259174/150‑1 (MQ=255)
                                                                                                                     ccccAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCAt                                 >  2:2088020/1‑150 (MQ=255)
                                                                                                                        cAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAAt                              <  2:1458788/150‑1 (MQ=255)
                                                                                                                         aGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAAtt                             >  2:2811864/1‑150 (MQ=255)
                                                                                                                         aGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAAtt                             >  1:4334913/1‑150 (MQ=255)
                                                                                                                          ggATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTa                            <  1:1121296/150‑1 (MQ=255)
                                                                                                                             tGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACt                         <  2:3015838/150‑1 (MQ=255)
                                                                                                                             tGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACt                         <  2:4441939/150‑1 (MQ=255)
                                                                                                                             tGCGCGGGTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACt                         <  1:736927/150‑1 (MQ=255)
                                                                                                                              gcgcGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTc                        <  2:1522308/150‑1 (MQ=255)
                                                                                                                                gcgGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCt                      <  1:4087917/150‑1 (MQ=255)
                                                                                                                                 cgGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCtt                     <  1:2380363/150‑1 (MQ=255)
                                                                                                                                   gTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCtttt                   <  1:1548030/150‑1 (MQ=255)
                                                                                                                                   gTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCtttt                   >  1:2662481/1‑150 (MQ=255)
                                                                                                                                   gTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCtttt                   <  1:1327171/150‑1 (MQ=255)
                                                                                                                                   gTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCtttt                   >  1:3609352/1‑150 (MQ=255)
                                                                                                                                    ttCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTa                  <  2:1237557/150‑1 (MQ=255)
                                                                                                                                          ttttCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTa            <  2:1113373/150‑1 (MQ=255)
                                                                                                                                            ttCCTATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGc          <  1:4007410/150‑1 (MQ=255)
                                                                                                                                                 aTTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAatg     >  1:3587758/1‑150 (MQ=255)
                                                                                                                                                   tACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCGTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAatgat   <  2:4351179/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AATTCGCGCCCAGGCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATG  >  CP009273/3896795‑3897092

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: