Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 54,392 G→T 94.5% A104E (GCG→GAG)  surA ← peptidyl‑prolyl cis‑trans isomerase (PPIase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927354,3920GT94.5% 527.5 / 20.6 181A104E (GCG→GAG) surApeptidyl‑prolyl cis‑trans isomerase (PPIase)
Reads supporting (aligned to +/- strand):  ref base G (3/7);  new base T (90/81);  total (93/88)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 181 total aligned reads displayed.

TGGTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCT  >  CP009273/54249‑54534
                                                                                                                                               |                                                                                                                                              
tGGTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGAt                                                                                                                                            <  2:2133189/148‑1 (MQ=255)
tGGTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGAt                                                                                                                                            >  1:298964/1‑148 (MQ=255)
 ggTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTGTCTCGATg                                                                                                                                           <  2:68202/148‑1 (MQ=255)
   tGATGCGACGACGCACCTCGGTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGtt                                                                                                                                         <  1:3145399/148‑1 (MQ=255)
     aTGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAg                                                                                                                                       <  1:3404627/148‑1 (MQ=255)
             acgCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGcc                                                                                                                               <  2:3352040/148‑1 (MQ=255)
              cgCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCt                                                                                                                              <  1:1199802/148‑1 (MQ=255)
              cgCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCt                                                                                                                              >  1:2195011/1‑148 (MQ=255)
                cACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGa                                                                                                                            <  2:977361/148‑1 (MQ=255)
                cACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGa                                                                                                                            <  2:2600857/148‑1 (MQ=255)
                  ccTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATc                                                                                                                          >  2:305074/1‑148 (MQ=255)
                  ccTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATc                                                                                                                          >  2:2760498/1‑148 (MQ=255)
                   cTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATcc                                                                                                                         >  2:2577264/1‑148 (MQ=255)
                   cTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATcc                                                                                                                         >  2:2564349/1‑148 (MQ=255)
                    tCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCa                                                                                                                        >  1:2613862/1‑148 (MQ=255)
                     cGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAg                                                                                                                       <  2:2693557/148‑1 (MQ=255)
                     cGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAg                                                                                                                       >  1:783658/1‑148 (MQ=255)
                        tgttACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAgctg                                                                                                                    <  1:2686343/148‑1 (MQ=255)
                        tgttACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAgctg                                                                                                                    <  1:2682847/148‑1 (MQ=255)
                          ttACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAgctgct                                                                                                                  >  1:1595803/1‑148 (MQ=255)
                          ttACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAgctgct                                                                                                                  >  1:2817483/1‑148 (MQ=255)
                           tACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTc                                                                                                                 >  1:599975/1‑148 (MQ=255)
                             cGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCAt                                                                                                               <  2:885705/148‑1 (MQ=255)
                             cGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCAt                                                                                                               <  1:2316841/148‑1 (MQ=255)
                               cACTTCAGATATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCg                                                                                                             <  1:2981580/148‑1 (MQ=255)
                                aCTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCgg                                                                                                            <  2:1960432/148‑1 (MQ=255)
                                 cTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGga                                                                                                           >  2:3171410/1‑148 (MQ=255)
                                 cTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGga                                                                                                           <  1:3497652/148‑1 (MQ=255)
                                  ttCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGgag                                                                                                          >  1:707540/1‑148 (MQ=255)
                                  ttCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGgag                                                                                                          >  1:3265857/1‑148 (MQ=255)
                                     agagATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGAtt                                                                                                       <  1:2527524/148‑1 (MQ=255)
                                         aTAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCa                                                                                                   >  2:2319001/1‑148 (MQ=255)
                                         aTAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCa                                                                                                   <  1:110663/148‑1 (MQ=255)
                                         aTAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCa                                                                                                   >  2:400234/1‑148 (MQ=255)
                                          tAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCAc                                                                                                  >  2:3104940/1‑148 (MQ=255)
                                           aaTCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACt                                                                                                 <  1:940746/148‑1 (MQ=255)
                                            atcatcTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTc                                                                                                <  1:3017221/148‑1 (MQ=255)
                                            atcatcTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTc                                                                                                >  1:2112872/1‑148 (MQ=255)
                                            atcatcTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTc                                                                                                >  1:2112772/1‑148 (MQ=255)
                                             tcatcTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTcc                                                                                               <  1:553246/148‑1 (MQ=255)
                                             tcatcTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTcc                                                                                               <  2:1059949/148‑1 (MQ=255)
                                               atcTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCa                                                                                             <  2:932611/148‑1 (MQ=255)
                                                 ctctTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCAtt                                                                                           <  2:1942416/148‑1 (MQ=255)
                                                   ctTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCAtttt                                                                                         <  1:1632254/148‑1 (MQ=255)
                                                      tGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTg                                                                                      <  2:2212488/148‑1 (MQ=255)
                                                         ggATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGccc                                                                                   <  1:1977765/148‑1 (MQ=255)
                                                            tCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCAt                                                                                <  2:623269/148‑1 (MQ=255)
                                                            tCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCAt                                                                                >  2:910334/1‑148 (MQ=255)
                                                             cTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATc                                                                               <  2:303808/148‑1 (MQ=255)
                                                             cTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATc                                                                               >  2:2940778/1‑148 (MQ=255)
                                                             cTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATc                                                                               <  1:1824550/148‑1 (MQ=255)
                                                                   aCGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAgg                                                                         <  1:2356956/148‑1 (MQ=255)
                                                                      aTAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGgatg                                                                      >  2:869920/1‑148 (MQ=255)
                                                                       tAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGgatga                                                                     <  2:1158390/148‑1 (MQ=255)
                                                                        aGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGgatgat                                                                    >  2:2913187/1‑148 (MQ=255)
                                                                          gtgtTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGAttt                                                                  >  2:817982/1‑148 (MQ=255)
                                                                           tgttgtAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTg                                                                 >  1:1499568/1‑148 (MQ=255)
                                                                           tgttgtAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTg                                                                 >  1:1499641/1‑148 (MQ=255)
                                                                             ttgtAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGAt                                                               <  1:2480531/148‑1 (MQ=255)
                                                                                tAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCa                                                            <  1:1338518/148‑1 (MQ=255)
                                                                                  gTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATg                                                          <  2:2188108/148‑1 (MQ=255)
                                                                                   ttCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGa                                                         >  2:716730/1‑148 (MQ=255)
                                                                                   ttCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGa                                                         <  2:1950868/148‑1 (MQ=255)
                                                                                     cAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATc                                                       >  1:1819244/1‑148 (MQ=255)
                                                                                       gTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCaa                                                     >  1:3053439/1‑148 (MQ=255)
                                                                                        tCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGAAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCaaa                                                    >  2:2398966/1‑148 (MQ=255)
                                                                                           aTCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGt                                                 >  1:1044293/1‑148 (MQ=255)
                                                                                             cGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTc                                               >  2:1503985/1‑148 (MQ=255)
                                                                                              gTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTcc                                              <  2:1204419/148‑1 (MQ=255)
                                                                                                aaGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCAt                                            <  2:783658/148‑1 (MQ=255)
                                                                                                aaGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCAt                                            >  1:989413/1‑148 (MQ=255)
                                                                                                   ccAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGAt                                         >  1:311574/1‑148 (MQ=255)
                                                                                                      gACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTtg                                      <  2:1667628/148‑1 (MQ=255)
                                                                                                            tGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGgcg                                <  1:3545546/148‑1 (MQ=255)
                                                                                                                cATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgc                            >  1:2627982/1‑148 (MQ=255)
                                                                                                                cATCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgc                            >  1:3410116/1‑148 (MQ=255)
                                                                                                                 aTCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcg                           <  2:994724/148‑1 (MQ=255)
                                                                                                                 aTCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcg                           <  1:2940778/148‑1 (MQ=255)
                                                                                                                  tCTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgt                          >  2:3283109/1‑148 (MQ=255)
                                                                                                                   cTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgtc                         <  1:61884/148‑1 (MQ=255)
                                                                                                                   cTGATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgtc                         >  1:1909939/1‑148 (MQ=255)
                                                                                                                   cTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgtc                         >  2:826989/1‑148 (MQ=255)
                                                                                                                     gATCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgtcgc                       >  1:2888100/1‑148 (MQ=255)
                                                                                                                      aTCCAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgtcgcg                      >  1:2586679/1‑148 (MQ=255)
                                                                                                                        ccAGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAgcgtcgcgtc                    >  2:590360/1‑148 (MQ=255)
                                                                                                                          aGCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGtcat                  >  1:405625/1‑148 (MQ=255)
                                                                                                                           gCGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGtcatc                 >  2:3458793/1‑148 (MQ=255)
                                                                                                                            cGTCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATTATTTGGTGGCGCAGCGTCGCGtcatca                >  1:2839352/1‑148 (MQ=255)
                                                                                                                              tCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAgg              <  1:496908/148‑1 (MQ=255)
                                                                                                                              tCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAgg              >  1:154361/1‑148 (MQ=255)
                                                                                                                              tCATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAgg              >  1:2463084/1‑148 (MQ=255)
                                                                                                                               cATGTTGTTCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGa             >  2:1007639/1‑148 (MQ=255)
                                                                                                                                  gttgttCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGt          >  2:1478395/1‑148 (MQ=255)
                                                                                                                                  gttgttCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGt          >  1:1984881/1‑148 (MQ=255)
                                                                                                                                  gttgttCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGt          >  2:2210553/1‑148 (MQ=255)
                                                                                                                                  gttgttCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTCGGTGGCGCAGCGTCGCGTCATCAGGAAGt          >  1:3020867/1‑148 (MQ=255)
                                                                                                                                       tCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTtgctg     >  2:3394956/1‑148 (MQ=255)
                                                                                                                                       tCTGTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTtgctg     >  1:1028156/1‑148 (MQ=255)
                                                                                                                                          gTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCtgcct  <  2:1316972/148‑1 (MQ=255)
                                                                                                                                          gTTTCTCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCtgcct  <  2:1316952/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                              
TGGTGATGCGACGACGCACCTCGTTGTTACGCACTTCAGAGATAATCATCTCTTTGCGGATCTGGTTACGATAGGTGTTGTAGTTCAGTCCATCGTAAGCCAGACGGCTGCGCATCTGATCCAGCGTCATGTTGTTCTGTTTCGCGATGTTAGCAATCGCCTGATCCAGCTGCTCATCGGAGATTTTCACTCCCATTTTCTGCCCCATCTGCAGGATGATTTGATCCATGATCAAACGTTCCATGATTTGGTGGCGCAGCGTCGCGTCATCAGGAAGTTGCTGCCT  >  CP009273/54249‑54534

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: ‑