Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,896,944 C→A 15.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,9440CA15.3% 961.7 / 103.0 293V40L (GTA→TTA) bglBcryptic phospho‑beta‑glucosidase B
Reads supporting (aligned to +/- strand):  ref base C (94/154);  new base A (22/23);  total (116/177)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 296 total aligned reads displayed.

AATTCGCGCCCAGGCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGG  >  CP009273/3896795‑3897093
                                                                                                                                                     |                                                                                                                                                     
 aTTCGCGCCCAGGCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACg                                                                                                                                                      >  1:4311484/1‑150 (MQ=255)
         ccAGGCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAgg                                                                                                                                              >  1:4498420/1‑150 (MQ=255)
         ccAGGCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAgg                                                                                                                                              >  1:3610385/1‑150 (MQ=255)
           aGGCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCt                                                                                                                                            <  2:926390/150‑1 (MQ=255)
           aGGCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCt                                                                                                                                            <  1:79838/150‑1 (MQ=255)
           aGGCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCt                                                                                                                                            <  1:1880575/150‑1 (MQ=255)
            ggCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTg                                                                                                                                           <  1:4112472/150‑1 (MQ=255)
             gCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGt                                                                                                                                          >  2:3682962/1‑150 (MQ=255)
              cAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTa                                                                                                                                         >  1:3106801/1‑150 (MQ=255)
                         cGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTc                                                                                                                              >  2:1274630/1‑150 (MQ=255)
                            aGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCgag                                                                                                                           >  2:2948615/1‑150 (MQ=255)
                                aGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATcc                                                                                                                       >  1:2984528/1‑150 (MQ=255)
                                   tGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCtt                                                                                                                    >  2:521715/1‑150 (MQ=255)
                                   tGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCtt                                                                                                                    <  2:3170694/150‑1 (MQ=255)
                                    gAAGCCCATCTCGGCAAATCACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCttt                                                                                                                   <  1:3863868/150‑1 (MQ=255)
                                        cccATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCa                                                                                                               >  1:3461424/1‑150 (MQ=255)
                                         ccATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCAt                                                                                                              >  2:4055107/1‑150 (MQ=255)
                                           aTCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCt                                                                                                            <  2:4139569/150‑1 (MQ=255)
                                           aTCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCt                                                                                                            <  2:3712300/150‑1 (MQ=255)
                                            tctcGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCtt                                                                                                           <  2:2182250/150‑1 (MQ=255)
                                             ctcGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTc                                                                                                          >  2:871956/1‑150 (MQ=255)
                                              tcGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTcc                                                                                                         <  2:3975911/150‑1 (MQ=255)
                                              tcGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTcc                                                                                                         <  2:563858/150‑1 (MQ=255)
                                                   aaaTAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCa                                                                                                    >  1:2668596/1‑150 (MQ=255)
                                                   aaaTAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCa                                                                                                    <  2:3931185/150‑1 (MQ=255)
                                                   aaaTAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCa                                                                                                    <  1:2268778/150‑1 (MQ=255)
                                                      tAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGc                                                                                                 >  1:2608818/1‑150 (MQ=255)
                                                       aaCGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCa                                                                                                >  2:4278386/1‑150 (MQ=255)
                                                        aCGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCAc                                                                                               <  1:100673/150‑1 (MQ=255)
                                                         cgcgATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCAcc                                                                                              >  1:1221891/1‑150 (MQ=255)
                                                           cgATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCtt                                                                                            <  2:4310846/150‑1 (MQ=255)
                                                           cgATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCtt                                                                                            >  1:2691863/1‑150 (MQ=255)
                                                            gATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTc                                                                                           >  2:1753208/1‑150 (MQ=255)
                                                              tatCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCa                                                                                         >  1:3289195/1‑150 (MQ=255)
                                                                 cTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCt                                                                                      >  1:2336089/1‑150 (MQ=255)
                                                                    ccGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGAt                                                                                   <  2:2619640/150‑1 (MQ=255)
                                                                      gggTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTg                                                                                 <  1:3554130/150‑1 (MQ=255)
                                                                      gggTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTg                                                                                 <  2:144198/150‑1 (MQ=255)
                                                                           aCGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGc                                                                            <  1:1140576/150‑1 (MQ=255)
                                                                             ggTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTg                                                                          >  2:1477295/1‑150 (MQ=255)
                                                                              gTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGt                                                                         >  1:594872/1‑150 (MQ=255)
                                                                                gATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTg                                                                       >  2:1124240/1‑150 (MQ=255)
                                                                                  tAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTgcg                                                                     <  2:656233/150‑1 (MQ=255)
                                                                                  tAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTgcg                                                                     >  2:3858157/1‑150 (MQ=255)
                                                                                     aaaTCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCa                                                                  <  1:1112643/150‑1 (MQ=255)
                                                                                        tCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCAccc                                                               <  2:758667/150‑1 (MQ=255)
                                                                                         cGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCAcccc                                                              <  1:3476108/150‑1 (MQ=255)
                                                                                         cGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCAcccc                                                              >  2:2566709/1‑150 (MQ=255)
                                                                                         cGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCAcccc                                                              >  1:447374/1‑150 (MQ=255)
                                                                                          gATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCa                                                             >  2:2100225/1‑150 (MQ=255)
                                                                                          gATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCa                                                             >  1:761885/1‑150 (MQ=255)
                                                                                          gATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCa                                                             <  1:4018987/150‑1 (MQ=255)
                                                                                          gATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCa                                                             <  2:3592893/150‑1 (MQ=255)
                                                                                            tGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCaaa                                                           >  2:388893/1‑150 (MQ=255)
                                                                                             ggCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAg                                                          <  1:818632/150‑1 (MQ=255)
                                                                                             ggCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAg                                                          <  1:2946572/150‑1 (MQ=255)
                                                                                               cGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGaa                                                        >  1:3217413/1‑150 (MQ=255)
                                                                                               cGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGaa                                                        <  1:1366635/150‑1 (MQ=255)
                                                                                                gACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGaaa                                                       >  2:1231779/1‑150 (MQ=255)
                                                                                                gACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCAAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGaaa                                                       >  2:2215192/1‑150 (MQ=255)
                                                                                                  cATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGt                                                    >  1:3364953/1‑150 (MQ=255)
                                                                                                  cATCTTTGATATTCTCTTTACCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGaaagg                                                     >  2:2077107/1‑148 (MQ=255)
                                                                                                   aTCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGt                                                    <  1:3692168/150‑1 (MQ=255)
                                                                                                    tCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATAGGCAGCTGTTGAGCCACCCCAAAGAAATGtt                                                   >  2:2502756/1‑150 (MQ=255)
                                                                                                       ttGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCt                                                <  2:2777278/150‑1 (MQ=255)
                                                                                                       ttGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCt                                                <  2:1649822/150‑1 (MQ=255)
                                                                                                       ttGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCt                                                <  1:3974622/150‑1 (MQ=255)
                                                                                                               ctctTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCGGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGc                                        >  1:1604420/1‑150 (MQ=255)
                                                                                                                tctTTCCCCAGGATGCGCGGTTCCTTTTTTCCCATTACGACCTGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCt                                       <  2:1185092/150‑1 (MQ=255)
                                                                                                                tctTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCt                                       <  2:252712/150‑1 (MQ=255)
                                                                                                                  tttCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCttt                                     >  2:3023384/1‑150 (MQ=255)
                                                                                                                  tttCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCttt                                     <  1:2771201/150‑1 (MQ=255)
                                                                                                                   ttCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTc                                    <  2:4465935/150‑1 (MQ=255)
                                                                                                                      cccAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATa                                 <  2:426241/150‑1 (MQ=255)
                                                                                                                      accAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATa                                 <  2:4066663/149‑1 (MQ=255)
                                                                                                                       ccAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATaa                                >  1:1834097/1‑150 (MQ=255)
                                                                                                                       ccAGGATGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATaa                                >  1:3051288/1‑150 (MQ=255)
                                                                                                                         aGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAAtt                              >  2:2753680/1‑150 (MQ=255)
                                                                                                                         aGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAAtt                              <  1:2456271/150‑1 (MQ=255)
                                                                                                                          ggATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTa                             <  2:4549033/150‑1 (MQ=255)
                                                                                                                            aTGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAAc                           <  2:2026693/150‑1 (MQ=255)
                                                                                                                            aTGCGCGGTTCCATTTTTCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAAc                           <  1:4073409/150‑1 (MQ=255)
                                                                                                                                gcgGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCt                       >  1:3836813/1‑150 (MQ=255)
                                                                                                                                gcgGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCt                       >  1:1511235/1‑150 (MQ=255)
                                                                                                                                 cgGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCtt                      <  1:3106063/150‑1 (MQ=255)
                                                                                                                                  ggTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCttt                     <  1:3732734/150‑1 (MQ=255)
                                                                                                                                  ggTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCttt                     >  2:1525331/1‑150 (MQ=255)
                                                                                                                                   gTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCtttt                    <  1:1948720/150‑1 (MQ=255)
                                                                                                                                   gTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCtttt                    >  1:1497980/1‑150 (MQ=255)
                                                                                                                                      ccATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATc                 <  1:871956/150‑1 (MQ=255)
                                                                                                                                       cATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCg                >  2:4636472/1‑150 (MQ=255)
                                                                                                                                         tttttCCCATTACGACATGAGGCTGTAAATCTGAGGTCGAGATACCTTTGCCATCTTCCTGCAAGGGACCTTCCACCTGATTGGCATCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATATTTAAATCCTTTTATCGtt              <  2:2783570/150‑1 (MQ=255)
                                                                                                                                          ttttCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTa             <  1:2110313/150‑1 (MQ=255)
                                                                                                                                          ttttCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTa             <  2:334497/150‑1 (MQ=255)
                                                                                                                                          ttttCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTa             <  2:1404317/150‑1 (MQ=255)
                                                                                                                                             tCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCg          >  1:568863/1‑150 (MQ=255)
                                                                                                                                             tCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCg          <  1:1569435/150‑1 (MQ=255)
                                                                                                                                             tCCCATTAAGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCg          >  2:3520776/1‑150 (MQ=255)
                                                                                                                                                cATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAat       <  2:1738893/150‑1 (MQ=255)
                                                                                                                                                    aCGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAatgatg   <  1:453736/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATTCGCGCCCAGGCAATGGAAATACGCAGACAGGTGAAGCCCATCTCGGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGG  >  CP009273/3896795‑3897093

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: