Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,808,907 C→A 7.5% E298* (GAA→TAA)  gshA ← glutamate‑cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,808,9070CA 7.5% 1060.3 / 38.9 280E298* (GAA→TAA) gshAglutamate‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base C (128/131);  new base A (14/7);  total (142/138)
Fisher's exact test for biased strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 283 total aligned reads displayed.

CACCTGCTGTTCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGAT  >  CP009273/2808758‑2809056
                                                                                                                                                     |                                                                                                                                                     
  ccTGCTGTTCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGt                                                                                                                                                     <  1:4230423/150‑1 (MQ=255)
    tgctgTTCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGttt                                                                                                                                                   <  2:4019863/150‑1 (MQ=255)
     gctgTTCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGtttt                                                                                                                                                  <  1:3389967/150‑1 (MQ=255)
     gctgTTCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGtttt                                                                                                                                                  >  2:2753342/1‑150 (MQ=255)
     gctgTTCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGtttt                                                                                                                                                  <  2:3147807/150‑1 (MQ=255)
        gTTCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCaa                                                                                                                                               >  2:1810042/1‑150 (MQ=255)
         ttCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAAt                                                                                                                                              <  1:2784655/150‑1 (MQ=255)
            aTCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTg                                                                                                                                           >  1:1294365/1‑150 (MQ=255)
               tACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCaa                                                                                                                                        >  2:490885/1‑150 (MQ=255)
               tACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCaa                                                                                                                                        >  2:2121946/1‑150 (MQ=255)
                 cacCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAaca                                                                                                                                      >  1:3521318/1‑150 (MQ=255)
                 cacCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAaca                                                                                                                                      >  1:918184/1‑150 (MQ=255)
                 cacCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAaca                                                                                                                                      >  2:829823/1‑150 (MQ=255)
                  acCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAacac                                                                                                                                     >  1:4526842/1‑150 (MQ=255)
                        cGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCt                                                                                                                               <  2:2199226/150‑1 (MQ=255)
                           cGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGtt                                                                                                                            >  1:1970418/1‑150 (MQ=255)
                            gagaACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTg                                                                                                                           <  1:1055801/150‑1 (MQ=255)
                             agaACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGa                                                                                                                          >  2:3536043/1‑150 (MQ=255)
                               aaCGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGAtt                                                                                                                        >  2:4234389/1‑150 (MQ=255)
                               aaCGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTAGCTGTTGAtt                                                                                                                        >  2:4070142/1‑150 (MQ=255)
                                aCGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGAttt                                                                                                                       <  2:4188876/150‑1 (MQ=255)
                                  gggTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGc                                                                                                                     >  2:2525481/1‑150 (MQ=255)
                                   ggTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCa                                                                                                                    <  1:2488850/150‑1 (MQ=255)
                                         tGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCctct                                                                                                              <  2:1669917/150‑1 (MQ=255)
                                               gCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGt                                                                                                        <  1:137156/150‑1 (MQ=255)
                                                 gAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCt                                                                                                      <  2:1196029/150‑1 (MQ=255)
                                                 gAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCt                                                                                                      >  1:4573973/1‑150 (MQ=255)
                                                 gAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCt                                                                                                      >  2:901594/1‑150 (MQ=255)
                                                  aaCGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCtt                                                                                                     >  2:1158645/1‑150 (MQ=255)
                                                  aaCGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCtt                                                                                                     >  2:519669/1‑150 (MQ=255)
                                                   aCGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCttt                                                                                                    <  1:3217827/150‑1 (MQ=255)
                                                     gCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTATTACCGTCTTtct                                                                                                  >  2:4321504/1‑150 (MQ=255)
                                                       aCTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCa                                                                                                >  1:4384169/1‑150 (MQ=255)
                                                        cTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCaa                                                                                               >  1:2756509/1‑150 (MQ=255)
                                                          tCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATa                                                                                             <  1:3239329/150‑1 (MQ=255)
                                                            aaTATTTTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGTATCGGCGCGTACAGTTCGTTTTCATTCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATAcc                                                                                           <  1:752053/150‑1 (MQ=255)
                                                            aaTATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATAcc                                                                                           >  1:3904965/1‑150 (MQ=255)
                                                              tatatTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCaa                                                                                         >  1:1637701/1‑150 (MQ=255)
                                                                tatTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATc                                                                                       <  1:259405/150‑1 (MQ=255)
                                                                tatTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATc                                                                                       >  1:1861454/1‑150 (MQ=255)
                                                                 atTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCt                                                                                      <  2:1407373/150‑1 (MQ=255)
                                                                 atTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCt                                                                                      >  2:1361686/1‑150 (MQ=255)
                                                                  ttCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCtt                                                                                     <  2:1944133/150‑1 (MQ=255)
                                                                  ttCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCtt                                                                                     >  2:463574/1‑150 (MQ=255)
                                                                   tCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTc                                                                                    >  2:3709875/1‑150 (MQ=255)
                                                                    cAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTcg                                                                                   >  2:4160093/1‑150 (MQ=255)
                                                                      aTGCCGCCACGTAACAGCGCATCATAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTATCAATACCAATCTTcgcg                                                                                 <  2:2911622/150‑1 (MQ=255)
                                                                      aTGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTcgcg                                                                                 <  1:2899749/150‑1 (MQ=255)
                                                                        gccgccACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTa                                                                               >  2:3020605/1‑150 (MQ=255)
                                                                          cgccACGTAACAGCGCATCAGAAGGCGACTCGTCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTAct                                                                             >  2:3620157/1‑150 (MQ=255)
                                                                          cgccACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTAct                                                                             >  1:1552446/1‑150 (MQ=255)
                                                                           gccACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTActc                                                                            >  1:1718654/1‑150 (MQ=255)
                                                                           gccACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTActc                                                                            >  1:4349453/1‑150 (MQ=255)
                                                                             cACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCtt                                                                          <  2:1093272/150‑1 (MQ=255)
                                                                             cACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCtt                                                                          <  2:30623/150‑1 (MQ=255)
                                                                                gTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCg                                                                       <  2:4457928/150‑1 (MQ=255)
                                                                                gTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTAACGTCTTTCTCAATACCAATCTACGCGTACTCTTcct                                                                       >  1:226187/1‑149 (MQ=255)
                                                                                 tAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTccta                                                                      <  2:2684071/150‑3 (MQ=255)
                                                                                  aaCAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTAt                                                                     >  1:3162153/1‑150 (MQ=255)
                                                                                     aGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGc                                                                  >  2:3247360/1‑150 (MQ=255)
                                                                                        gcATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGtt                                                               >  2:65346/1‑150 (MQ=255)
                                                                                             aGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGAt                                                          <  1:4486698/150‑1 (MQ=255)
                                                                                              gAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGAtt                                                         >  2:2485754/1‑150 (MQ=255)
                                                                                              gAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGAtt                                                         <  1:267310/150‑1 (MQ=255)
                                                                                               aaGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTg                                                        <  2:1880676/150‑1 (MQ=255)
                                                                                               aaGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTg                                                        <  1:4354592/150‑1 (MQ=255)
                                                                                                aGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGc                                                       >  1:1147300/1‑150 (MQ=255)
                                                                                                   cGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTg                                                    <  1:4397223/150‑1 (MQ=255)
                                                                                                    gACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGt                                                   >  1:3629738/1‑150 (MQ=255)
                                                                                                      cTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGttt                                                 >  1:1490910/1‑150 (MQ=255)
                                                                                                       tCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTa                                                >  1:4059842/1‑150 (MQ=255)
                                                                                                        cgccgcTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTaa                                               <  2:3347429/150‑1 (MQ=255)
                                                                                                         gccgcTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAg                                              <  1:2229665/150‑1 (MQ=255)
                                                                                                            gcTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGcc                                           <  1:768165/150‑1 (MQ=255)
                                                                                                             cTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGccc                                          <  1:2365624/150‑1 (MQ=255)
                                                                                                             cTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGccc                                          >  2:2271714/1‑150 (MQ=255)
                                                                                                                 gggTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTa                                      <  1:2483602/150‑1 (MQ=255)
                                                                                                                     aaCGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTa                                  >  2:1099516/1‑150 (MQ=255)
                                                                                                                              tGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAg                         >  2:1303708/1‑150 (MQ=255)
                                                                                                                                gACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATTGCTTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGAt                       <  1:4290121/150‑1 (MQ=255)
                                                                                                                                 aCGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATc                      <  2:1575587/150‑1 (MQ=255)
                                                                                                                                 aCGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATc                      >  2:1999304/1‑150 (MQ=255)
                                                                                                                                   gAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGt                    <  1:1311442/150‑1 (MQ=255)
                                                                                                                                    aaTCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGtt                   <  2:3287250/150‑1 (MQ=255)
                                                                                                                                    aaTCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGtt                   <  2:3097102/150‑1 (MQ=255)
                                                                                                                                    aaTCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGtt                   >  2:3871520/1‑150 (MQ=255)
                                                                                                                                        ggCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAg               <  1:3027085/150‑1 (MQ=255)
                                                                                                                                            cgTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTaa           >  1:1052701/1‑150 (MQ=255)
                                                                                                                                            cgTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTaa           <  1:1263252/150‑1 (MQ=255)
                                                                                                                                            cgTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTaa           <  1:3299491/150‑1 (MQ=255)
                                                                                                                                             gTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAAt          >  2:3442462/1‑150 (MQ=255)
                                                                                                                                             gTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAAt          <  1:528356/150‑1 (MQ=255)
                                                                                                                                             gTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAAt          >  1:3821630/1‑150 (MQ=255)
                                                                                                                                               aCAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATAc        <  1:2935282/150‑1 (MQ=255)
                                                                                                                                                 aGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCa      <  2:2597444/150‑1 (MQ=255)
                                                                                                                                                 aGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCa      <  1:64313/150‑1 (MQ=255)
                                                                                                                                                 aGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCa      <  2:578928/150‑1 (MQ=255)
                                                                                                                                                 aGTTCGTTTTCAATCTGCAACACGGTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCa      <  2:1777558/150‑1 (MQ=255)
                                                                                                                                                  gTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCaa     >  2:319213/1‑150 (MQ=255)
                                                                                                                                                     cGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGAt  <  1:2545186/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACCTGCTGTTCATCTACACCAATCGGCGAGAACGGGTTGATGTCCAGCGAACGCACTTCAATATATTCAATGCCGCCACGTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGAT  >  CP009273/2808758‑2809056

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: