Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 667,899 C→A 100% D781Y (GAT→TAT)  leuS ← leucyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273667,8990CA96.7% 927.6 / 18.5 302D781Y (GAT→TAT) leuSleucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (3/7);  new base A (142/150);  total (145/157)
Fisher's exact test for biased strand distribution p-value = 3.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 323 total aligned reads displayed.

GCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGTTTGTTCATCA  >  CP009273/667750‑668045
                                                                                                                                                     |                                                                                                                                                  
     cACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATAtagcctt                                                                                                                                               <  2:1898383/150‑1 (MQ=255)
     cACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATAtagcctt                                                                                                                                               >  2:1225222/1‑150 (MQ=255)
       cGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATAtagccttcg                                                                                                                                             <  2:2689819/150‑1 (MQ=255)
        gTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATAtagccttcgc                                                                                                                                            >  2:3343936/1‑150 (MQ=255)
        gTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATAtagccttcgc                                                                                                                                            >  2:4193515/1‑150 (MQ=255)
        gTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATAtagccttcgc                                                                                                                                            >  1:4114934/1‑150 (MQ=255)
        gTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATAtagccttcgc                                                                                                                                            >  1:2310995/1‑150 (MQ=255)
         ttCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATAtagccttcgcc                                                                                                                                           >  1:2830946/1‑150 (MQ=255)
           cgcgAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATagccttcgcctt                                                                                                                                         <  2:2559260/150‑1 (MQ=255)
             cgAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTc                                                                                                                                       >  2:3900688/1‑150 (MQ=255)
             cgAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTc                                                                                                                                       <  1:4452221/150‑1 (MQ=255)
                aCCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGt                                                                                                                                    >  1:4856242/1‑150 (MQ=255)
                  cTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTTGATATAGCCTTCGCCTTTCAGTTc                                                                                                                                  >  2:3943277/1‑150 (MQ=255)
                  cTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTc                                                                                                                                  >  2:2913403/1‑150 (MQ=255)
                   tGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTcc                                                                                                                                 >  1:42618/1‑150 (MQ=255)
                     ttcttcCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTg                                                                                                                               <  1:3489653/150‑1 (MQ=255)
                     ttcttcCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTg                                                                                                                               >  1:171426/1‑150 (MQ=255)
                       cttcCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGcc                                                                                                                             >  1:1518958/1‑150 (MQ=255)
                           cGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAg                                                                                                                         <  2:1472046/150‑1 (MQ=255)
                           cGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAg                                                                                                                         >  2:2462926/1‑150 (MQ=255)
                             ttGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCg                                                                                                                       >  2:2312085/1‑150 (MQ=255)
                              tGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGt                                                                                                                      <  2:4114617/150‑1 (MQ=255)
                                   ccACCGGAACGGTGATTTGGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGc                                                                                                                 <  2:2744664/150‑1 (MQ=255)
                                    cACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCa                                                                                                                <  1:4654503/150‑1 (MQ=255)
                                    cACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCa                                                                                                                <  2:4767608/150‑1 (MQ=255)
                                       cGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGAt                                                                                                             >  1:4885362/1‑150 (MQ=255)
                                        ggAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACGGCGTGAAGCAGagg                                                                                                            >  2:3027611/1‑148 (MQ=255)
                                        ggAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGAtg                                                                                                            >  1:3733123/1‑150 (MQ=255)
                                        ggAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGAtg                                                                                                            >  1:146472/1‑150 (MQ=255)
                                          aaCGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGAtgtg                                                                                                          <  2:1670190/150‑1 (MQ=255)
                                             ggTGTTTTTGGCACGGACTTTACCGTTAACCTGCACCTCGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGAAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGATGGAGATGTGCgg                                                                                                       <  2:1520351/150‑1 (MQ=255)
                                             ggTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCgg                                                                                                       <  1:2695460/150‑1 (MQ=255)
                                               tGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCgggg                                                                                                     >  1:3542596/1‑150 (MQ=255)
                                                gATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGt                                                                                                    <  2:5132072/150‑1 (MQ=255)
                                                gATTTTGGCACGGACTTGACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGt                                                                                                    <  2:3662265/150‑1 (MQ=255)
                                                 aTTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTg                                                                                                   >  1:3770355/1‑150 (MQ=255)
                                                  ttttGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGa                                                                                                  <  1:2486937/150‑1 (MQ=255)
                                                     tGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGACCg                                                                                               >  2:2031951/1‑150 (MQ=255)
                                                      ggCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACgg                                                                                              <  1:3338507/150‑1 (MQ=255)
                                                       gCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCCGCCACAGCGTGAAGCAGATGTGCGGGGTGAACggg                                                                                             <  1:1102266/150‑1 (MQ=255)
                                                            gACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAg                                                                                        <  1:3113465/150‑1 (MQ=255)
                                                            gACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAg                                                                                        <  1:4089665/150‑1 (MQ=255)
                                                              cTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCa                                                                                      >  1:3953241/1‑150 (MQ=255)
                                                                 tACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATAc                                                                                   <  1:2047778/150‑1 (MQ=255)
                                                                    cGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGa                                                                                >  1:480746/1‑150 (MQ=255)
                                                                    cGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGa                                                                                <  1:5156833/150‑1 (MQ=255)
                                                                    cGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGa                                                                                <  1:113947/150‑1 (MQ=255)
                                                                     gTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGAc                                                                               >  1:799590/1‑150 (MQ=255)
                                                                      ttAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACa                                                                              <  1:1279462/150‑1 (MQ=255)
                                                                        aaCCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAAc                                                                            <  1:5161806/150‑1 (MQ=255)
                                                                        aaCCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAAc                                                                            >  1:1693794/1‑150 (MQ=255)
                                                                            tGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGcc                                                                        >  2:1031362/1‑150 (MQ=255)
                                                                             gCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGACGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCca                                                                       >  2:5066143/1‑150 (MQ=255)
                                                                             gCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCca                                                                       <  2:239543/150‑1 (MQ=255)
                                                                              caccacGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCcag                                                                      >  1:1831105/1‑150 (MQ=255)
                                                                                ccacGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCcagca                                                                    <  1:1961133/150‑1 (MQ=255)
                                                                                ccacGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCcagca                                                                    >  1:4304055/1‑150 (MQ=255)
                                                                                   cGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTg                                                                 <  1:3293910/150‑1 (MQ=255)
                                                                                       cAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTc                                                             <  1:933684/150‑1 (MQ=255)
                                                                                        aGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCCTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTcc                                                            >  2:3443477/1‑150 (MQ=255)
                                                                                         gCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCt                                                           >  2:4475643/1‑150 (MQ=255)
                                                                                           gTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGc                                                         <  1:3035989/150‑1 (MQ=255)
                                                                                            tGGAGTCTTCCACCAGCGCTTTTTCGTCAGCAACCGGCCCCGGCGCGGTGTCGATATAGCCTTCGCCTTTCCGTTCCTGCCACAGCGTGACGCAGATGTGCGGGGTGAACGGGTTAAGCAGACGGACAACGGCCAGCAGTGCTTCCGGCa                                                        >  2:2289220/1‑150 (MQ=255)
                                                                                               aGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCa                                                     >  2:3187474/1‑150 (MQ=255)
                                                                                               aGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCa                                                     <  2:3317751/150‑1 (MQ=255)
                                                                                                    tCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAgcg                                                >  2:4529156/1‑150 (MQ=255)
                                                                                                    tCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAgcg                                                <  1:2462926/150‑1 (MQ=255)
                                                                                                      caccaTCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGcgcg                                              >  1:3835093/1‑150 (MQ=255)
                                                                                                      caccaTCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGcgcg                                              >  2:3502315/1‑150 (MQ=255)
                                                                                                           tCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCt                                         <  1:2749224/150‑1 (MQ=255)
                                                                                                           tCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCt                                         <  1:630727/150‑1 (MQ=255)
                                                                                                                   tCGTCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCGTCCTGCATCAGAGCGCGATCCTGCTCGccc                                 >  1:2276607/1‑149 (MQ=255)
                                                                                                                    cttCAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCAt                                >  1:1621470/3‑150 (MQ=255)
                                                                                                                       cAGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCgg                             >  1:4977398/1‑150 (MQ=255)
                                                                                                                        aGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGt                            >  1:3673982/1‑150 (MQ=255)
                                                                                                                        aGCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGt                            <  2:2132450/150‑1 (MQ=255)
                                                                                                                         gCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGtt                           <  1:4488947/150‑1 (MQ=255)
                                                                                                                         gCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGtt                           <  1:996512/150‑1 (MQ=255)
                                                                                                                         gCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGtt                           <  2:5004316/150‑1 (MQ=255)
                                                                                                                         gCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGtt                           <  1:100736/150‑1 (MQ=255)
                                                                                                                         gCAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGtt                           <  1:1031362/150‑1 (MQ=255)
                                                                                                                          cAACCGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTg                          >  2:931235/1‑150 (MQ=255)
                                                                                                                             ccGGCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTg                       <  2:4346330/150‑1 (MQ=255)
                                                                                                                               ggCCACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCt                     >  1:689107/1‑150 (MQ=255)
                                                                                                                                 ccACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCttt                   >  1:4534178/1‑150 (MQ=255)
                                                                                                                                 ccACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCttt                   <  2:4671956/150‑1 (MQ=255)
                                                                                                                                  cACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTc                  <  1:3162460/150‑1 (MQ=255)
                                                                                                                                  cACGGCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTc                  <  2:4856242/150‑1 (MQ=255)
                                                                                                                                     ggCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATCCGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGcc               >  2:665466/1‑150 (MQ=255)
                                                                                                                                     ggCGCGTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCCGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGcc               >  2:3251307/1‑150 (MQ=255)
                                                                                                                                       cgcgTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAg             >  1:609333/1‑150 (MQ=255)
                                                                                                                                       cgcgTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAg             >  1:2576145/1‑150 (MQ=255)
                                                                                                                                          gTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGttt          <  1:2913403/150‑1 (MQ=255)
                                                                                                                                          gTTGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGttt          >  2:724052/1‑150 (MQ=255)
                                                                                                                                           ttGTCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGTTTg         <  2:845655/150‑1 (MQ=255)
                                                                                                                                              tCGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGTTTGTtc      >  2:4289426/1‑150 (MQ=255)
                                                                                                                                               cGATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGTTTGTtca     >  2:4090110/1‑150 (MQ=255)
                                                                                                                                               agATATAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGTTTGTtca     >  2:3614583/2‑150 (MQ=255)
                                                                                                                                                 atatAGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGTTTGTtcatc   >  1:266028/1‑150 (MQ=255)
                                                                                                                                                  tataGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGTTTGTtcatca  >  1:2029580/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCCAGCACGTTCGCGAACCTGTTCTTCCGTTGCGTCCACCGGAACGGTGATTTTGGCACGGACTTTACCGTTAACCTGCACCACGACCAGCGTGGAGTCTTCCACCATCGCTTTTTCGTCAGCAACCGGCCACGGCGCGTTGTCGATATCGCCTTCGCCTTTCAGTTCCTGCCACAGCGTGAAGCAGATGTGCGGGGTGAACGGGTTAAGCATACGGACAACGGCCAGCAGTGCTTCCTGCATCAGAGCGCGATCCTGCTCGCCATCGGTTGGTGCTTTCGCCAGTTTGTTCATCA  >  CP009273/667750‑668045

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: